Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,578,1430AC42.4% 127.0 / 12.7 66L79V (TTA→GTA) maeBmalic enzyme: putative oxidoreductase/phosphotransacetylase
Reads supporting (aligned to +/- strand):  ref base A (32/6);  new base C (0/28);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 4.05e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCC  >  NC_000913/2577999‑2578291
                                                                                                                                                      |                                                                                                                                                     
cGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCa                                                                                                                                                 <  2:858620/150‑1 (MQ=255)
 gAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGc                                                                                                                                               <  2:973557/151‑1 (MQ=255)
   gcgcGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCAc                                                                                                                                             <  2:14358/151‑1 (MQ=255)
      cggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCAccgc                                                                                                                                          >  1:6959/1‑151 (MQ=255)
        gcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ACCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCAccgccg                                                                                                                                        <  2:357509/151‑1 (MQ=255)
         cgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGt                                                                                                                                       >  1:249408/1‑151 (MQ=255)
            cAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGC‑GATTTTTTTAAACAGACCGCCCTTGCCTCCC‑ATCCCCGTTTT‑GCCT‑G‑CAAGCGCGCAATTG‑TTGCCTAC‑CCCCAGCACCGCCGt                                                                                                                                       <  2:214219/148‑1 (MQ=255)
               cTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGtt                                                                                                                                 >  1:708059/1‑151 (MQ=255)
                   aaTAAATTTGTCGGGGTCGATTTCGTCAACTTCAATGTAAAATCCATCAACCCCGGC‑AAATTTCTTAAAAAGACCGCCTTTGCCTTCC‑ATCACCGGTTT‑GCCG‑G‑CCAGCGCGCCAATG‑TGGCCAAC‑CCCCAGCACCGCCGTACCGTTagag                                                                                                                             <  1:795328/151‑1 (MQ=255)
                    aTAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTag                                                                                                                               >  1:546142/1‑148 (MQ=255)
                         ttggccggggTCGAGTTCGTAAACTTCAATGTCAAATACATCAACCCCGGC‑GAATTTCTTAAACAGACCTCCCTTCCCTTCC‑ACCCCCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAC‑CCCCAGCCCCGCCGTACCGTTAGAGATCAc                                                                                                                        <  2:749652/143‑1 (MQ=255)
                          ttGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCg                                                                                                                      >  2:767707/1‑151 (MQ=255)
                           tGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGc                                                                                                                     >  2:515096/1‑151 (MQ=255)
                                gggTCGAGTGCTACAAATAAAATATCAAATACACCAACCCCGGG‑AAATTTCTAAAAAAAACCGCCTTTTCCTTCC‑CCCCCGGCTTT‑CCCA‑C‑CGAGCCCCCAAATT‑TGCCCGAC‑CCCCAGAACCGCCGTACCGTTAGAGATCACCGccacca                                                                                                                <  2:100824/151‑1 (MQ=255)
                                         tCGTCAACTTCAATGTCAAATACATCAATCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTc                                                                                                       >  1:179665/1‑151 (MQ=255)
                                           gTCAACTTCAATGTCAAATACATCAATCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCgg                                                                                                     >  2:1050988/1‑151 (MQ=255)
                                                   cAATGTCAAATACATCAATCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtat                                                                                              >  1:29598/1‑150 (MQ=255)
                                                   cAATGTCAAATACATCAATCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtat                                                                                              >  1:244438/1‑150 (MQ=255)
                                                    aaTGTAAAATACATCAACCCGGGC‑AATTTTCTAAAAAAGAACGCCTTTCCCTTCC‑ATCCCCGTTTT‑GCCT‑C‑CCGCCGCCCCAATG‑TTGCCTAC‑CCCCGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGt                                                                                                <  2:602311/147‑1 (MQ=255)
                                                      tGTCAAATACATCAATCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATAttt                                                                                          >  1:866317/1‑151 (MQ=255)
                                                          aaaTCCATCAACCCGGGC‑GAATTTCTAAAACAGAACCCCTTTGCCTCCC‑ATCCCCGGTTT‑GCCT‑G‑CCAGCGCGCCATTT‑TTCCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTa                                                                                       <  1:279255/150‑1 (MQ=255)
                                                          aaaTCCACCAACCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTCCC‑ACCCCCGTTTG‑GCCT‑G‑CCAGCGCGCAAATG‑TTGCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTa                                                                                       <  1:78225/150‑1 (MQ=255)
                                                           aaatcATCAAACCCCGC‑CAATTTCTTAAAAAGGACCCCCCTTCCTTCCAACCCCCGTTTT‑CCCT‑C‑CAACCCCCCCATTG‑TCGCCTAC‑CCCCGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATAttt                                                                                          <  2:15224/143‑1 (MQ=255)
                                                               cATCAATCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCttt                                                                                 >  2:180871/1‑151 (MQ=255)
                                                               cATCAATCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCttt                                                                                 >  2:1004312/1‑151 (MQ=255)
                                                                 accaacccGGCGGATTTTCTTAAACAGAAC‑CCCCTGCCCTCC‑CCCCCCGGTTTCGCCG‑G‑CAAGCGCCCAAATG‑TT‑CCTAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCttt                                                                                 <  1:920141/145‑1 (MQ=255)
                                                                  caaacccggcaatttttCTAAAAAAGCACG‑CCTTTCCCTCC‑CTCCCCGGTTT‑GCCC‑GCCCGGCCCGCCAGTT‑TGCCCTAC‑CCCCGGGACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt                                                                                <  2:256322/136‑1 (MQ=255)
                                                                   aaccccGGC‑AATTTTCTAAAAAAAAACCCCCTTCCCTCCC‑TCCCCCGGTTT‑GCCG‑G‑CCGGCGCCCCAATG‑TTGCCTAC‑CCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAAc                                                                             <  1:931794/148‑1 (MQ=255)
                                                                    accccGGC‑AATTTTCTTAAACAGACCCCCTTTGCCTCCC‑ATCCCCGGTTT‑CCCT‑C‑CCAGCGCCCCAATT‑TTGCCTAC‑CCCCGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCttt                                                                                 <  1:1002797/144‑1 (MQ=255)
                                                                      cccGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACggg                                                                          >  2:531299/1‑151 (MQ=255)
                                                                       ccGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGt                                                                         >  2:987481/1‑151 (MQ=255)
                                                                          gC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtt                                                                      >  1:1032219/1‑151 (MQ=255)
                                                                             gAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtttt                                                                    >  1:81698/1‑151 (MQ=255)
                                                                             gAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtttt                                                                    >  1:936199/1‑151 (MQ=255)
                                                                                tttttaaaaCAGAACCCCCTTCCCTCCC‑ATCACCGTTTT‑GCCG‑C‑CCGGCCCCCCATTT‑TTCCCTAC‑CCCCGGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTc                                                                  <  2:408403/144‑1 (MQ=255)
                                                                                tttCTAAAACAGACCCCCCTTCCCTCCC‑A‑CCCCGGTTT‑GCCTGG‑CCGGCGCGCCGATT‑TTGCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttttt                                                                   <  2:749744/149‑1 (MQ=255)
                                                                                      aaaCAGAACCCCCTTGCCTCCA‑ATCCCCGTTTT‑CCCT‑C‑CCACCGCGCCAATT‑TTCCCTAC‑CCCCGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTc                                                            <  2:616473/150‑1 (MQ=255)
                                                                                        aCAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa                                                           >  1:609975/1‑149 (MQ=255)
                                                                                          aGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGAc                                                       >  1:906514/1‑151 (MQ=255)
                                                                                          aGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGAc                                                       >  2:764867/1‑151 (MQ=255)
                                                                                           gAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGAc                                                       >  1:283300/1‑150 (MQ=255)
                                                                                            accccccTTGCCTTCC‑ATCCCCGTTTT‑CCCT‑C‑CCAGCGCCCAAATT‑TGCCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACaa                                                     <  1:78206/147‑1 (MQ=255)
                                                                                              cccccTTCCCTCCC‑A‑CCCCGGTTTGGCCT‑G‑CCGGCGCGCAAATG‑TTGCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAgg                                                   <  2:226780/149‑1 (MQ=255)
                                                                                              cGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAgg                                                   >  2:122210/1‑151 (MQ=255)
                                                                                               gCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGt                                                  >  2:267888/1‑151 (MQ=255)
                                                                                                cctttGCCTTCC‑ATCCCCGGTTT‑GCCC‑G‑ACGGCCCGCCAATT‑TTGCCAAC‑CCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGt                                                  <  1:72499/147‑1 (MQ=255)
                                                                                                  cTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCg                                               >  2:738061/1‑151 (MQ=255)
                                                                                                      ccTCCC‑ACCCCCGGTTT‑GCCT‑G‑CCGGCGCGCCAATT‑TGGCCAAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCg                                               <  1:833/147‑1 (MQ=255)
                                                                                                       cTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAAc                                          >  1:794456/1‑151 (MQ=255)
                                                                                                         tCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCAACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGc                                        >  2:747010/1‑151 (MQ=255)
                                                                                                           cctcccccGTTTG‑GCCT‑G‑CGAGCGCGCCATTT‑TTCCCACC‑CCCGACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGc                                        <  2:509283/144‑1 (MQ=255)
                                                                                                           ccacccccGTTTT‑GCCT‑G‑CCGGCGCCCAAATT‑TTCCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAAc                                          <  1:39912/142‑1 (MQ=255)
                                                                                                              cccccGTTTT‑GCCG‑G‑CAAGCGCCCAATTG‑TTCCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCt                                      <  2:673236/146‑1 (MQ=255)
                                                                                                               cACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGt                                   >  1:798165/1‑151 (MQ=255)
                                                                                                                     ttt‑GCCT‑G‑CCAGCGCGCAATTG‑TTGCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTa                              <  1:126425/150‑1 (MQ=255)
                                                                                                                           ccgg‑ccGGCGCCCAAATG‑TTCCCTAC‑CCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAg                             <  2:873723/144‑1 (MQ=255)
                                                                                                                            cgtgccgccgcgcCGTTG‑TTGCCTAC‑CCCGGAAACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGt                               <  2:884396/137‑1 (MQ=255)
                                                                                                                             cctcccggcgcccAATGTTTCCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAgg                            <  2:640839/137‑1 (MQ=255)
                                                                                                                                 cAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGcc                    >  2:531826/1‑151 (MQ=255)
                                                                                                                                    cgcccCAATT‑TTCCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGcc                    <  1:609540/148‑1 (MQ=255)
                                                                                                                                     gcgcCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGcc                    >  2:448057/1‑147 (MQ=255)
                                                                                                                                      cgcCAATG‑TGCCCTAC‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGAt                <  1:473556/150‑1 (MQ=255)
                                                                                                                                            tG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCt          >  2:561869/1‑150 (MQ=255)
                                                                                                                                            tG‑TTCCCCAC‑CCCCAACCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcgc             <  2:197221/147‑1 (MQ=255)
                                                                                                                                               ttGCCTAC‑CCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGACAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCt          <  1:511499/148‑1 (MQ=255)
                                                                                                                                                tGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTgtgt      >  2:689665/1‑151 (MQ=255)
                                                                                                                                                      a‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGcc  >  2:833173/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGC‑GAATTTCTTAAACAGAACGCCCTTGCCTTCC‑ATCACCGGTTT‑GCCT‑G‑CCAGCGCGCCAATG‑TTGCCTAA‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCC  >  NC_000913/2577999‑2578291

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