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

TTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAGGG  >  NC_000913/29347‑29646
                                                                                                                                                      |                                                                                                                                                      
ttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAatat                                                                                                                                                        >  1:1025867/1‑150 (MQ=255)
  aaTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATaa                                                                                                                                                      >  2:724582/1‑150 (MQ=255)
    tGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAGAATAATATaaaa                                                                                                                                                    >  1:986554/1‑150 (MQ=255)
     gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATaaaaa                                                                                                                                                   <  2:334908/150‑1 (MQ=255)
         aTTTTGACCATTTGGTCCACTTTTCTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATccc                                                                                                                                              >  1:377416/1‑151 (MQ=255)
          ttttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCg                                                                                                                                             >  1:700119/1‑151 (MQ=255)
            ttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGcc                                                                                                                                           >  1:887756/1‑151 (MQ=255)
             tGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGcc                                                                                                                                           <  2:993493/150‑1 (MQ=255)
                    ttGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGtt                                                                                                                                   >  1:782770/1‑151 (MQ=255)
                        tCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACt                                                                                                                               >  2:49643/1‑151 (MQ=255)
                         ccttttttttCTGCTCTTTTTTATTTCATGCAATCTTCTTGCTGCGCAACCTTTTCCAAAAAAAGTTTAAATAATTTTTTTGCCCTTAAATCCCTAAAAAAAATGAAGGAC‑CAAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGa                                                                                                                                 <  2:178401/144‑1 (MQ=255)
                         ccATTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCTTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGAAAAACATGAATGAG‑CCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACt                                                                                                                               <  1:56369/150‑1 (MQ=255)
                              ttttttTGCTCTTTTTTATTTAATAAATTCTTTTTGCTGCGCAAGCTTTTCCAAAAACAGTTTAAATAATTTTTTTGTCTTTTAATGCCTAAAAAACATGCATGAG‑CCAAAAAATAATAAAAAAACTCCCGCCATTAAGTTGACTTTTAg                                                                                                                          <  2:238249/150‑1 (MQ=255)
                              tttttCTGCTTTTTTTTTTTTAATGCAATTTTTTTGCTGCGCAAGCGTTTCCAAGAACAGGTTAAATTATTTTTTTTTCGTTAATTGCCTGAAAAACATGAATGAC‑CCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAg                                                                                                                          <  2:608052/150‑1 (MQ=255)
                                 ttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGcc                                                                                                                      >  2:985746/1‑151 (MQ=255)
                                  tCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGccc                                                                                                                     >  2:565538/1‑151 (MQ=255)
                                      cTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCatat                                                                                                                 >  2:365987/1‑151 (MQ=255)
                                        cGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtct                                                                                                               >  1:596854/1‑151 (MQ=255)
                                        cGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtct                                                                                                               >  2:877124/1‑151 (MQ=255)
                                         gTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                              >  2:671397/1‑151 (MQ=255)
                                          tttttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                             >  2:500409/1‑151 (MQ=255)
                                          tttttATCTCATGCAATCTTTTAGCTGCAAAAGGTTTTTCAAAAACAGGTGAGATAATTTTTTGGTCGCTATGTGCATAAAAACTATGCAGGAC‑CCAAAAAAAAATAAAAAAATGCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                              <  1:1073195/150‑1 (MQ=255)
                                                tttCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATg                                                                                                       >  1:660706/1‑151 (MQ=255)
                                                tttCATGAATTCTTCTTGCTGCGCAAGCTTTTCCCAAAACAGTTTAGATTTTTTTTTTGTCCTTTATTCCCTAAAAAACATGCATGAA‑CCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGaa                                                                                                         <  2:528006/149‑1 (MQ=255)
                                                   cATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccg                                                                                                    >  2:845652/1‑151 (MQ=255)
                                                        aaTCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGtt                                                                                               >  1:59393/1‑151 (MQ=255)
                                                            tttttGCGGCGAAGGCTTTTCCCAAAACAGGTTAGATGTTTTTTTTGTCGCTAAAGCCCTGAAAAACATGCATGAC‑CCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGc                                                                                            <  1:448968/147‑1 (MQ=255)
                                                            ttcttGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                           >  2:43979/1‑151 (MQ=255)
                                                               ttGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                           >  2:311700/1‑148 (MQ=255)
                                                                  ctgcGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAt                                                                                     >  2:137020/1‑151 (MQ=255)
                                                                   tgcGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAtt                                                                                    >  1:1058087/1‑151 (MQ=255)
                                                                    gcgcAAGCGTTTCCCAGAACAGTTTAGATGATTTTTTTGTCGCTTAATGCCTGAAAAACATGCATGAG‑CCACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtc                                                                                   <  1:817856/151‑1 (MQ=255)
                                                                     cgcAAGCTTTTTCAAGAACAGGTTAGATGATTTTTTTGTCGCTTAATGCCTTAAAAACATGCATGAC‑CCACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAtt                                                                                    <  2:602723/149‑1 (MQ=255)
                                                                        aaGCTTTTCCAAAAACAGTTTAGATGTTCTTTTTGTCTCTTAATCCTTGTAAACCATGCTTGAC‑CCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAtt                                                                                    <  2:812684/146‑1 (MQ=255)
                                                                        aaGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                               >  2:406479/1‑151 (MQ=255)
                                                                        aaGCGTTTTCAAGAACAGGTTAGATGATTTTTTTGTCGTTAAATGCCTGAAAAACATGCATGAC‑CCAAAAAAAAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGa                                                                                        <  1:839297/142‑1 (MQ=255)
                                                                              tttCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                         >  2:731112/1‑151 (MQ=255)
                                                                              tgcccAGAAAGGGTTATATAATTTTTTTTTCGTTAAACTCTTGAAAACCATCTAGAAC‑CCAAAAAATAATAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAg                                                                          <  1:216652/147‑1 (MQ=255)
                                                                                   aaaaCAGTTTAGATGATTTTTTTGTCTCTTAATGCCTAAAAAACATGCATGAC‑CCACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAt                                                                     <  1:594014/148‑1 (MQ=255)
                                                                                   aGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAtt                                                                    >  1:603253/1‑151 (MQ=255)
                                                                                     aaCAGTTTAGATTTTTTTTTTGTCGTTTAATGCTTAAAAAACATGCATGAC‑CAAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAttt                                                                   <  1:186869/150‑1 (MQ=255)
                                                                                       ctgtttAGTTGTTTTTTTTGTCTCTTAACGCCTAAAAAACATGCA‑GAGACCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGa                                                                      <  2:668843/141‑1 (MQ=255)
                                                                                           ttAAATGTTTTTTTTTTCGTTTAATGCCTAAAAAACATGCATGAC‑CCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAt                                                               <  2:335914/148‑1 (MQ=255)
                                                                                                  atCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGt                                                     >  1:100973/1‑151 (MQ=255)
                                                                                                   tttttttGTCGTTTAATGCCTAAAAAACATGCATGCC‑CCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAtt                                                              <  2:1078219/139‑1 (MQ=255)
                                                                                                   tCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                                    >  2:811992/1‑151 (MQ=255)
                                                                                                       tttGTCGCTTAATGCCTGTAAAACATGTATGAC‑CCCCAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatc                                                 <  1:351730/150‑1 (MQ=255)
                                                                                                         tGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                                              >  2:142490/1‑151 (MQ=255)
                                                                                                          gTCGCTTAATGCCTGTAAACCATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGAGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAACTTCGTCGGTAAGCAGATTTGCATTGATTTACGTcatcat                                               >  2:555470/1‑149 (MQ=255)
                                                                                                          gTCGCTTAATGCCTAAAAACCATAAACAAC‑CCAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                                              <  2:570680/150‑1 (MQ=255)
                                                                                                            cGCTTAATGCCTGAAAAACATGCATGAC‑CCACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgt                                            <  1:433528/150‑1 (MQ=255)
                                                                                                                      ctaaaaaaCATGAATGAG‑CCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAatat                                  <  1:28108/146‑1 (MQ=255)
                                                                                                                      cTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                                 >  1:921685/1‑151 (MQ=255)
                                                                                                                      cTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                                 >  1:266744/1‑151 (MQ=255)
                                                                                                                          aaaaCATGCATGAG‑CCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaaa                             <  2:86267/151‑1 (MQ=255)
                                                                                                                          aaaaCATGCATGAC‑CCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTaa                                     <  1:223824/143‑1 (MQ=255)
                                                                                                                          aaaaCATGCATGAC‑CCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaaa                             <  2:754699/151‑1 (MQ=255)
                                                                                                                                       accaaaaaaaaaATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagt                   <  2:161283/143‑1 (MQ=255)
                                                                                                                                           aaaaaaaaATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat               <  2:616078/147‑1 (MQ=255)
                                                                                                                                               aataatATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtctc         >  1:652234/1‑151 (MQ=255)
                                                                                                                                               aataatATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtctc         >  2:192329/1‑151 (MQ=255)
                                                                                                                                                  aataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg      >  1:32613/1‑151 (MQ=255)
                                                                                                                                                  aataAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTg       <  1:419721/150‑1 (MQ=255)
                                                                                                                                                    tataAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg      >  2:590036/1‑149 (MQ=255)
                                                                                                                                                      taAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAggg  >  1:740346/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
TTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAG‑CCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAGGG  >  NC_000913/29347‑29646

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