Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091329,4960TA43.3% 86.4 / 45.3 67intergenic (+301/‑155)dapB/carAdihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base T (31/7);  new base A (0/29);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.08e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGG  >  NC_000913/29349‑29642
                                                                                                                                                       |                                                                                                                                                  
aaTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATaaa                                                                                                                                                 <  2:132765/151‑1 (MQ=255)
aaTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATaaa                                                                                                                                                 <  1:191231/151‑1 (MQ=255)
      aaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATcc                                                                                                                                           <  1:182826/151‑1 (MQ=255)
      aaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATcc                                                                                                                                           <  1:377844/151‑1 (MQ=255)
          ttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGcc                                                                                                                                       >  1:23225/1‑151 (MQ=255)
          ttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGcc                                                                                                                                       >  2:3207/1‑151 (MQ=255)
          ttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGc                                                                                                                                        <  2:855897/150‑1 (MQ=255)
           tGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGACACAAAATAATATAAAAAATCCCGcc                                                                                                                                       >  2:473685/1‑150 (MQ=255)
             aCCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCAtt                                                                                                                                    >  1:615986/1‑151 (MQ=255)
                  ttGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGtt                                                                                                                               >  1:53458/1‑151 (MQ=255)
                  ttGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGtt                                                                                                                               >  2:20297/1‑151 (MQ=255)
                        cACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACttt                                                                                                                         >  2:853419/1‑151 (MQ=255)
                           ttttttttGCTTTTTTTTTTTGCATGCATTCTTCTTCCTGCGCAGGCTTTTCCA‑AA‑AACAGTTTAAATTATTTTTTTGTC‑TTTAATTCCTTATAAAA‑ATTGCATGACCAAAAAAAAAATAAAAAAATCCCCGCCATTAAGTTGACTTTTa                                                                                                                       <  1:611218/150‑1 (MQ=255)
                           tttttcTTGCTTGTTTTTATTTCATGAAATTTTCTTGCTGCGCAAGCTTTTCCA‑AA‑AACAGTTAAAATAATTTTTTTTTCGCTAAATCCC‑TGAAAAA‑CATGCATGTGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACtttt                                                                                                                        <  1:171103/149‑1 (MQ=255)
                              tttCTCCTTTTTTTTTTTTCATGAACTTTTTTTGCGGCGAAAGCTTTTCCC‑AA‑AACGGGTAAGAGTTTTTTTTTGTCTTTAATTGCT‑TAAAAAA‑CATGAATGACCAAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAgcg                                                                                                                    <  2:877700/150‑1 (MQ=255)
                                  tGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                               >  2:406690/1‑151 (MQ=255)
                                        ttttttttTCATGCAATCTTCTTGCTGCGCAGGCTTTTCCA‑AA‑AAAAGTTTAGATGTTTTTTTTGTCGTTTAATGCC‑TAAAAAA‑CATGCATGAGCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCata                                                                                                              <  2:720943/146‑1 (MQ=255)
                                        ttttttttTCAGCAAATTTTCTTGCGCCGAAAGGGTTTCCA‑AG‑AACAGTTAAGATAATCTTTTTGCCGTTAAATGCT‑TGAAAAA‑CATGCAGGGCCCAAAAAAAAAAATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                          <  2:319612/150‑1 (MQ=255)
                                         ttttttttCATGCATTCTTCTTGCTGCGCAAGCTTTTTCA‑AG‑AACAGGTTAGATGTTTTTTTTGTCGTTTAATGCC‑TGAAAAA‑CATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                         <  2:775491/150‑1 (MQ=255)
                                             aTTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAAt                                                                                                    >  1:137913/1‑151 (MQ=255)
                                                  aTGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgc                                                                                               >  1:714851/1‑151 (MQ=255)
                                                     cAATCTTCTTGCTGCGCAAGCGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGt                                                                                            >  1:277380/1‑151 (MQ=255)
                                                          ttcttGCTGCGCAAGCGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                       >  1:252632/1‑151 (MQ=255)
                                                                      aaGCTTTTACC‑AG‑AACAGGTTAGATGTTTTTTTTGTCGTTTAATCCC‑TGAAAAA‑CATGCATGACCCACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                           <  2:396669/151‑1 (MQ=255)
                                                                       aagtttttCC‑AAAAACAGTTTAGTTGTTCTTTTTGCCGTTTAATCCT‑CATAAAA‑TTCCCATGCCCCAAAAAACAATAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                           <  1:390024/147‑1 (MQ=255)
                                                                         cGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGt                                                                         >  1:8998/1‑150 (MQ=255)
                                                                           ttttCA‑AA‑AAAAGGTTAGATGTTTTTTTTGCCGTTTAATGCC‑TGAAAAA‑CATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCgg                                                                          <  2:155253/147‑1 (MQ=255)
                                                                            tttCCAAG‑AACAGTTTAAATGATTTTTTTGTCGTTTATGGCT‑TAAAAAA‑CATGCATGACCAAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgt                                                                             <  1:242950/144‑1 (MQ=255)
                                                                             tttccaaaaacggtttAGATGTTTTTTTTTTCTTTTATTCCT‑TAAAAAA‑CTTGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTaa                                                                       <  1:760753/136‑1 (MQ=255)
                                                                             ttCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                     >  1:451926/1‑150 (MQ=255)
                                                                              tCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                     >  2:83893/1‑149 (MQ=255)
                                                                                c‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                     >  1:616487/1‑147 (MQ=255)
                                                                                c‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                     >  1:226956/1‑147 (MQ=255)
                                                                                  aG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAtt                                                                >  2:860361/1‑151 (MQ=255)
                                                                                         ggTTAGATGATTTTTTTGTCGTTTAATGCT‑TGAAAAA‑CATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAt                                                           <  2:8855/150‑1 (MQ=255)
                                                                                           ttAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGTATTGCAtt                                                          >  1:639392/1‑149 (MQ=255)
                                                                                                attttttttttTCGCTTAAAACC‑TTAAAAACCAGCCAGAAGCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTg                                                         <  1:720236/141‑1 (MQ=255)
                                                                                                   tCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                                >  1:277762/1‑151 (MQ=255)
                                                                                                            cgtttAATGCT‑TAAAAAA‑CTTGCATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                                          <  1:435918/145‑1 (MQ=255)
                                                                                                             cgtttaatccttaaaaaa‑aaTGCATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgtg                                       <  2:794357/137‑1 (MQ=255)
                                                                                                                ttattcccttaaaaa‑CATGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGaa                                     <  1:244462/139‑1 (MQ=255)
                                                                                                                  aTGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTa                                  >  1:820998/1‑150 (MQ=255)
                                                                                                                      cctaaaaaa‑CAGGCATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAatat                              <  1:293740/146‑1 (MQ=255)
                                                                                                                      c‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                             >  2:835258/1‑151 (MQ=255)
                                                                                                                        tGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGc                            >  1:18808/1‑151 (MQ=255)
                                                                                                                        tGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGc                            >  1:181949/1‑151 (MQ=255)
                                                                                                                         gTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGc                            >  2:404374/1‑150 (MQ=255)
                                                                                                                         gTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGc                            >  1:419127/1‑150 (MQ=255)
                                                                                                                         gTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGc                            >  2:545298/1‑150 (MQ=255)
                                                                                                                           aaaa‑aaTGCATGAGCAAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCa                           <  2:697142/149‑1 (MQ=255)
                                                                                                                           aaaa‑CATGCATGCCCAAAAAAAAAATAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCa                           <  1:828742/149‑1 (MQ=255)
                                                                                                                           aaaa‑CATGCAGGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                             <  2:621468/147‑1 (MQ=255)
                                                                                                                             aa‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                        >  2:822807/1‑150 (MQ=255)
                                                                                                                                   gaatgaccaaaaaaaTAATAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                  <  1:880724/140‑1 (MQ=255)
                                                                                                                                     cgggcacaaaaaaTAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtga                 <  1:3206/141‑1 (MQ=255)
                                                                                                                                     atgAGCCACAAAATAATATAAAAAATCTCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGCTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                  >  1:749342/1‑149 (MQ=255)
                                                                                                                                      tgAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAag                    <  2:419178/146‑1 (MQ=255)
                                                                                                                                       gaccaaaaaaaaaATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaaa                     <  1:98707/137‑1 (MQ=255)
                                                                                                                                            aaaaaaaaATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGaa            <  2:100433/146‑1 (MQ=255)
                                                                                                                                            aaaaaaTAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat         <  2:740836/149‑1 (MQ=255)
                                                                                                                                             caAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAtt        >  1:46995/1‑151 (MQ=255)
                                                                                                                                             aaaaaaaATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat           <  1:835149/147‑1 (MQ=255)
                                                                                                                                              aaaaTAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat         <  1:178132/149‑1 (MQ=255)
                                                                                                                                                 ataataTAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTctct    <  2:451980/151‑1 (MQ=255)
                                                                                                                                                  taataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTg   >  2:227495/1‑151 (MQ=255)
                                                                                                                                                   aataAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTctct    <  2:23224/149‑1 (MQ=255)
                                                                                                                                                   aataAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg  <  1:692209/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
AATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCC‑TGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGG  >  NC_000913/29349‑29642

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