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

AATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCT  >  NC_000913/29349‑29640
                                                                                                                                                       |                                                                                                                                                
aaTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaa                                                                                                                                                <  2:1151852/150‑1 (MQ=255)
aaTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaa                                                                                                                                                <  1:631527/150‑1 (MQ=255)
 aTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaaa                                                                                                                                              <  2:842026/151‑1 (MQ=255)
    tAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaaaaa                                                                                                                                            <  2:168338/150‑1 (MQ=255)
     aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAt                                                                                                                                           <  1:284140/150‑1 (MQ=255)
     aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATc                                                                                                                                          <  2:832907/151‑1 (MQ=255)
      aaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATcc                                                                                                                                         <  2:675465/151‑1 (MQ=255)
         tttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCg                                                                                                                                       <  1:964078/150‑1 (MQ=255)
           tGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGcc                                                                                                                                     <  2:19129/150‑1 (MQ=255)
                   tGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTg                                                                                                                            >  1:623587/1‑151 (MQ=255)
                    ggTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTg                                                                                                                            >  1:300923/1‑150 (MQ=255)
                                 cTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGccc                                                                                                               >  1:869241/1‑150 (MQ=255)
                                      cGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtct                                                                                                         >  2:394919/1‑151 (MQ=255)
                                       ttttttATTTCATGCAATCTTTTTGCTGCGCAAGCTTTTCCAAAAACGGTTTAGATGATTTTTTTGTCGTTTAATGCTTGAAAAAAATGCAGGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGcc                                                                                                                <  1:241977/142‑1 (MQ=255)
                                        ttttttttTCATGCATTCTTCTTGCTGCGAAAGCTTTTCCAAAAAAAGTTAAGATAATTTTTTTGTCTTTAATTGCTTGAAAAACATGAATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                        <  1:457473/150‑1 (MQ=255)
                                        tttttATTTCATTCATTTTTCTTGCTGCGCAAGCGTTTTCCAAAACAGGTAAGATTATTTTTTTGTCTCTTAATGCCTGAAAAACATGCATGAGCCAAAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCata                                                                                                            <  2:353239/146‑1 (MQ=255)
                                        tttttATTTCATGCATTCTTCTTGCTGCGCAAGCTTTTTCAAGAACGGTTTAGATGATTTTTTTGCCGCTTAATGCCTAAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCatat                                                                                                           <  1:569715/147‑1 (MQ=255)
                                        tttttATTTCATGCAATCTTCTTGCTGCGAAAGCTTTTTCAAAAAAAGGTTAGATGATTTTTTTGTCGCTAAAGGCCTGAAAAACATGCTTGGCCCACAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                        <  2:371466/150‑1 (MQ=255)
                                        tttttATTTCATGAATTCTTCTTGCTGCGAAAGCGTTTTCAAGAACAGTTTAGATGATTTTTTTGTCGTTTAATGCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCatat                                                                                                           <  1:443476/147‑1 (MQ=255)
                                         ttttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                       >  1:1153747/1‑150 (MQ=255)
                                            tATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                                    >  2:1138008/1‑150 (MQ=255)
                                               ttCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgc                                                                                                >  2:378050/1‑151 (MQ=255)
                                                           tcttGTTGCGCAAGCGTTTCCCAAAACAGTTTAGATGATTTTTTTGTCGTTTATTCCCTTAAAAACATGCATGAGCCAAAAAATAATATAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTg                                                                                       <  1:959123/148‑1 (MQ=255)
                                                            cttGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                     >  2:889360/1‑149 (MQ=255)
                                                            cttGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                     >  2:436835/1‑149 (MQ=255)
                                                             ttGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                     >  2:734818/1‑148 (MQ=255)
                                                              ttccggcgaaggcttttCCAAGAAGAGTTTAGATGATTTTTTTGTCGCTTAATCCTTGAAAAACTTCCATGCCCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGttt                                                                                        <  2:936741/132‑1 (MQ=255)
                                                              tgctgcGCAGTCTTTTTCAAGAACGGGTTAGATGATTTTTTTGTCGTTAAATCCCTAAAAAACATGCAGGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGa                                                                                  <  2:608612/150‑1 (MQ=255)
                                                                  gcgcAAGCTTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGTTTAATGCCTGTAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtc                                                                             <  1:228028/151‑1 (MQ=255)
                                                                   cgaaaGCTTTTCCAAAAACAGGTTAGATGTTTTTTTTGCCGTTTAATGCCTGAAAACCATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGaaa                                                                                <  1:113735/144‑1 (MQ=255)
                                                                   cccAAGCTTTTTCCAAAAGATTTTAGATTATTTTTTTGTCGTTTAATGCCTGAAAAACATGCTTGACCCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAt                                                                               <  2:604008/146‑1 (MQ=255)
                                                                      aagttttttCAAGAACAGGTTAGATGTTTTTTTTGTCGTTTAATGCTTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGa                                                                                  <  1:906627/137‑1 (MQ=255)
                                                                      aaGCTTTTTCCAAAACAGTTAAGATGATTTTTTTGTCGCTTAATGCCTGAAAAACATGCAGGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCa                                                                                    <  1:53818/140‑1 (MQ=255)
                                                                      aaGCTTTTTCAAGAACGGTTAAGATATTTTTTTTGTCCTTTAATGCTTGAAAAACATGAATGCGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGa                                                                                  <  2:762243/142‑1 (MQ=255)
                                                                            tttccaaaaaaaGGTTAAATATTTTTTTTGTCTTTTATTCCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                                      <  2:223533/138‑1 (MQ=255)
                                                                            tttCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                   >  1:128783/1‑151 (MQ=255)
                                                                            tttCCAGAACAGGTTAATTGATTTTTTTTTCCCTTAATGCTTGAAAAACATGCATGACCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAg                                                                    <  1:272761/150‑1 (MQ=255)
                                                                             tcccAGAACGGTTTAGATTATTTTTTTGTCTCTTAACGCCTGAAAAACATGCGTGACCAAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                   <  2:71119/148‑1 (MQ=255)
                                                                              tCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                   >  1:228358/1‑149 (MQ=255)
                                                                               ccAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                   >  2:199716/1‑148 (MQ=255)
                                                                                  aGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAtt                                                              >  2:200429/1‑151 (MQ=255)
                                                                                    accgggTTAGATGATTTTTTTGGCCCTAAAGGCTTCAAAAAATTGCATGACCCAAAAAAAAAAAAAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTg                                                            <  1:494616/147‑1 (MQ=255)
                                                                                    aaCAGGTTAGATAATTTTTTTGTCGTTTAATGCCGGTAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAttt                                                             <  2:598672/150‑1 (MQ=255)
                                                                                        ggTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAtt                                                        >  1:987893/1‑151 (MQ=255)
                                                                                         gTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTg                                                       >  1:946441/1‑151 (MQ=255)
                                                                                          ttAAATGTTTTTTTTGTCGTTAAAGCCTTAAAAAACATGCAGGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTg                                                       <  1:751218/150‑1 (MQ=255)
                                                                                           taagaatattttttttTCGTTAAATGCCTGAAAAACAGGCATAGGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAtt                                                        <  1:602416/139‑1 (MQ=255)
                                                                                                 gatCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTCTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACg                                                >  1:813457/1‑151 (MQ=255)
                                                                                                 catCTTTTTGTCGCTTAATGCCTAAAAAACATGCATGCCCAAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAtt                                                    <  2:617235/146‑1 (MQ=255)
                                                                                                   tttttttGTCGTTTAATGCCTTTAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACg                                                <  2:513832/147‑1 (MQ=255)
                                                                                                    cTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                              >  1:374854/1‑150 (MQ=255)
                                                                                                    cTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                              >  2:956562/1‑150 (MQ=255)
                                                                                                     tttttGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcat                                            >  1:110185/1‑151 (MQ=255)
                                                                                                     tttttGTCGCTAATTGCCTTTAAAACATGCATGGCCAAAAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGt                                               <  2:305105/148‑1 (MQ=255)
                                                                                                     tttttGTCGTTAAAGCCTAAAAAAACATAAAGGAACCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                                 <  2:693513/146‑1 (MQ=255)
                                                                                                      tttgcccactaaaTGCCTGAAAACCATGCACGGCCCAAAAAATAAAATAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtca                                             <  1:1175005/138‑1 (MQ=255)
                                                                                                         tGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                                        >  1:288466/1‑151 (MQ=255)
                                                                                                           tCGCTTAATGCCTATAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtg                                       <  1:520109/150‑1 (MQ=255)
                                                                                                            cctttAATGCCTAAAAAACATGAATGAGCAAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtg                                       <  2:615578/146‑1 (MQ=255)
                                                                                                               ttAATGCCTGAAAAACATGAATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                                        <  1:618703/145‑1 (MQ=255)
                                                                                                                 aaTGCCTAAAAAACATGCATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                                 <  1:746600/150‑1 (MQ=255)
                                                                                                                 aaTCCCTGAAAAACATGAAGGACCCAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgtg                                     <  2:954338/146‑1 (MQ=255)
                                                                                                                 aaTCCCTGAAAAACATCAATGACCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                                 <  1:63353/150‑1 (MQ=255)
                                                                                                                  aTGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTa                                >  2:597605/1‑150 (MQ=255)
                                                                                                                   aatccttaaaaacCATGCATGAGCAAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                                 <  2:707331/141‑1 (MQ=255)
                                                                                                                   aaggcctaaaaaaCTTGAATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                                 <  1:769697/141‑1 (MQ=255)
                                                                                                                    gCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAat                              >  2:606905/1‑150 (MQ=255)
                                                                                                                      cTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                           >  1:402886/1‑151 (MQ=255)
                                                                                                                        ttaaaaaaCTTCCACGAGCCAAAAAATAATATAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAata                             <  1:380359/145‑1 (MQ=255)
                                                                                                                         tgaaaaaCATGCATGACCAAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTaa                               <  1:392318/143‑1 (MQ=255)
                                                                                                                         tgaaaaaCATCCATGACCAAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAatat                            <  2:337425/146‑1 (MQ=255)
                                                                                                                         tatAAAACAGGAAGGACCCACAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAatat                            <  1:1139066/147‑1 (MQ=255)
                                                                                                                           aaaaaCAGGCAGGACCAAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                           <  2:641530/147‑1 (MQ=255)
                                                                                                                            aaaaatggattgAGCAAAAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCa                         <  1:689906/139‑1 (MQ=255)
                                                                                                                             aaaCATGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                      <  2:341813/151‑1 (MQ=255)
                                                                                                                              aaCATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                      >  1:1164232/1‑150 (MQ=255)
                                                                                                                              aaCATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                      >  2:559805/1‑150 (MQ=255)
                                                                                                                              aaCATGCATGAGCCACAAAATAATATAAAAAATACCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                      >  2:311024/1‑150 (MQ=255)
                                                                                                                               aCATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa                     >  1:1052601/1‑150 (MQ=255)
                                                                                                                                catgcatgAGCAAAAAAAAAAAAAAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa                     <  2:956818/149‑1 (MQ=255)
                                                                                                                                    catgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                >  2:694047/1‑150 (MQ=255)
                                                                                                                                     atgaccaaaaaaaaaaaaaaaaaaaTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAag                  <  1:1113010/138‑1 (MQ=255)
                                                                                                                                     atgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                >  2:22173/1‑149 (MQ=255)
                                                                                                                                       gacccAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgag              <  1:1187947/146‑1 (MQ=255)
                                                                                                                                       gAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagt             <  2:445438/150‑1 (MQ=255)
                                                                                                                                         cccAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagt             <  2:505761/147‑1 (MQ=255)
                                                                                                                                          caaaaaaaTAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGaa          <  1:604006/146‑1 (MQ=255)
                                                                                                                                           cacaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat         >  1:1105758/1‑150 (MQ=255)
                                                                                                                                            aaaaaaaaaaaaaaaaaaCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtg            <  2:612182/144‑1 (MQ=255)
                                                                                                                                             aaaaaaaATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtg            <  1:194068/144‑1 (MQ=255)
                                                                                                                                             aaaaaaaATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAata        <  2:982700/148‑1 (MQ=255)
                                                                                                                                               aaaTAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtct    >  2:678141/1‑151 (MQ=255)
                                                                                                                                                aataatATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtctc   >  1:144522/1‑151 (MQ=255)
                                                                                                                                                aataatATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtctc   >  2:1015795/1‑151 (MQ=255)
                                                                                                                                                 ataataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTctct  >  2:360747/1‑151 (MQ=255)
                                                                                                                                                 ataataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTctct  >  2:340250/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                
AATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCT  >  NC_000913/29349‑29640

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