Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,448,424 T→A 5.2% intergenic (+189/‑23) ytfP → / → chpS GGCT‑like protein/antitoxin of the ChpBS toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,448,4240TA 5.2% 549.0 / 14.8 153intergenic (+189/‑23)ytfP/chpSGGCT‑like protein/antitoxin of the ChpBS toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (75/70);  new base A (6/2);  total (81/72)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

Only 100 of 158 total aligned reads displayed.

AGACGACTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCAACTGATT  >  NC_000913/4448277‑4448572
                                                                                                                                                   |                                                                                                                                                    
agaCGACTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTg                                                                                                                                                    <  2:2024946/150‑1 (MQ=255)
agaCGACTCGCATTCTGTTTTGGAATTCCCTCACCTTGTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAACAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTg                                                                                                                                                    <  1:1270656/150‑1 (MQ=255)
  acgacTCGCATTCTGTGTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCGATTAACGCAGAAAATACTCTGCTGGCCTTCACAAATGCGCAGGGGTAAAACGTGTCCTGTCGCACCGTGAGTTCTAGTTTGta                                                                                                                                                  >  2:457729/1‑150 (MQ=255)
    gacTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGtata                                                                                                                                                >  1:98291/1‑150 (MQ=255)
    gacTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCAATTAACGCATAAAAAACTCAGCTGGCATTCACAAAATCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGtaaa                                                                                                                                                >  1:98172/1‑148 (MQ=255)
     acTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATaa                                                                                                                                               <  1:1135067/150‑1 (MQ=255)
       tCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACt                                                                                                                                             <  1:776857/150‑1 (MQ=255)
         gCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACAAAGTATAACtca                                                                                                                                           <  2:1412840/150‑3 (MQ=255)
          cATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTaa                                                                                                                                          <  1:975604/150‑1 (MQ=255)
          cATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTaa                                                                                                                                          <  1:268805/150‑1 (MQ=255)
          cATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTaa                                                                                                                                          <  2:2089420/150‑1 (MQ=255)
           aTTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATTCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAag                                                                                                                                         >  1:1521154/1‑150 (MQ=255)
              cTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAaggag                                                                                                                                      <  1:1854080/150‑1 (MQ=255)
                 ttttGTAATTCCCTCACCTTTTGCTTTTCTCACCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTg                                                                                                                                   >  2:140140/1‑150 (MQ=255)
                         ttCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCg                                                                                                                           >  2:78770/1‑150 (MQ=255)
                         ttCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCg                                                                                                                           >  2:655613/1‑150 (MQ=255)
                            ccTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTAt                                                                                                                        >  2:847257/1‑150 (MQ=255)
                             cTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTAtt                                                                                                                       >  2:1125174/1‑150 (MQ=255)
                             cTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTAtt                                                                                                                       >  2:1431723/1‑150 (MQ=255)
                             cTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTAtt                                                                                                                       >  2:1749347/1‑150 (MQ=255)
                                 ccTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCa                                                                                                                   <  1:402514/150‑1 (MQ=255)
                                 ccTTTTGCTTTTCTCTCCGAGCCGCTTTCCATAACTATTAAAGCATAAAAAAATCAGCTGGCATTAACAAATGCGCAGGGGTAAAACGTTTCCGGTAGCACCGTGAGTTATACGTTTTATAAATTAAGGAGGTGCAGATGCGTATTAACa                                                                                                                   >  1:23389/1‑150 (MQ=255)
                                        cTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAg                                                                                                            >  2:1354899/1‑150 (MQ=255)
                                        cTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAg                                                                                                            >  2:2272290/1‑150 (MQ=255)
                                        cTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAg                                                                                                            >  2:1856218/1‑150 (MQ=255)
                                         ttttCTCTCCGAGCCGCTTTCCATATCTATTAAAGCATAAAAAAATATGCTGGCATTCACAAATGCGCAGGGGTAAAAAGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTAACAAAAAAAGa                                                                                                           >  1:511890/1‑150 (MQ=255)
                                            actcacCGAGCCGCTTTCCATAACTAATAACCCAAAAAAAACTCTGCTGGCAATCACAAATCCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGAGCAGATGCGTATTACCATAAAAAGATgg                                                                                                        <  2:21683/149‑1 (MQ=255)
                                             ctctcCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATggg                                                                                                       <  1:906364/150‑1 (MQ=255)
                                              tctcCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATgggg                                                                                                      <  1:1980478/150‑1 (MQ=255)
                                              tctcCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATgggg                                                                                                      <  2:1892759/150‑1 (MQ=255)
                                              tctcCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATgggg                                                                                                      <  1:1980308/150‑1 (MQ=255)
                                                tcCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGa                                                                                                    >  2:1115218/1‑150 (MQ=255)
                                                 ccGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGaa                                                                                                   <  2:629963/150‑1 (MQ=255)
                                                     gCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATTGGGGAACAGt                                                                                               >  1:553259/1‑150 (MQ=255)
                                                       cGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGc                                                                                             <  1:2226113/150‑1 (MQ=255)
                                                        gCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCa                                                                                            <  2:1283393/150‑1 (MQ=255)
                                                        gCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCa                                                                                            <  2:743274/150‑1 (MQ=255)
                                                        gCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCa                                                                                            <  2:244687/150‑1 (MQ=255)
                                                                 atCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTc                                                                                   >  2:696206/1‑150 (MQ=255)
                                                                     aTTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTc                                                                               <  1:2167634/150‑1 (MQ=255)
                                                                     aTTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTc                                                                               <  2:603456/150‑1 (MQ=255)
                                                                     aTTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTc                                                                               <  2:1744429/150‑1 (MQ=255)
                                                                     aTTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTc                                                                               <  1:1320993/150‑1 (MQ=255)
                                                                     aTTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTc                                                                               <  1:1941748/150‑1 (MQ=255)
                                                                      ttAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTcc                                                                              >  2:2102895/1‑150 (MQ=255)
                                                                      ttAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTcc                                                                              >  1:1241508/1‑150 (MQ=255)
                                                                      ttAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTcc                                                                              >  1:119448/1‑150 (MQ=255)
                                                                           gCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAata                                                                         >  2:2208627/1‑150 (MQ=255)
                                                                             aTAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATc                                                                       >  2:1679934/1‑150 (MQ=255)
                                                                             aTAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATc                                                                       >  2:1178536/1‑150 (MQ=255)
                                                                             aTAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGAGCAGGAATGGTCATTCCCAATATc                                                                       >  2:1178539/1‑150 (MQ=255)
                                                                              tAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCg                                                                      <  1:266872/150‑1 (MQ=255)
                                                                              tAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCg                                                                      <  1:266897/150‑1 (MQ=255)
                                                                               aaaaaaCTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGt                                                                     <  1:1963578/150‑1 (MQ=255)
                                                                                    aCTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGa                                                                <  2:412335/150‑1 (MQ=255)
                                                                                    aCTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGa                                                                <  2:412345/150‑1 (MQ=255)
                                                                                      tctGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGaaa                                                              <  2:379880/150‑1 (MQ=255)
                                                                                      tctGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGaaa                                                              <  2:379876/150‑1 (MQ=255)
                                                                                      tctGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGaaa                                                              >  1:1606119/1‑150 (MQ=255)
                                                                                          ctgGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGaac                                                          >  2:909944/1‑150 (MQ=255)
                                                                                           tgGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGaact                                                         >  1:359641/1‑150 (MQ=255)
                                                                                           tgGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGaact                                                         >  1:1649418/1‑150 (MQ=255)
                                                                                             gCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGaactta                                                       >  2:731426/1‑150 (MQ=255)
                                                                                                 tCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGaacttaactt                                                   >  1:1554962/1‑150 (MQ=255)
                                                                                                 tCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGaacttaactt                                                   >  1:219483/1‑150 (MQ=255)
                                                                                                  cacaAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAacttaactta                                                  >  1:1552990/1‑150 (MQ=255)
                                                                                                   acaAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTCCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTAc                                                 >  2:1827065/1‑150 (MQ=255)
                                                                                                   acaAATGCGCAGGGGTAAAACGTTTCCGGTAGCACCGTGAGTTATACGTTGTATACCTTAAGGAGGAGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGAAATGAAAGAACTTAACTTAc                                                 >  2:1827057/1‑150 (MQ=255)
                                                                                                    caAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACa                                                <  1:1685392/150‑1 (MQ=255)
                                                                                                        tGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCg                                            <  1:2030405/150‑1 (MQ=255)
                                                                                                         gcgcAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCgg                                           <  1:1924339/150‑1 (MQ=255)
                                                                                                           gcAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATCTCGTAATGAAAGAACTTAACTTACAGCCgggg                                         >  2:783099/1‑150 (MQ=255)
                                                                                                           gcAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCgggg                                         <  2:1559873/150‑1 (MQ=255)
                                                                                                           gcAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCgggg                                         <  2:1559910/150‑1 (MQ=255)
                                                                                                             aGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCa                                       <  1:1837352/150‑1 (MQ=255)
                                                                                                                      aCGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGa                              <  2:1489651/150‑1 (MQ=255)
                                                                                                                        ggttCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGCTGCGTAGTACCATCACAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGCAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAgg                            <  2:1410922/148‑1 (MQ=255)
                                                                                                                         tttCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGgc                           >  1:2102016/1‑150 (MQ=255)
                                                                                                                         tttCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGgc                           >  1:557289/1‑150 (MQ=255)
                                                                                                                         tttCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCAGAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGgc                           >  2:933893/1‑150 (MQ=255)
                                                                                                                             cTGTAGCACCGTGATTTATACTTTGTATTACTTAATGAGGTGCAGATGCGTATTACATTAAAAAGTTGTGGGAACAGTGCAGGTATGGTCTTTACCTTTATCGTAATGAAAGAACTTAACTTTCAGCCGGGGCTGAGCGTGGAggatcta                       <  2:2383317/150‑6 (MQ=255)
                                                                                                                             cTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCaa                       <  2:2009977/150‑1 (MQ=255)
                                                                                                                             cTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCaa                       <  1:2360473/150‑1 (MQ=255)
                                                                                                                             cTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCaa                       <  1:1250309/150‑1 (MQ=255)
                                                                                                                                     ccGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGcaa               <  2:1097446/150‑1 (MQ=255)
                                                                                                                                      cGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAGAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGcaac              >  1:2046615/1‑150 (MQ=255)
                                                                                                                                      cGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGcaac              >  1:2046404/1‑150 (MQ=255)
                                                                                                                                       ttaaGTTATACAAAGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGcaaca             >  2:1391033/4‑150 (MQ=255)
                                                                                                                                           gTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCa         >  2:2074068/1‑150 (MQ=255)
                                                                                                                                           gTTATACAAAGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCa         >  1:1580667/1‑150 (MQ=255)
                                                                                                                                           gTTATACAAAGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCa         >  1:607320/1‑150 (MQ=255)
                                                                                                                                           gTTATACAAAGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCa         >  1:854279/1‑150 (MQ=255)
                                                                                                                                            ttATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCaa        >  2:2384543/1‑150 (MQ=255)
                                                                                                                                            ttATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCaa        >  1:60589/1‑150 (MQ=255)
                                                                                                                                               taCTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCAACTg     >  1:1014062/1‑150 (MQ=255)
                                                                                                                                                 cTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCAACTGAt   >  2:2313438/1‑150 (MQ=255)
                                                                                                                                                  tttGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCAACTGAtt  <  2:1886900/150‑1 (MQ=255)
                                                                                                                                                  tttGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCAACTGAtt  <  1:794304/150‑1 (MQ=255)
                                                                                                                                                  tttGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCAACTGAtt  <  2:1226180/150‑1 (MQ=255)
                                                                                                                                                  tttGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCAACTGAtt  <  2:952033/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGACGACTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCAACTGATT  >  NC_000913/4448277‑4448572

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