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

TAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGA  >  NC_000913/29348‑29643
                                                                                                                                                     |                                                                                                                                                   
taagggaaaaTTTTAACATTTGGTGCGCTGTTGTCTGCTCGTGTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAAATGATCTTTTTGTCGCTTGATGCCTGTAAAACATGGATGAGCCGCAAAATAAtata                                                                                                                                                    <  1:586521/144‑1 (MQ=255)
tAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAAtata                                                                                                                                                    <  2:1037090/150‑1 (MQ=255)
tAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAAtata                                                                                                                                                    <  2:1037395/150‑1 (MQ=255)
 aaTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaa                                                                                                                                                  <  2:415761/151‑1 (MQ=255)
     tAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaaaaa                                                                                                                                               <  1:18690/150‑1 (MQ=255)
     tAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAt                                                                                                                                              >  2:1046868/1‑151 (MQ=255)
      aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATc                                                                                                                                             >  2:1027196/1‑151 (MQ=255)
      aaaTTTTGACCATTTGGTCCACTTTTTTCTACTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATGATATaaaaaa                                                                                                                                               <  1:111321/149‑1 (MQ=255)
       aaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCGGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATcc                                                                                                                                            <  1:1117086/151‑1 (MQ=255)
        aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATccc                                                                                                                                           >  2:639749/1‑151 (MQ=255)
             gACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCa                                                                                                                                       <  2:247425/150‑1 (MQ=255)
               ccATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCAtt                                                                                                                                     >  2:185152/1‑150 (MQ=255)
                        ccACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACtt                                                                                                                           >  1:173772/1‑151 (MQ=255)
                           cTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                                         >  2:365979/1‑150 (MQ=255)
                                ttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGcc                                                                                                                   >  2:359128/1‑151 (MQ=255)
                                  cTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGccc                                                                                                                  >  1:903024/1‑150 (MQ=255)
                                    gCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAATCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATCTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                                >  1:1071336/1‑150 (MQ=255)
                                    gCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                                >  2:988091/1‑150 (MQ=255)
                                     cTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGATAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATCAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCatat                                                                                                              >  2:666032/1‑151 (MQ=255)
                                      tcttttttATTTCATGCATTCTTCTTGCTGCGCAAGCGTTTTCAAGAAC‑AGTTTAGATGATTTTTTTGTCGTTTAATGCTTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                                <  2:435722/145‑1 (MQ=255)
                                        ttttttATTTCATAAATTTTTCTTGCTGCGCAATCTTTTTCCAAAAC‑AGTTAAGATAATTTTTTTTTCGTTTAATCCTTGAAAAACATGCTTGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCata                                                                                                               <  1:577225/146‑1 (MQ=255)
                                        gTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                           >  1:903075/1‑151 (MQ=255)
                                         tttttATTTCATGAAATCTTCTTGCTGCGCAAGTGTTTCCCAGAAC‑GGGTTAGATGTTTTTTTTGTCGTTAAATGCCTAAAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtct                                                                                                            <  2:1099375/149‑1 (MQ=255)
                                          ttttATTTCATCCAATCTTTTTGCTGCGCAAGCGTTTCCCAGAAC‑GGTTTAAATGATTTTTTTGTCGCTAAATGCCTTAAAAACATGCAGGACCCAAAAAAAAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                          <  2:1054544/150‑1 (MQ=255)
                                           tttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                           >  1:581533/1‑148 (MQ=255)
                                             tATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                                       >  2:724555/1‑150 (MQ=255)
                                                    taaattcttctTGCCGCGCAAGCTTTTCACAAAAC‑AGTTTAAATTATTTTTTTGTCGTTTAAGGCCTAAAAAACATCAAGGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccg                                                                                                 <  1:628651/144‑1 (MQ=255)
                                                         tcttctTGCTGCGCAAGCTTTTTCCAAAAC‑AGGTTAGATGATTTTTTTGTCGTTTAATGCCTGAAAAACATGAATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgc                                                                                                <  2:167449/145‑1 (MQ=255)
                                                          cttcttGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGc                                                                                         >  1:280568/1‑151 (MQ=255)
                                                            tcttGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                        >  2:264514/1‑150 (MQ=255)
                                                            tcttGCTGCGCAAGC‑TTTTCCAAAACGAGTTTAGATGATTTTTTTGTCGTTAAATGCCTAAAAAACAGGCAGGAGCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGttt                                                                                           <  2:131349/147‑1 (MQ=255)
                                                                gctgcGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGaa                                                                                    >  2:106752/1‑150 (MQ=255)
                                                                  tgcGCAAGCTTTTTCAAGAAG‑AGGTAAAATAATTTTTTTTTCGTTTAATCCCTTAAAAACATGCTTGGCCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGa                                                                                     <  1:288061/147‑1 (MQ=255)
                                                                      cAAGCTTTTTCCAGAAC‑AGGGTAGATGATCTTTTTGTCGCTTAATGCCTGTAAGACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGCTCGCCAGAAATtcg                                                                               >  1:497186/1‑149 (MQ=255)
                                                                       aaGCGTTTTCCAGAAC‑AGGTTAGATGTTTTTTTTGTCGTTTAATGCCTGAAAAACAGGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc                                                                             <  2:989244/150‑1 (MQ=255)
                                                                       aaGCGTTTTCCAAAAC‑GGTTTAGATGTTTTTTTTGTCGCTTAATGCCTGAAAAACATGAATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgt                                                                              <  2:798690/149‑1 (MQ=255)
                                                                       aaGCGTTTCCCAAACC‑AGGTTAGATGTTTTTTTTGTCCCTTAATGCCTAAAAAACATGCGTGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc                                                                             <  2:660645/150‑1 (MQ=255)
                                                                            ttttCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                                         >  1:766541/1‑149 (MQ=255)
                                                                            ttttCCAAAAC‑AGTTTAGATGTTTTTTTTGTCCCTAAATCCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                                         <  2:818274/149‑1 (MQ=255)
                                                                            ttttCAAGAAC‑AGTTTAATTTATATTTTTGTCTTTAAATCCCTTAAAAAATTGCATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGt                                                                          <  1:131779/148‑1 (MQ=255)
                                                                            tttcCCAGAAC‑AGGTTAGATGATTTTTTTGTCGTTTATTGCTTGAAAAACAGGCATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTaa                                                                        <  2:464588/150‑1 (MQ=255)
                                                                             tttCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAg                                                                       >  2:200480/1‑150 (MQ=255)
                                                                              ttCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                      >  2:449266/1‑150 (MQ=255)
                                                                              tccaaGAAC‑AGTTTAGATGATTTTTTTGTCGCTTAATGCCTGAAAAACATGCATGAGCCAAAAAAAAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                      <  2:231382/146‑1 (MQ=255)
                                                                                  aaaaC‑AGTTTAAATGATTTTTTTGTCTCTTAATCCCTGTAAAACATGCATGACCCAAAAAAAAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                      <  1:419581/144‑1 (MQ=255)
                                                                                  aGAAC‑AGTTTAGATGTTTTTTTTGTCGTTTAATGCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAg                                                                    <  2:594781/148‑1 (MQ=255)
                                                                                  aGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAtt                                                                 >  2:642107/1‑151 (MQ=255)
                                                                                  aGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAtt                                                                 >  2:357625/1‑151 (MQ=255)
                                                                                    aaC‑AGGTTAGATGTTTTTTTTGTCGTTTAATGCCTGAAAAACATCCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGa                                                                   <  2:74954/147‑1 (MQ=255)
                                                                                        aGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCAAATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAt                                                            >  1:789628/1‑151 (MQ=255)
                                                                                           ttAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTg                                                          >  1:1126842/1‑150 (MQ=255)
                                                                                                     tttttGTCGCTTAATGCTTAAAAAACATGCATGACCAAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                                 <  2:72490/149‑1 (MQ=255)
                                                                                                     tgtttGTCGTTTATTCCTTGAAAAACATGCATTACCAACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                                 <  1:567845/147‑1 (MQ=255)
                                                                                                      ttttttcaattAATCCTTGTAAAACATGCTAGCGCAACAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGt                                                  <  1:146936/138‑1 (MQ=255)
                                                                                                         tttCGCTAAATGCTTAAAAAACATGTAGGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTcatcat                                            <  2:418340/148‑1 (MQ=255)
                                                                                                          gTCGTTTAATGCCTATAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtg                                          <  2:1110358/151‑1 (MQ=255)
                                                                                                               ttAATGACTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAt                                     >  1:469845/1‑151 (MQ=255)
                                                                                                                 attcccTGAAAAATATGAAGGACCCACAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAt                                     <  1:709528/145‑1 (MQ=255)
                                                                                                                   tgcataaaaaaTATGCATGTGCAAAAAAAAAAAAAAAAAAAGCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTTATTTACGTCATCATTGTGAATTaa                                  <  1:13028/143‑1 (MQ=255)
                                                                                                                    gCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTa                                   >  1:413207/1‑148 (MQ=255)
                                                                                                                    ccctaaaaaaCAGGCAAGACCAAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTa                                   <  1:357367/142‑1 (MQ=255)
                                                                                                                      cTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                              >  2:281345/1‑151 (MQ=255)
                                                                                                                          aaacaatgaatgaACCAAAAAAAAAATAAAAAAAACCCGCCAGTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaaa                          <  1:1025189/146‑1 (MQ=255)
                                                                                                                             aCATGCATGACCAAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa                       <  2:763328/151‑1 (MQ=255)
                                                                                                                                tgcatgAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaaa                      <  2:24626/149‑1 (MQ=255)
                                                                                                                                   ctgATCCACAAACTAATATAAAAAATCCAGCCCTTAAGTTGACTTTTAGCGCCAATATATCCAGAATGCCGCCGTTTGCAAGCAATTCTCCGTTACGCAGATTTGCATTGCTTTCCGTCATCCTTGTGAATGAATATGCAAATAAagtg                   >  1:809417/2‑149 (MQ=255)
                                                                                                                                    tgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgag                 >  1:907277/1‑150 (MQ=255)
                                                                                                                                        ccAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa                       <  1:185099/140‑1 (MQ=255)
                                                                                                                                          acaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat          >  1:459601/1‑151 (MQ=255)
                                                                                                                                           caAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAtt         >  2:176531/1‑151 (MQ=255)
                                                                                                                                           aaaaaTAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAtt         >  1:374961/2‑151 (MQ=255)
                                                                                                                                           aaaaaTAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAtt         <  2:98393/150‑1 (MQ=255)
                                                                                                                                            aaaaaaaaaaaaaaaaTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc        <  1:492045/151‑1 (MQ=255)
                                                                                                                                            aaaaaaaaaaaaaaATCCCCGCCATTAAGTTGACTTTTAGCG‑CCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtct       <  2:1103025/151‑1 (MQ=255)
                                                                                                                                              aataatATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtct       >  1:553293/1‑150 (MQ=255)
                                                                                                                                              aataatAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtctc      <  2:478732/151‑1 (MQ=255)
                                                                                                                                                 aaaataAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg   <  1:13321/148‑1 (MQ=255)
                                                                                                                                                  atatAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGa  <  1:553621/151‑1 (MQ=255)
                                                                                                                                                   tataAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg   >  2:998188/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAAC‑AGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGA  >  NC_000913/29348‑29643

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