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

TTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAGG  >  NC_000913/29347‑29645
                                                                                                                                                        |                                                                                                                                                     
ttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAatat                                                                                                                                                       <  2:929370/150‑1 (MQ=255)
     gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATaaaaaa                                                                                                                                                 <  2:1033964/151‑1 (MQ=255)
       aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAAt                                                                                                                                                <  2:847522/150‑1 (MQ=255)
         aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATccc                                                                                                                                             <  1:753721/151‑1 (MQ=255)
         aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATccc                                                                                                                                             <  1:104283/151‑1 (MQ=255)
              gACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCAt                                                                                                                                        >  1:836438/1‑151 (MQ=255)
                 cATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                                      >  2:708452/1‑150 (MQ=255)
                      gtcccttttttttCTGCTCTTTTTTTTTTCATGAAATTTTTTTGCTGCGAAAGCTTTTCCA‑AAAAGAGGTT‑AATTGTTTTTTTTGTCGTTTAA‑GCCTTAAAAACCATACATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGa                                                                                                                                <  1:1151743/144‑1 (MQ=255)
                            tttttttCTGCTCTTTTTTTTTTCAGGCAATTTTTTTGCTGCGAAGGCTTTTCCC‑AGAACGGGTT‑AGATGATTTTTTTGTCGCTTATTCCTTATAAAA‑CATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGAc                                                                                                                               <  2:1022410/145‑1 (MQ=255)
                             ttttttCTGCTCTTTTTTTTTTCATGCAATCTTCTTGCTGCGCAAGCTTTTTCA‑AGAACAGGTT‑AGATAATTTTTTTTTCGCTTAATCCCTTAAAAA‑CATGAATGACCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                                           <  1:1086345/149‑1 (MQ=255)
                              tttttCTGCTCGTTTTTATTTTATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgc                                                                                                                        >  1:83550/1‑151 (MQ=255)
                                 ttCTGCTCTTTTTTATTTCATGCAATCTTCTTGCTGCGCAGGCTTTTTCC‑AAAACAGGTT‑AGATGATTTTTTTGTCGTTTAATGCTTGAAAAA‑CATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAgcgc                                                                                                                      <  2:704210/150‑1 (MQ=255)
                                       tcttttttATTTCATGCATTCTTCTTGCTGCGCAAGCGTTTCCA‑AAAACGGGTT‑AGATGATTTTTTTGTCGTTTAATGCCTGAAAAA‑AATGCATGACCAAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGccc                                                                                                                    <  2:1045200/143‑1 (MQ=255)
                                       tCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtc                                                                                                               >  2:742152/1‑151 (MQ=255)
                                        cGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtct                                                                                                              >  2:493082/1‑151 (MQ=255)
                                            tttATTTCATGCATTCTTCTTGCTGCGAAAGCTTTTTCC‑AGAACAGTTT‑AAATGATTTTTTTGTCGTTTAATGCCTGAAAAA‑CATGCATGGGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCa                                                                                                           <  1:1180003/150‑1 (MQ=255)
                                             ttttttAATGAATTCTTCTTGCTCCGCAACCTTTTTCA‑AGAACAGTTT‑AGATGATTTTTTTTTCGTTTAATCCCTAAAAAA‑CTTGCATGCGCAAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAg                                                                                                          <  1:535950/147‑1 (MQ=255)
                                               aTTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAAt                                                                                                       >  2:202665/1‑151 (MQ=255)
                                                tttCATGCAATTTTCTTGCTGCGCAAGCGTTTTCA‑AGAACAGGTT‑AGATGATTTTTTTGTCGTTTAATGCCTGAAAAA‑CATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGaa                                                                                                        <  2:408089/149‑1 (MQ=255)
                                                tttCATGAAATCTTCTTGCTGCGAAAGCTTTTTCC‑AGAACAGGTA‑AGATGATCTTTTTGTCGCTTATTGCTTGTAAAA‑CATGCATGACCAACAAAAAAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                                         <  2:1170695/148‑1 (MQ=255)
                                                tttCAGGAAATCTTTTGGCTGCGCAAGCGTTTTCA‑AGAACAGGTA‑AGTTGTTTTTTTTGTCCTTTAAGGTCTGAAAAA‑CATGCATGAACCAAAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAAt                                                                                                       <  2:167927/150‑1 (MQ=255)
                                                          tttttttGCTGCGAAAGC‑TTTTCC‑AAAAAGGGTTAAAATGTTTTTTTTTCCGTTTACGGCCTGAAAAA‑CATGAAGGGCCAAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGtt                                                                                              <  2:370221/144‑1 (MQ=255)
                                                          ttttcttGCTGCGCAAGCTTTTCCC‑AGAACAGTTT‑AGATGATTTTTTTGTCGCTTAATCCTTATAAAA‑CATGCTGGAGCAAAAAAAAAATAAAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgcc                                                                                                 <  1:1164313/144‑1 (MQ=255)
                                                            tttttGCGGCGCAGGCTTTTCCC‑AAAACAGTTA‑AGATGTTTTTTTTGTCGTTAAATGCCTGAAAAA‑CAGGCAGGCCCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGtt                                                                                              <  2:811269/144‑1 (MQ=255)
                                                              cttGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                          >  1:1152413/1‑149 (MQ=255)
                                                                tgctgcGCAAGCTTTTCCC‑AGAACAGTTA‑AGATTTTTTTTTTGTCGCTTAATGCCTAAAAAA‑CATGCATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGa                                                                                       <  1:29026/150‑1 (MQ=255)
                                                                tgctgcGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGa                                                                                       >  2:1062858/1‑150 (MQ=255)
                                                                      cgaaagtttttCCAAGAACAGTTT‑AAATGATTTTTTTGTCGCTTAATCCTTAAAAAA‑CATCCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtc                                                                                  <  1:339696/143‑1 (MQ=255)
                                                                        aaGCTTTTCCA‑AAAACAGTTT‑AGATGTTTTTTTTGTCCTTTAATGCCTGAAAAA‑CATGAATGAGCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                              <  2:158686/151‑1 (MQ=255)
                                                                        aaGCGTTTTCA‑AAAACAGGTT‑AGATGATTTTTTTGTCGCTTAATGCCTGAAAAA‑CATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                              <  2:1203658/151‑1 (MQ=255)
                                                                            gTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCgg                                                                             >  2:467532/1‑148 (MQ=255)
                                                                             ttttCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                                           >  2:275139/1‑149 (MQ=255)
                                                                              tttCA‑AGAACAGGTT‑AGATGATTTTTTTGTCGCTAAATGCCTGAAAAA‑CATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAg                                                                         <  2:45928/150‑1 (MQ=255)
                                                                               ttCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                        >  1:236882/1‑150 (MQ=255)
                                                                                tCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                        >  2:247667/1‑149 (MQ=255)
                                                                                     gAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAttt                                                                  >  1:901668/1‑151 (MQ=255)
                                                                                      aaCAGGTT‑AGATGATTTTTTTGTCTTTTAATCCCTAAAAAA‑CATGCATGAGCCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTg                                                                 <  1:734283/151‑1 (MQ=255)
                                                                                           gttaaGATGTTTTTTTTGCCCTTTAATGCCTGAAAAA‑AATGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCa                                                               <  2:303062/145‑1 (MQ=255)
                                                                                             t‑AGATCATTTTTTTGTCGCTAAATGCTTGAAAAA‑CATGCAGGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGa                                                           <  2:1023706/150‑1 (MQ=255)
                                                                                               aGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAtt                                                         >  1:841939/1‑151 (MQ=255)
                                                                                                      cTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                                   >  2:322055/1‑150 (MQ=255)
                                                                                                       tttttGTCGTTTAATCCCTGAAAAA‑CATGCATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACg                                                     <  1:609241/147‑1 (MQ=255)
                                                                                                       tttttGGCGTTAAATGCCTAAAAAA‑CATGAATGAGCAAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                                   <  2:63435/149‑1 (MQ=255)
                                                                                                          ttGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTcatcat                                              >  1:521641/1‑151 (MQ=255)
                                                                                                               gCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgtg                                          >  1:394790/1‑150 (MQ=255)
                                                                                                                 ttAATGCCTATAAAA‑CATGCATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGa                                         <  2:1192712/149‑1 (MQ=255)
                                                                                                                    atcccTTTAAAA‑AATGCATGGCCCACAAAATAATAAAAAAAACTCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTaa                                    <  2:939329/148‑1 (MQ=255)
                                                                                                                       ccTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAatat                                 >  1:1174232/1‑151 (MQ=255)
                                                                                                                         tGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGc                               >  2:785023/1‑151 (MQ=255)
                                                                                                                            aaaa‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaa                             >  2:1053708/1‑150 (MQ=255)
                                                                                                                            aaaa‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaa                             >  2:536620/1‑150 (MQ=255)
                                                                                                                             aaaacttgaagggcccAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAata                                  <  2:675026/131‑1 (MQ=255)
                                                                                                                             aaa‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                           >  2:221908/1‑151 (MQ=255)
                                                                                                                             aaa‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                           >  2:444484/1‑151 (MQ=255)
                                                                                                                              aa‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                           >  1:1123108/1‑150 (MQ=255)
                                                                                                                                 catgcatgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa                         >  2:740565/1‑150 (MQ=255)
                                                                                                                                      atgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                     >  2:85527/1‑149 (MQ=255)
                                                                                                                                       tgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgag                   >  2:936953/1‑150 (MQ=255)
                                                                                                                                         acccAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgag                   <  1:1162399/146‑1 (MQ=255)
                                                                                                                                               aaaaaaATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat              <  2:56405/147‑1 (MQ=255)
                                                                                                                                               aaaaTAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc          >  1:633264/1‑151 (MQ=255)
                                                                                                                                               aaaaTAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc          <  1:551328/151‑1 (MQ=255)
                                                                                                                                                  ataataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTctct       >  2:702847/1‑151 (MQ=255)
                                                                                                                                                   taataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTg      >  2:1124737/1‑151 (MQ=255)
                                                                                                                                                    aataAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg     <  1:843709/151‑1 (MQ=255)
                                                                                                                                                       ataAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAgg  >  2:1149197/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
TTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAA‑CATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAGG  >  NC_000913/29347‑29645

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: ‑