Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913650,4580AG47.7% 114.1 / 24.0 90T11T (ACT→ACC) citEcitrate lyase, citryl‑ACP lyase (beta) subunit
Reads supporting (aligned to +/- strand):  ref base A (7/39);  new base G (42/0);  total (51/39)
Fisher's exact test for biased strand distribution p-value = 3.93e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCG  >  NC_000913/650320‑650609
                                                                                                                                            |                                                                                                                                                          
gcggcgGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGC‑GGGgtttttggttt                                                                                                                                                >  1:1480822/1‑144 (MQ=255)
      ggcggTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGG‑GAGGTTTTCGttgttggcg                                                                                                                                           >  1:384029/1‑147 (MQ=255)
           tgtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAg                                                                                                                                     <  2:892044/150‑1 (MQ=255)
           tgtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGG‑GGGCGTTTCGTTGTTGCAGC‑GAAg                                                                                                                                     >  1:256207/1‑150 (MQ=255)
              cTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGCGGGGTTTTAC‑TTGGTGCCGC‑GACGccgg                                                                                                                                 >  1:845743/1‑148 (MQ=255)
                ttttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGG‑GGGGTTTTCGTTGTTGCAgggcaacgggaa                                                                                                                                >  1:942810/1‑140 (MQ=255)
                ttttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGG‑GGGGTTTTCGTTTTTTCAGC‑GGAAC‑GGgaa                                                                                                                              >  2:1616934/1‑151 (MQ=255)
                   tCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCa                                                                                                                           <  2:707021/151‑1 (MQ=255)
                      cGCAATGCTACGGAGTCTTCGGGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTg                                                                                                                        <  2:838731/151‑1 (MQ=255)
                        cattGCTACGGAGTCTTCGAGGTCAAACAGCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGc                                                                                                                      <  1:525850/148‑1 (MQ=255)
                         aaTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCa                                                                                                                     <  2:398721/151‑1 (MQ=255)
                          aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGTTGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGG‑GGGTTTTTCTGTGTTGCGGC‑GCAGC‑GGCAAACAATTGCaa                                                                                                                    >  2:362705/1‑151 (MQ=255)
                          aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCGGGCACAAACAACATGCTGCGGC‑GGGGGG‑GGGTTTTTCGTTTTTGCCGC‑GAAAC‑GGGAAACATTTGcca                                                                                                                    >  2:1018849/1‑149 (MQ=255)
                          aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGG‑GGGTTTTTGGTTGTT‑CAGC‑GAAAC‑GGAGATCAtttggaa                                                                                                                    >  1:1153955/1‑146 (MQ=255)
                          aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGG‑GGGGTTTTTGGTTTTGGAGA‑GAAAC‑GGGAAACAtttggaa                                                                                                                    >  2:1518903/1‑146 (MQ=255)
                          aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGG‑GGGGTTTTCCTTGTTGCAGC‑GAAAC‑GGGAAACAtttggca                                                                                                                    >  1:339108/1‑146 (MQ=255)
                          aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGG‑GGGGTTTACGGTGTTGGAGG‑GAAGA‑GGGAAACaattggaa                                                                                                                    >  2:203127/1‑144 (MQ=255)
                          aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GGGGTTTTCGGTGTTGCAGC‑GAAGC‑GGAAAACATTGGGaa                                                                                                                    >  2:1665169/1‑151 (MQ=255)
                          aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGG‑GGGTTTTTCGTTGTTTCCGC‑GAAGA‑GGAAAACATTTGcga                                                                                                                    >  2:1573113/1‑149 (MQ=255)
                          aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGG‑GGGTTTTTCGTTGTTGTAGC‑GAAGC‑GGAAAACATTGGCaa                                                                                                                    >  2:1048732/1‑151 (MQ=255)
                          aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGG‑GGGGTTTTCGTTTTTGCAGC‑GAAGA‑GGGAAAAATTTGCaa                                                                                                                    >  2:596092/1‑151 (MQ=255)
                          aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGG‑GGGGTTTTCGTTGTTGCAGC‑GAAAC‑GGAAAACAtttggaa                                                                                                                    >  2:1487048/1‑146 (MQ=255)
                          aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGG‑GGGGTTTTCGGTTTTGGAGC‑GAAGC‑GGGAAACAtttggaa                                                                                                                    >  2:1228160/1‑146 (MQ=255)
                          aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGG‑GGGGTTTTCGGTGTTGGAGC‑GAAAA‑GGGAATAATTTGCaa                                                                                                                    >  1:186254/1‑151 (MQ=255)
                          aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGC‑GGGGTTTTTGTTGTTGGCGC‑GGAAC‑GGAAAACAATGGCaa                                                                                                                    >  2:616240/1‑151 (MQ=255)
                            gCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGG‑GGGGTTTTCGTTGTTTCCGC‑GAAAC‑GGAAAACATTTGccaa                                                                                                                   >  2:1253149/1‑147 (MQ=255)
                            gCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GGGGTTTTCGTTGTTGCTGC‑GAAAC‑GGAAATCATTGGCCAt                                                                                                                   >  1:340725/1‑150 (MQ=255)
                            gCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGC‑GGGGTTTTCGTTTTTGCAGC‑GAAAC‑GGGAAACATTGGCaaa                                                                                                                   >  2:292123/1‑149 (MQ=255)
                            acTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGG‑GGGGTTTTCGTTGTTGCAGC‑GAAGC‑GGAAAAAATTGGGAAt                                                                                                                   >  1:346035/2‑150 (MQ=255)
                                    gTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GGGGTTTTCGTTTTTGCAGC‑GGAGC‑GGGAAACATTGGCAAAccc‑ccccag                                                                                                         >  2:1150339/1‑148 (MQ=255)
                                    gTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGG‑GGGGTTTTTGTTGTTTCAGC‑GAAGC‑GGGAAAAATTGGGAATCCT‑CCCATg                                                                                                         >  1:498155/1‑151 (MQ=255)
                                     tCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATg                                                                                                         <  1:1370385/150‑1 (MQ=255)
                                      cTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGG‑GGGGTTTTCGTTGTTGCAGC‑GAAGC‑GGGAATCATTTGCAAaccccgtctctt                                                                                                        >  1:968880/1‑142 (MQ=255)
                                      cTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAAtcc                                                                                                                 <  2:969397/142‑1 (MQ=255)
                                      cTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGc                                                                                                       <  2:731195/151‑1 (MQ=255)
                                      cTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGc                                                                                                       <  2:1324773/151‑1 (MQ=255)
                                      cTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGc                                                                                                       <  2:1541866/151‑1 (MQ=255)
                                      cTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGc                                                                                                       <  1:129038/151‑1 (MQ=255)
                                      cTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGc                                                                                                       <  2:427588/151‑1 (MQ=255)
                                      cTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGc                                                                                                       <  1:447316/151‑1 (MQ=255)
                                      cTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGc                                                                                                       <  2:1768442/151‑1 (MQ=255)
                                      cTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGc                                                                                                       <  1:1071819/151‑1 (MQ=255)
                                            ggTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGC‑GGGGTTTTTGTTTTTGCCGC‑GAAAC‑GGAAAACATTTGCAAACCC‑CCCCTGGGAGAGc                                                                                                  >  2:289972/1‑150 (MQ=255)
                                                  aCATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GGGGTTTTTGGTTTTGCAGC‑GCAAC‑GGGAAACATTGGCCATCCT‑CCCCTTGGAGAGCCcgggga                                                                                           >  1:38756/1‑148 (MQ=255)
                                                   cATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGG‑GCGGGTTTTGTTGTTGGAGG‑CAACC‑GGAAAAAAGTGGGAAGCCT‑CCCCTAGGagaggtgggcga                                                                                           >  1:809106/1‑142 (MQ=255)
                                                     tCAGGGCATCAGCCGGGGAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGCTGA‑GCGGTTTGTGGTGTTGGTGCGGAAGC‑‑GAAACCATTGGCCAACCTGCCCATGGCGAAGCTgggtgac                                                                                          >  1:1201118/1‑146 (MQ=255)
                                                        gggCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGG‑GGGTTTTTTTTTTTTTTTGC‑GAAAC‑GGAAATTCTTTGCAATCCC‑CCCCCGGGCGGGCCGGGGTAccc                                                                                        >  2:1479963/1‑147 (MQ=255)
                                                        gggCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGC‑GGGGTTTTAGTTGTTGTAGC‑GAAAA‑GGGAAACATTTGGAATACC‑CCCCTGGCCGGGGCgggggacccccg                                                                                     >  2:1744564/1‑142 (MQ=255)
                                                         ggCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGG‑GAGGTTTTCGTTGTTTCAGC‑GAAAC‑GGAAAACATTGGGAATCCT‑CCCCTGGGAGAGCCGGGGTAcccccgg                                                                                    >  2:1471171/1‑146 (MQ=255)
                                                         ggCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGC‑GAGGTTTTCGTTGCTGCTGC‑CAAGC‑GGAGAACCTTGTCAATCCT‑CCCCTGGCCCAGCCgcgggtaccccgg                                                                                    >  1:1527719/1‑139 (MQ=255)
                                                           cATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGgcg                                                                                  <  2:756244/151‑1 (MQ=255)
                                                                gCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGG‑GGGGTTTTCGTTGTTTCAGC‑GGAAC‑GGGAAACAATGGGAAACCC‑CCCCCCGCCGAGGCGGGGGACCGCTGGcgcggg                                                                              >  1:1747578/1‑148 (MQ=255)
                                                                gCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGG‑GGGGTTTTCGTTGTTGCAGC‑GAAAC‑GGGAAACATTTGCACTCCT‑CCCCTGGCAGAGCCCGGGGACCGCCGGcgcggg                                                                              >  1:1668571/1‑148 (MQ=255)
                                                                    ggTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGG‑GGGGTTTTCGTTTTTGCCGC‑GAAGC‑GGGAAAAATTTGCAAACCC‑CCCCTGGCAGAGCCGGGGGACCCCCGGCGCGGGcccg                                                                          >  2:1193495/1‑148 (MQ=255)
                                                                     gTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCCCCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagc                                                                         <  1:406506/150‑1 (MQ=255)
                                                                      tAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagcag                                                                       <  1:913516/151‑1 (MQ=255)
                                                                      tAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagcag                                                                       <  2:161056/151‑1 (MQ=255)
                                                                        gATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGG‑GGGGTTTTAGTTGTTGTGGC‑GAAGC‑GGAAAACATTGGCAAACCC‑CCCCTTGGAGAGGCGGGGTAACCCTGGcgcggggcaccaagg                                                                     >  1:1658901/1‑140 (MQ=255)
                                                                         aTGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAggg                                                                     <  1:4563/150‑1 (MQ=255)
                                                                         aTGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGC‑GGGTTTTTCGTTGTTGCAGC‑GCAGC‑GGAAAACATTGGCAATCCT‑CCCATGGCAGAGCCGGGGTACCGCCGGCGCGCGCCCGCCgggg                                                                    >  2:1013487/1‑150 (MQ=255)
                                                                               gAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGG‑GGGGTTTTCGTTTTTTCCGC‑GAAACGGGAAAAATTTGGCATCCCT‑CCC‑TGGCCGAGCCCGGGGACCCCGGGGGCGGGCCACCCGGGCTTCCa                                                               >  1:710996/1‑150 (MQ=255)
                                                                                aGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGC‑GGGGTTTTTGTTGTTTGCGC‑GGAAC‑GGGAAATATTTGGAATACC‑CCCCTGGGCGAGGCGGGAGACCCCCGGCGCGGGCCACCCGGGGCTCCAg                                                              >  1:936669/1‑150 (MQ=255)
                                                                                  ttGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCGTTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGt                                                             <  2:943323/149‑1 (MQ=255)
                                                                                  ttGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGt                                                             <  2:346152/149‑1 (MQ=255)
                                                                                  ttGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGGGG‑GAGGTTTTCGCTGTTTCAGC‑GAAAC‑CGGAAACCTTGGGAACCCT‑CCCCCCGCAGAAGCGGGGCATCCCGGGCGCCTGCCCGCAGGGGtttcccgcg                                                           >  2:260678/1‑144 (MQ=255)
                                                                                    gCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GGGGTTTTTGGTGTTTCGGC‑GAAAC‑GGAAATTATTTGCAATCCT‑CCCCCGGCAGAGCCGGGGGAACGCCGGCGCGGGCCCGCAGGGCCTcccggccgg                                                         >  1:1381095/1‑144 (MQ=255)
                                                                                                 gCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGc                                             <  2:614530/150‑1 (MQ=255)
                                                                                                 gCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGc                                             <  2:657567/150‑1 (MQ=255)
                                                                                                     tGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTcc                                        <  2:984061/151‑1 (MQ=255)
                                                                                                     tGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTcc                                        <  2:268695/151‑1 (MQ=255)
                                                                                                     tGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTcc                                        <  1:1065647/151‑1 (MQ=255)
                                                                                                     tGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGCTTTATGTTTTTGCACC‑GAAGC‑GGAAAATAATGTCAAACCC‑CCCATGGGCGGGGCGGGGAACCCCCGGGCCGTGCCAGCCGGGCTTCCCGTCGTGCCCGAaaaaaacagac                                         >  2:419661/1‑143 (MQ=255)
                                                                                                      ggCACCAGGCACAAACAACATGCTGCGGC‑GGGGGC‑GGGTTTTATGTTGTTGCAGC‑GAAGC‑GGAAATCATTTGCCATCCT‑CCCCTGGCCGAGCCGGGATACCCCCGGCGCGTGCCCGCAGGGGTTCCAGTCGTGCACGGAAAATGCAGTcc                                        >  1:1482327/1‑150 (MQ=255)
                                                                                                        cACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAg                                      <  2:960038/150‑1 (MQ=255)
                                                                                                          ccAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCG‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTgc                                   <  2:238982/151‑1 (MQ=255)
                                                                                                             ggCACAAACAACATGCTGCGGC‑GGGGGG‑GGGGTTTTCGTTTTTTCCGC‑GAAAC‑GGAAAACATTTGCAATCCC‑CCCCTGGCAGAGCCGGGGTAACGCCGGGGCGGGGCCGCAGGGCTTCCCGGCCGGCCCGTAAAAAGCAATCCAgggggc                                 >  2:989361/1‑145 (MQ=255)
                                                                                                             ggCACAAACAACATGCTGCGGC‑GGGGGG‑GGGGTTTTAGTTGTTGCAGG‑GGAAC‑GGAAAAAATTTGCAAACCT‑CCCCTGGCCAAGGCGGGGGAGCGCGGGGGCGGGGCAACAGGGGTTCCCCTCGGGCGCGGAAAAAACAGGCCAGAgcgg                                 >  2:694369/1‑149 (MQ=255)
                                                                                                             ggCACAAACAACATGCTGCGGC‑GGGGGG‑GGGGTTTACGTTGTTGGAGC‑GAAAC‑GGAAAACATTGGGAAAACC‑CCCCAGGCAGAGCCCGGATAACGCCGGCGCGTGCCAGCAGGGCTTCCAGTCCGGCACGGAAAAATCAGGCCAGTgcgc                                 >  1:1363060/1‑150 (MQ=255)
                                                                                                                   aacaacATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTc                          <  2:1426073/151‑1 (MQ=255)
                                                                                                                    acaacaTGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCa                         <  2:38770/151‑1 (MQ=255)
                                                                                                                    acaacaTGCTGCGGC‑GGGTGC‑GAGGTTTTCGGTGTTTTAGC‑GAAAA‑GGAAAACATTGGGAATCCT‑CCCACGGGAGAGGCGGGGTACCGCTGGCGCGTGGCCGCAGGGGTTCCAGTCGTGCACGGAAAAAGCAGTGCAGTGGGCCTTTGTCa                         >  1:1622304/1‑151 (MQ=255)
                                                                                                                                 gc‑ggGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGc             <  1:1347202/150‑1 (MQ=255)
                                                                                                                                 gc‑ggGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGc             <  1:682280/150‑1 (MQ=255)
                                                                                                                                      nngN‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAc        <  2:503890/149‑1 (MQ=255)
                                                                                                                                       tGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCGTGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAc        <  1:1404208/150‑1 (MQ=255)
                                                                                                                                       tGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAc        <  1:285319/150‑1 (MQ=255)
                                                                                                                                         c‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAcgc      <  1:624100/150‑1 (MQ=255)
                                                                                                                                         c‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAcgc      <  1:70747/150‑1 (MQ=255)
                                                                                                                                         c‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAcgc      <  2:1718204/150‑1 (MQ=255)
                                                                                                                                           gAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGGGCGCCTTTGTCATCGACATTCAGCTGTAcgccg    >  2:427508/1‑151 (MQ=255)
                                                                                                                                            aGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGCGCGCCTTTGTCATCGACATTTAGCTGTAcgccg    >  2:327401/1‑150 (MQ=255)
                                                                                                                                            aGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGC‑TCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCcgcg  >  1:1647523/1‑151 (MQ=255)
                                                                                                                                            |                                                                                                                                                          
GCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGC‑GGGTGC‑GAGTTTTACGTTGTTGCAGC‑GAAGC‑GGAAATCATTGGCAATCCT‑CCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCG  >  NC_000913/650320‑650609

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑