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

TAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGG  >  NC_000913/29348‑29642
                                                                                                                                                       |                                                                                                                                                  
tAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAAtata                                                                                                                                                   >  1:703269/1‑150 (MQ=255)
 aaTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATaaa                                                                                                                                                 <  2:1335643/151‑1 (MQ=255)
 aaTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATaaa                                                                                                                                                 <  2:208052/151‑1 (MQ=255)
  aTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATaaaa                                                                                                                                                <  2:999748/151‑1 (MQ=255)
      aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATAAAAAAAATc                                                                                                                                            <  2:95637/151‑1 (MQ=255)
      aaaTTTTGAACATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGATTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATaaaaa                                                                                                                                               <  2:816842/148‑1 (MQ=255)
        aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATccc                                                                                                                                          <  2:837067/151‑1 (MQ=255)
           ttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGcc                                                                                                                                       >  2:751171/1‑151 (MQ=255)
             gACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCa                                                                                                                                      <  1:1008365/150‑1 (MQ=255)
               ccATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCAtt                                                                                                                                    >  2:614765/1‑150 (MQ=255)
                 aTTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                                   >  2:82815/1‑149 (MQ=255)
                 aTTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                                   >  1:560113/1‑149 (MQ=255)
                 aTTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                                   >  1:914293/1‑149 (MQ=255)
                        ccttttttttCTTCTTTTTTTTTTTTTATCCATTTTTCTTGCTGCGCAAGCTTTTTCAAAAACAGT‑TAAAATGATTTTTTTGT‑CGTTTAATCCCTAAAAAACA‑TGCAGGGCCCAAAAAATAAAATAAAAAATCCCGCCATTAAGTTg                                                                                                                              <  1:1094961/143‑1 (MQ=255)
                         cttttttttCTGCTTTTTTTTTTTTCATGCATTCTTCTTGCTGCGAAAGCTTTTTCCAGAACAGG‑TTAGATGATTTTTTTGT‑CGTTTATTGCCTGAAAAACA‑TGCAGGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACt                                                                                                                           <  2:186513/146‑1 (MQ=255)
                           cTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                                        >  2:791879/1‑150 (MQ=255)
                            ttttttCTGCTCTTTTTTATTAAACGAAATCTTCTTCCTGCGCAAGCTTTTCCCAAAAAAGT‑TTAAATTATTTTTTTGC‑CGTTAAACGCCGGAAAAACA‑TGTATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTa                                                                                                                       <  1:1135088/150‑1 (MQ=255)
                            ttttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                                        >  2:1329533/1‑149 (MQ=255)
                             tttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgc                                                                                                                     >  1:695001/1‑151 (MQ=255)
                               ttgctgttttttttttttttATGCATTCTTCTTGCTGCGAAAGCTTTTTCAAGATGAGT‑TTAGATGTTTTTTTTTT‑CGTTTATTCCCTAAAAAACA‑TGCAGGAGCCAAAAAATAATAAAAAAATTCCCGCCATTAAGTTGACTTTTAgcg                                                                                                                    <  1:330479/140‑1 (MQ=255)
                                ttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGcc                                                                                                                  >  2:1171601/1‑151 (MQ=255)
                                 tCTGGTCGTTTATATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGGTTAATACCTGTAAAACA‑TTCATGCGCCACAAAATAATATAAAAAATCCCGCCATTAATTTGACTTTTAGCGccc                                                                                                                 >  1:86526/1‑151 (MQ=255)
                                           tttATTTCATGAATTTTTCTTGCTGCGCAAGCGTTTCCCAGAAAAGT‑TTAGATGATTTTTTTGT‑CGTTTAATCCCTAAAAAACA‑TGCATGATCAAAAAAAAAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCa                                                                                                        <  1:156497/150‑1 (MQ=255)
                                           tttATTTCATACATTTTTTTTGCGGCGCAAGCTTTTCCCAAAAAAGG‑TAAGATATTTTTTTTGT‑CGTTTAAGGCCTAAAAAACA‑TGCAGGACCAAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCa                                                                                                        <  1:942686/150‑1 (MQ=255)
                                            ttATTTCATGCAATCTTCTTGCTGCCCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTCT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                          >  1:1116326/1‑147 (MQ=255)
                                             tATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                                      >  1:445685/1‑150 (MQ=255)
                                              aTTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAAt                                                                                                    >  2:656629/1‑151 (MQ=255)
                                               tttCATGCAATTTTCTTCCTGCCCAAGCTTTTTCAAGAACAGG‑TTAGATTTTTTTTTTGT‑CGCTTAATGCTTGTAAAACA‑TGAAGGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAAt                                                                                                    <  2:1039739/150‑1 (MQ=255)
                                                ttCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgc                                                                                                  >  1:746486/1‑151 (MQ=255)
                                                   aTGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgc                                                                                               >  2:133709/1‑151 (MQ=255)
                                                        aTCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGttt                                                                                          >  2:1339767/1‑151 (MQ=255)
                                                         tcttttGGCTGCGAAAGCGTTTCCCAAAACAGT‑TAAAATAATTTTTTTGT‑CCCTTAATCCCTGTAAAAAA‑TGCATGAGCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGttt                                                                                          <  1:376999/150‑1 (MQ=255)
                                                         tcttctTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTg                                                                                         >  1:1224186/1‑151 (MQ=255)
                                                         tcttctTCCCCCGCAGGCTTTCTCCAAACCGGT‑TTAGATGTTTTTTTTTC‑CGTTTAAT‑CCTTAAAAACATTGCAGGACCAAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgc                                                                                               <  2:819687/145‑1 (MQ=255)
                                                            ttgcttggctaattcGTTTCCAATAACTGT‑TTAGAT‑ATTTTTTTGTGCGATAAGTGCCAGAAAAAAA‑GTCAAGAGCCAAAAAATAAAAAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                       <  1:884768/136‑1 (MQ=255)
                                                                gctgcGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGaa                                                                                   >  2:4758/1‑150 (MQ=255)
                                                                 ctgcGCAAGCTTTTCCCAGAACAGG‑TTAGATGATTTTTTTTT‑CGCTTAATGCCTGAAAAACA‑TGCATGAGCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGaaa                                                                                  <  2:1274373/150‑1 (MQ=255)
                                                                   gcgcAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtc                                                                               >  1:691576/1‑151 (MQ=255)
                                                                       aaGCTTTTTCAAGAACAGT‑TTAGATGATTTTTTTGT‑CGTTTAATGCTTGAAAAACA‑TGCAGGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcg                                                                              <  2:13614/148‑1 (MQ=255)
                                                                           gTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                                        >  1:947521/1‑150 (MQ=255)
                                                                            ttttCAAAAACAGT‑TAAGATGATTTTTTTGC‑CTTTTAATGCCTGAAAAACA‑TGCATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                                        <  1:664285/149‑1 (MQ=255)
                                                                             tttCCAAAAAAGGTTTAGATGATTTTTTTTT‑CTTTAAATGCGTGAAAAACT‑TGCAGGAGCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGt                                                                         <  2:1027514/148‑1 (MQ=255)
                                                                                 cAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                     >  2:1169657/1‑147 (MQ=255)
                                                                                    aaCAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTg                                                              >  2:587455/1‑151 (MQ=255)
                                                                                     aCAGT‑TTAGATTTTTTTTTTGT‑CTTTTATTGCCTAAAAAACA‑TGCAGGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGc                                                             <  1:1233643/151‑1 (MQ=255)
                                                                                            tAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAt                                                       >  2:382684/1‑151 (MQ=255)
                                                                                             aaattttttttttCT‑CGCTTACCGCCGAAAAACCA‑TGCTAGGCCCACAAAAAAAAAAAAAAATGCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAtt                                                      <  1:902054/143‑1 (MQ=255)
                                                                                              gatgatCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAttt                                                     >  1:342673/1‑151 (MQ=255)
                                                                                               atgatCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAttt                                                     >  2:1022313/1‑150 (MQ=255)
                                                                                                tgatCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                                   >  2:1139695/1‑151 (MQ=255)
                                                                                                  atCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGt                                                 >  1:1199879/1‑151 (MQ=255)
                                                                                                   tttttttGT‑TGTTTAATCCCTAAAAAAAA‑GTAATGACCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTa                                                    <  1:59182/145‑1 (MQ=255)
                                                                                                   tttttttGT‑CTTTTAATGCTTGAAAAACA‑TGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                                   <  2:1015443/146‑1 (MQ=255)
                                                                                                   tttttttGT‑CGTTTAATGCTTGTAAAACA‑TGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACg                                                  <  2:3833/147‑1 (MQ=255)
                                                                                                       tttGT‑CGTTTAATGCCTGAAAAACA‑TGCATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatc                                             <  2:332040/150‑1 (MQ=255)
                                                                                                        tttt‑CGTTTAATGCCTGAAAAACA‑TGTAGGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatca                                            <  1:1151923/147‑1 (MQ=255)
                                                                                                          gT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtg                                         >  1:339821/1‑151 (MQ=255)
                                                                                                             ctcTTAATGCCTAAAAAACA‑TGCATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtg                                         <  2:422362/147‑1 (MQ=255)
                                                                                                             cgtttAATGCTTGAAAAACA‑GCCATGGCCAAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtg                                         <  2:1179346/146‑1 (MQ=255)
                                                                                                                ttAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAt                                    >  2:772063/1‑151 (MQ=255)
                                                                                                                 tAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                                   >  1:9171/1‑151 (MQ=255)
                                                                                                                  aaTGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                                   >  1:1062664/1‑150 (MQ=255)
                                                                                                                  aaTCCCTGAAAAACA‑TGCTGGAGCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGaa                                     <  2:631519/148‑1 (MQ=255)
                                                                                                                    tGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTa                                  >  1:752167/1‑149 (MQ=255)
                                                                                                                      ccTGTAAAACA‑TGAATGCGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAata                               <  2:349243/150‑1 (MQ=255)
                                                                                                                       ctaaaaaaCA‑TGCATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAt                                    <  1:346080/140‑1 (MQ=255)
                                                                                                                       cTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCAGCCATTAAGTTGGCTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAata                               <  2:352500/149‑1 (MQ=255)
                                                                                                                          aaaaaCA‑TGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaa                          <  2:947400/150‑1 (MQ=255)
                                                                                                                           aaaaCC‑AGCATGACCCAAAAAATAACAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAatat                              <  2:1092843/146‑1 (MQ=255)
                                                                                                                           aaaaCA‑TGCATGAGCAAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                             <  2:305748/147‑1 (MQ=255)
                                                                                                                           aaaaCA‑TGCATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCa                           <  1:1009431/149‑1 (MQ=255)
                                                                                                                              aCA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa                       >  2:1223176/1‑150 (MQ=255)
                                                                                                                                aatgcgggaccaaaaaaaaaaaaaaaaaaaCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa                      <  1:1020548/136‑1 (MQ=255)
                                                                                                                                a‑tgcatgACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa                      <  2:359989/149‑1 (MQ=255)
                                                                                                                                     atgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                  >  2:1336804/1‑149 (MQ=255)
                                                                                                                                       gacccAAAAAATAATAAAAAAAACCCCGCCATTAAGTTAACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagt               <  2:499105/147‑1 (MQ=255)
                                                                                                                                       gAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtg              >  1:1066482/1‑151 (MQ=255)
                                                                                                                                        agcaaaaaaaaaaaaaaaaaaaTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagt               <  2:529591/143‑1 (MQ=255)
                                                                                                                                         gCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAgtgagtga             >  1:1060188/1‑150 (MQ=255)
                                                                                                                                         cccAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtga                 <  1:476552/145‑1 (MQ=255)
                                                                                                                                          ccAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat           <  2:493393/151‑1 (MQ=255)
                                                                                                                                            acaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat         >  1:1056420/1‑151 (MQ=255)
                                                                                                                                              aaaaaaATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGTTTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAtt        <  2:824871/150‑1 (MQ=255)
                                                                                                                                              aaaaTAATAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc       <  1:608316/151‑1 (MQ=255)
                                                                                                                                               aaaaaaTAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtct      <  2:860094/151‑1 (MQ=255)
                                                                                                                                                   aataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg  >  2:802532/1‑151 (MQ=255)
                                                                                                                                                   aataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTg   >  2:1291996/1‑150 (MQ=255)
                                                                                                                                                   aataTAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTg   <  2:220365/150‑1 (MQ=255)
                                                                                                                                                   aataAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg  <  2:884361/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
TAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGG‑TTAGATGATCTTTTTGT‑CGCTTAATGCCTGTAAAACA‑TGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGG  >  NC_000913/29348‑29642

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