Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,805,237 A→C 29.2% intergenic (+164/‑585) dsbG → / → CJ019_04513 Thiol:disulfide interchange protein DsbG/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,805,2370AC29.2% 71.8 / 7.4 98intergenic (+164/‑585)dsbG/CJ019_04513Thiol:disulfide interchange protein DsbG/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (8/46);  new base C (4/15);  total (29/69)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGAAGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACATTGGTGATTGG  >  ACB143/4805103‑4805368
                                                                                                                                          |                                                                                                                                         
gTGCAGCAGGCGACGTGGGGCAGGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGAg                                                                                                                                          <  1:395366/136‑1 (MQ=255)
   aaGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTT‑‑GGGGCgg                                                                                                                                    >  6:72325/1‑136 (MQ=255)
   aaGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTG‑‑GGGGCgg                                                                                                                                    >  5:220371/1‑136 (MQ=255)
   aaGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTG‑‑GGACCgg                                                                                                                                    >  6:243655/1‑136 (MQ=255)
    aGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCCTA‑‑GGCGCcgg                                                                                                                                   >  2:116256/1‑136 (MQ=255)
           gACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTACGGGCCGGGTTTAcc                                                                                                                             >  4:346898/1‑136 (MQ=255)
            aCGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑TGGGCGCGGGTTAccg                                                                                                                            >  1:14677/1‑135 (MQ=38)
            aCGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGCGCGGGGTTTAccc                                                                                                                            >  7:89914/1‑136 (MQ=255)
            aCGTGGGGCATGGCAAATTATACGCCGGTGTTCGCGGGTGACCCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGAGGTCGGCCGTACGGGGGGGGTTTGccc                                                                                                                            >  3:303070/1‑136 (MQ=255)
                ggggCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCCTA‑‑GGGGCGGGGTTTACCcccg                                                                                                                         >  1:384078/1‑135 (MQ=38)
                      tGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTT‑‑GGGGCGGGGGTTACCGCCGAAGGGGc                                                                                                                  >  8:304955/1‑136 (MQ=37)
                          aCCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGCGCCGGGTTTACCCCCGCATGGCCCgcc                                                                                                              >  6:162048/1‑134 (MQ=38)
                               cTACGCCGGTGTTCGCTGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTAGTGGTGGGGGTTTACCCCGGACAGGGGCgccgagg                                                                                                          >  4:234537/1‑130 (MQ=255)
                               cTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTTG‑‑GGGGCGGGTTTCCCGCCCAAgggccgggccaggg                                                                                                       >  6:189069/1‑132 (MQ=37)
                               cTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCCGCGCGGTCGGCCGTGG‑‑GGGGCGGGGTTGCCCCCCAAGGGCCCggccagcc                                                                                                         >  6:118411/1‑130 (MQ=37)
                                  cGCCGGTGTTCGCGGGGTAACCCGCTCCTACGGGGCTCGCGCAAACGGCTGGGTTTTGCGCCCGGCCCGTGAGCGCACGGCCGGCGCGGTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGc                                                                                                   <  5:293775/136‑1 (MQ=37)
                                          ttCGCGGGGGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTc                                                                                            <  7:379339/135‑1 (MQ=255)
                                              cgGGTTAACCCGCTCCTACGGGGCGCGCGCTAACGGCGGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTccacca                                                                                       <  7:195932/136‑1 (MQ=255)
                                                 gttaccccGCTCCTACGGGGCGCGGGATACCGGGTGGGTTTTGCGCCCGGCCCGTTATGGCACGGCGTGCGGGGTGGGCCTTT‑‑GGTGGGGTTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGcc                                                                                    <  7:73580/131‑1 (MQ=16)
                                                 gttaaCCCGCTCCTACGGGGCGCGGGCTAAGGGGTGGGTTTTGCGCCTGGCCCGTTATCGCACGGCCTGCGCGGTCGGCCGTT‑‑GGTGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGcc                                                                                    <  6:354757/133‑1 (MQ=21)
                                                             gtACGGGGCGCGCGCTAAGGGGGGGGTTTTGCGCCCGGCCCGTGATCGCAAGCCCTGCGCGGTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTgc                                                                        <  2:21900/135‑1 (MQ=21)
                                                             cTAGGGGGCTCGCGCTACCGGCTGGGTTTGGCGCCAGGCCCGTGAGCGCACGGCCGGCGCGGGCGGCCGTA‑‑GGCGGGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTgc                                                                        <  8:261215/136‑1 (MQ=21)
                                                                gggggCGCGCGCTACGGGCTGGGTTTTGCGCCTGCCCGGTTAGCGCACGGCGGGCGGGGTCGGCGGTA‑‑GGCGGGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTgccgc                                                                     <  3:355377/135‑1 (MQ=18)
                                                                  gggccggcgcTATCGGGTGGGGTTTGCGCCTGGCCCGTGATCGCATGGCCGGGGGGGCGGGCCTTA‑‑GTTGGGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTgccgccg                                                                   <  7:230265/130‑1 (MQ=17)
                                                                  gggCTCGGGCAAACGGCGGGGTTTTGCGCCAGGCCGGTCATCGCACGGCCTGCGCGGGCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTgccgccg                                                                   <  6:181743/136‑1 (MQ=18)
                                                                  gggCTCGCGCTACCGGCGGGGTTTTGCGCCAGGCCGGTGATCGCAAGGCCTGCGCGGGCGGCCGTA‑‑GGCGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTgccgccg                                                                   <  7:259190/136‑1 (MQ=21)
                                                                  gggCTCGCGCTACCGGCGGGGTTTTGCGCCAGGCCGGTGATCCCACGGCCTGCCGGGTCGGCCGTA‑‑GGAGGGGTTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTgccgccg                                                                   <  6:438/136‑1 (MQ=21)
                                                                  gggCGCGCGCTAACGGGTGGGTTTTGCGCCAGGCCCGTGATCGCACGGCCTGCGCGGTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTgccgccg                                                                   <  7:371438/136‑1 (MQ=255)
                                                                  gggCGCGCCCTAAGGGCGGGGTTTTGCGCCAGGCCCGGTAGCGCACGGCGGGCGGGGTGGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTgccgccg                                                                   <  8:43353/136‑1 (MQ=18)
                                                                   ggCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTG‑‑GGGGCCGGGGTTACCCCGCAAGGGGCGGGCCCGCGGGGCCAATTAACACCACCCGTTCCCCGGTgcccccc                                                                   >  5:111645/1‑135 (MQ=255)
                                                                    gcgcgcgCTAACGCCTGGGGTTGGCGCCAGGCCCGTGATCGCAAGGCCGGCGCGGGCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTgccgccgcc                                                                 <  7:128954/133‑1 (MQ=21)
                                                                    gCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGGGCGGGGTTAACCGCGCAAGGGCCGGGCCGGCGGGGCCAATTCAACCAGCCGGTCACCCGTgccgccccc                                                                 >  1:251729/1‑136 (MQ=255)
                                                                        gcgcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTTG‑‑GGGGGGGGTTTTCCCCGGAAGGGGCCGCCCGGGGGGGCCAATTTGCCCAGCCGGGGACCGGTggcgccggccccg                                                             >  5:80154/1‑130 (MQ=255)
                                                                        gcgcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTT‑‑GGGGCCGGGTTTACCGCCCAACGGGCGGGCCGGGGGGGCCAATTACCCCACCCGGGAACCCTTgccgcccccccag                                                             >  5:256827/1‑132 (MQ=255)
                                                                        gcgcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTT‑‑GGGGCCGGGTTTACCCCGCAATGGGCGGCCCAGGTGGCCAATTTACCCCACCCGTTCAACGGTgccgcccccccag                                                             >  5:297717/1‑132 (MQ=255)
                                                                        gcgcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGAGCGGGTTTACCCCCCCAAGGGGCGGACCTGCTGGCCCATTTACCCCCGCCGTTGCCCGTTGCCCCCCCCCCTg                                                             >  8:273346/1‑136 (MQ=21)
                                                                        gcgcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGAGCGGGGTTACCCCCGAAAGGGCCCGCCCAGCTGGCCCATTTACCCCAGCCGGTCCCCGGTTCCCCCGCCCCTg                                                             >  8:95058/1‑136 (MQ=255)
                                                                        gcgcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGAGCGGGGTTACCCCCCAAATGGGCGGCCCAGCGGGCCAATTTCCCCCCGCCGGTCCCCCTTgccgccgccccag                                                             >  7:307852/1‑132 (MQ=21)
                                                                          gcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTT‑‑TGAGCGGGGTTTACCCCGCAAGGGGCGGCCCATCGGGGCAATTTTCACCCGGCGGTGACCGGTggcgccgccccagtt                                                           >  7:220033/1‑130 (MQ=255)
                                                                          gcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTT‑‑GGGGCCGGGTTTACCCCGCAAGGGGCCGGCCGGGGGGGCCAATTTCCCCCGCCCGTGACCCGTTCCGCCGCCCCAGTa                                                           >  2:211403/1‑136 (MQ=255)
                                                                          gcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGCGCGGGGTTTACCCCGCAAGGGGCCGGCCCGCGGGGCCAATTGCCCCCGCCGGTGACCGGTgccgccgccccagtt                                                           >  4:304833/1‑130 (MQ=255)
                                                                          gcTAACGGCGGGGTTTTGCGCCAGGCCGGTTATCGCACGGCCGGCGCGGTCGGCGGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTa                                                           <  2:47952/136‑1 (MQ=18)
                                                                          gacagcGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTa                                                           <  2:322679/131‑1 (MQ=255)
                                                                           aaaaCGGCTGGGTTTTGCGCCAGGCCCGTGATCGCACGGCCGGCGGGGGCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTAt                                                          <  2:58528/134‑1 (MQ=21)
                                                                                     tgTTTTGCGCCAGGCCGTTTTGCGAAAGGCCTGCGCGGTGGGCCGTA‑‑GGTGGGGTTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTaa                                               <  6:85619/135‑1 (MQ=14)
                                                                                     ggTTTTGGGCCCGGCCCGTGATCCCAAGCCCTGGGCGGGCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTaa                                               <  7:40322/136‑1 (MQ=18)
                                                                                     ggTTTTGCGCCAGGCCGGTGAGCGCAAGGCCGGCGCGGTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTaa                                               <  4:287217/136‑1 (MQ=255)
                                                                                     ggTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTaa                                               <  8:89059/136‑1 (MQ=255)
                                                                                     ggTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTaa                                               <  8:312986/136‑1 (MQ=255)
                                                                                     ggTTTTGCGCCAGGCCCGTGATCGCAAGGCCGGCGGGGTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTaa                                               <  2:384078/136‑1 (MQ=255)
                                                                                     ggTTTTGCCCCCGGCAGTTGAGCCCAAGCCCTGCGCGGTGGGCCTTA‑‑GGAGCGGTTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTaa                                               <  5:219933/136‑1 (MQ=14)
                                                                                     ggTTTTCCGACCGGACCGTGATCGCACGGCCTGCGCGGTCGACCGTA‑‑GGGGCGGTTTTACCAGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCATCCGTTCACCGTTGACGCCGCCACTGTATTCGCGGTTaa                                               <  8:113411/136‑1 (MQ=17)
                                                                                     ggTTTGGCGCCAGGCCCGTGAGCGCAGGCCCTGCGGGGGCGGCCGTA‑‑GGAGCGGTTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTaa                                               <  6:93660/136‑1 (MQ=18)
                                                                                     ggTTTGGCCCCTGCCCCGTTATCCCACGGCCGGCGCGGTGGGCCTTA‑‑GGCGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTaa                                               <  2:307684/136‑1 (MQ=17)
                                                                                     ggTTTGCCGCCCGCCCGGTGACCGCACGGCCGGCGGGGTCGGCGGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTaa                                               <  6:314782/136‑1 (MQ=17)
                                                                                      gTTTTGCGCCTGGCCCGTGATCGCATGGCCTGCGCGGTCGGCCGTA‑‑GGAGGGGGTTTACCCGCGAAAGGGCGGGACGGGCTGGCAAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTaaa                                              <  7:188788/136‑1 (MQ=18)
                                                                                      gTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGAGCCGGTTTTCCCCCGAAAGGGCCCGGCCGGCGGGCCAATTTACACCCGCCCTTGACCCGTGCCGCCGCCACAGTATTAGCGGGGaa                                               >  8:312917/1‑135 (MQ=255)
                                                                                         ttGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTT‑‑GGGGCGGGGTTTACCCCGCAAGGGGCGGCCCCGGGGGGCCAATTGCCCCCGCCGGTGACCGGTTGCGCCGCCCCTGTTTTAGCGGGGTAAAcc                                            >  2:241283/1‑136 (MQ=255)
                                                                                           gccccAGGCCGTTTATCCCACGGCCGGCGGGGTCGCCTTTT‑‑GCCGCGGTTTAACCGGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCg                                          <  7:209060/132‑1 (MQ=14)
                                                                                             gcCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTa                                                           <  7:153932/117‑1 (MQ=255)
                                                                                             gcCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTa                                                           >  8:153932/1‑117 (MQ=255)
                                                                                               ccggCCGGTGATCTCACGGCCGGCGGGGCCGGCCGTA‑‑GGCGGGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCt                                     <  3:329533/134‑1 (MQ=17)
                                                                                               ccggCCCGTGATCGCAAGGCCGGCGCGGTCGGCGGTC‑‑GGAGCGGTTTAACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCt                                     <  2:372510/134‑1 (MQ=17)
                                                                                               cAGGCCCGTGATCGCCCGGCCTGCGCGGGCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCt                                     <  4:306195/136‑1 (MQ=255)
                                                                                               cAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTT‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCt                                     <  2:292270/136‑1 (MQ=255)
                                                                                               cAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCt                                     <  8:106918/136‑1 (MQ=255)
                                                                                               cAGGCCCGTGACCGCAAGCCCGCCGCGGGCGGCGTTC‑‑GCCGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGATCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCt                                     <  5:100578/136‑1 (MQ=17)
                                                                                               acgccccGTGACCGCACGCCCTGCGGGGCCGGCCTTA‑‑GGAGGGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCt                                     <  5:113202/132‑1 (MQ=17)
                                                                                                  gCCCGTGATCGCACGGCCTGCGCGGTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACg                                  <  8:195716/136‑1 (MQ=255)
                                                                                                    cggttaTCCCACGGCCTGCGCGTGCGGCGTTA‑‑GGATCGGTTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCTTTCACCGTTGCCTCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCt                                <  8:157855/131‑1 (MQ=14)
                                                                                                    ccGTGATCGCACGGCCGGCGCGGTCGGCCGTA‑‑GGCGGGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCt                                <  7:75202/136‑1 (MQ=255)
                                                                                                    ccGTCAGCGCACGGCCGGCGCGGTGGGCGGTA‑‑GCAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCt                                <  3:278435/136‑1 (MQ=17)
                                                                                                      gTGATCGCACGGCGGGCGGGGTCGGCCGTA‑‑GTAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGa                              <  4:267657/136‑1 (MQ=17)
                                                                                                         agcGCACGGCCGGCGCGGGCGCCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCg                           <  4:294835/134‑1 (MQ=14)
                                                                                                         aTCGCACGGCCGGCGGGGGCGGCGGTA‑‑GCCGGGGTTTAACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATc                            <  2:89742/135‑1 (MQ=17)
                                                                                                          tCGCATGGCGGGCGGGGGGGGCGGTT‑‑GGCGGGGTTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCg                           <  5:196628/135‑1 (MQ=14)
                                                                                                          gcGCACGGCCGGCGGGGGCGGCCTTT‑‑GGAGGGGGTTAACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCg                           <  6:122572/134‑1 (MQ=17)
                                                                                                           cGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGCGGGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCgg                          <  6:358066/135‑1 (MQ=255)
                                                                                                           cGCAAGGCCTGCGCGGGCGGCGGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCgg                          <  5:179880/135‑1 (MQ=255)
                                                                                                             catggCCGGCGGGGTGGGCCGTT‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCggggg                       <  2:6294/133‑1 (MQ=17)
                                                                                                             cacggCCGGCGGGGGCGGCGGAA‑‑GTTGGGGGTTACCCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCggggg                       <  4:41124/133‑1 (MQ=17)
                                                                                                             cAAGGCCTGCGCGGTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCggggg                       <  1:233371/136‑1 (MQ=255)
                                                                                                               tggCCTGCGGGGTCGGCCGTA‑‑GGTGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGGGTg                     <  2:245645/135‑1 (MQ=255)
                                                                                                               cgggcTGCGCGGTGGGCCGTA‑‑GGCGGGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGGGTg                     <  5:127484/132‑1 (MQ=255)
                                                                                                               aGGCCTGCGCGGTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGGGTg                     <  1:248870/136‑1 (MQ=255)
                                                                                                                ggCCTGCGCGGGCGGCCTTT‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGGGTgg                    <  6:161160/136‑1 (MQ=255)
                                                                                                                gcccTGCGGGGTCGGCCGTA‑‑GCAGGGGTTTAACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGGGTgg                    <  1:214868/134‑1 (MQ=17)
                                                                                                                     gggcgGTCGGCCGTA‑‑GGCGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGaca               <  3:246591/134‑1 (MQ=255)
                                                                                                                      cgcgGACGGCCGTA‑‑GTAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGaca               <  5:68428/135‑1 (MQ=255)
                                                                                                                        gggTGGGCGGTA‑‑GCCGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACAt            <  2:245430/135‑1 (MQ=255)
                                                                                                                        gggTCGGCGGAA‑‑GTCGGGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACAt            <  2:75684/135‑1 (MQ=17)
                                                                                                                        gggTCGGCCGTA‑‑GGAGCGGTTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACAt            <  2:111005/135‑1 (MQ=255)
                                                                                                                        cgGTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACAt            <  8:33648/136‑1 (MQ=255)
                                                                                                                         ggTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTAt                                                          <  8:374611/90‑1 (MQ=255)
                                                                                                                         ggTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTAt                                                          >  7:374611/1‑90 (MQ=255)
                                                                                                                            gggCGGTA‑‑GGTGCGGTTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGTGTGGAGACACATTggg        <  2:100838/135‑2 (MQ=17)
                                                                                                                                    gtttggAGGGGTTTAACCCGGGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACATTGGTGAtt    <  2:396059/132‑1 (MQ=17)
                                                                                                                                       tggAGGGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACATTGGTGATTgg  <  8:4488/134‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
GTGAAGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTA‑‑GGAGCGGGTTTACCCGCGAAAGGGCCGGACCTGCTGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGCGGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACATTGGTGATTGG  >  ACB143/4805103‑4805368

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: