Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,448,425 T→A 5.9% intergenic (+190/‑22) 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,4250TA 5.9% 378.7 / 6.4 127intergenic (+190/‑22)ytfP/chpSGGCT‑like protein/antitoxin of the ChpBS toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (59/60);  new base A (5/3);  total (64/63)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

GAGACGACTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCAACTGATTCT  >  NC_000913/4448276‑4448574
                                                                                                                                                     |                                                                                                                                                     
gagaCGACTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACttt                                                                                                                                                       <  2:1623027/150‑1 (MQ=255)
   acgacTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGATGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGta                                                                                                                                                    >  1:275625/1‑150 (MQ=255)
    cgacTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGtat                                                                                                                                                   >  2:254151/1‑150 (MQ=255)
    cgacTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGAAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGtat                                                                                                                                                   >  1:457890/1‑150 (MQ=255)
       cTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAAc                                                                                                                                                <  1:254151/150‑1 (MQ=255)
         cGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACtt                                                                                                                                              <  1:92269/150‑1 (MQ=255)
         cGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACtt                                                                                                                                              >  2:1388316/1‑150 (MQ=255)
               cTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAACACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGATTTATACTTTGTATAACTTAaggag                                                                                                                                        >  2:1089607/1‑150 (MQ=255)
                tGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAaggagg                                                                                                                                       <  1:333032/150‑1 (MQ=255)
                 gTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCGTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGt                                                                                                                                      >  1:1788856/1‑150 (MQ=255)
                 gTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGt                                                                                                                                      >  1:1392713/1‑150 (MQ=255)
                 gTTTTGTAAGTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAACCTCTGCTGGCATTCACAAAAGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGt                                                                                                                                      >  2:1012030/1‑150 (MQ=255)
                  ttttGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTg                                                                                                                                     >  1:1384543/1‑150 (MQ=255)
                  ttttGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTg                                                                                                                                     >  1:975201/1‑150 (MQ=255)
                     tGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAg                                                                                                                                  <  2:379448/150‑1 (MQ=255)
                      gTAATTCCCTCACCATTTGCTTTTCTCTCCGCGCCGCTTTGCATATCGATTAACCCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAACCGTTTCCTGTAGCACCGTGAGTTTTACTTTGTATAACTTAAGGAGGTGCAGa                                                                                                                                 >  2:845205/1‑150 (MQ=255)
                       tAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGAt                                                                                                                                <  2:1497603/150‑1 (MQ=255)
                         aTTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTAACAAATGCGCAGGGGAAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGc                                                                                                                              >  2:1433054/1‑150 (MQ=255)
                            cccTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGGGCAGATGCGTa                                                                                                                           <  1:1132056/150‑1 (MQ=255)
                             ccTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTAt                                                                                                                          <  1:608383/150‑1 (MQ=255)
                             ccTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTAt                                                                                                                          >  1:46294/1‑150 (MQ=255)
                              cTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTACCGCATAAAAAACTCTGCTGGCATTCCCAAAAGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGAGGCGTAtt                                                                                                                         >  1:345416/1‑150 (MQ=255)
                                cACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTAc                                                                                                                       <  1:126457/150‑1 (MQ=255)
                                cACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTAc                                                                                                                       <  1:1595043/150‑1 (MQ=255)
                                cACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTAc                                                                                                                       >  1:1838604/1‑150 (MQ=255)
                                cACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTAc                                                                                                                       <  1:48501/150‑1 (MQ=255)
                                cACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTCCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGAATAACTAAAGGAGGTGCAGATGCGTATTAc                                                                                                                       >  1:1186536/1‑150 (MQ=255)
                                 aCCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTAcc                                                                                                                      <  1:779226/150‑1 (MQ=255)
                                 aCCTTTAGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTAcc                                                                                                                      <  1:1744915/150‑1 (MQ=255)
                                   cTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCAt                                                                                                                    <  1:1490552/150‑1 (MQ=255)
                                         cTTTTCTCTCCGAGCCGATTTCCATATCTATTAACGCATAAAAAACTCTCCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTCTAGCACCGTGAGTTATACTTTGTAAAACTTAAGGAGGTGCAGATGCGTATAACCAGAAAAAg                                                                                                              >  1:285995/1‑150 (MQ=255)
                                               tctcCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTGTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATgggg                                                                                                        <  2:280472/150‑1 (MQ=255)
                                                  ccGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGaa                                                                                                     >  1:1167225/1‑150 (MQ=255)
                                                  ccGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGAATTCACAAATTAGCAGGGGTAAAACGTGTCCTGTAGCACCGTGAGTTATACTTAGTATAACTTAAGGAGGTGCAGATGCGTATTACAATAAAAAGATGGGGGaa                                                                                                     >  2:305786/1‑150 (MQ=255)
                                                    gAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACa                                                                                                   >  1:653676/1‑150 (MQ=255)
                                                     aGCCTCTTTCCAAATCTAATAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTTAGTTATACTTTGTATAACTTAAGGAGGTGCAGATTCGTATTACCATAAAAAGATGGGGGAACAg                                                                                                  <  1:1839858/150‑1 (MQ=255)
                                                     aGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAg                                                                                                  <  1:1847115/150‑1 (MQ=255)
                                                      gCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGt                                                                                                 >  1:643609/1‑150 (MQ=255)
                                                         gCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATTGGGGAACAGTGCa                                                                                              >  1:1696684/1‑150 (MQ=255)
                                                          cTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAg                                                                                             <  1:1567945/150‑1 (MQ=255)
                                                           tttCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAgg                                                                                            <  2:457890/150‑1 (MQ=255)
                                                             tCCATATCTATTAACGCATAAAAAACTCTGCTGGTATTCACAAATGCGCAGGGGTAAAACGTTTCATGTAGCACCGTGAGTTATAATTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTa                                                                                          >  1:917247/1‑150 (MQ=255)
                                                              ccATATCTATTAACGAATAAAAAACTCTGCTGGCATTCACAAAAGCGCAGGGGTAAAACGTTTCGTGTAGCGCCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGAGTATTACAATAAAAAGATTGTGGAACTGTGCAGGTAt                                                                                         >  2:1811894/1‑150 (MQ=255)
                                                                   tCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCa                                                                                    >  2:1423948/1‑150 (MQ=255)
                                                                      aTTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTc                                                                                 >  1:363724/1‑150 (MQ=255)
                                                                        tAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTccc                                                                               <  2:1310900/150‑1 (MQ=255)
                                                                         aaCGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCa                                                                              >  2:769552/1‑150 (MQ=255)
                                                                          aCGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCaa                                                                             <  2:1833237/150‑1 (MQ=255)
                                                                           cGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCaaa                                                                            <  2:1107559/150‑2 (MQ=255)
                                                                            gCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAACACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAATAGATGGGGGAACAGTGCAGGTATGGGCATTCCCAata                                                                           <  1:909849/150‑1 (MQ=255)
                                                                            gCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAata                                                                           >  1:1136559/1‑150 (MQ=255)
                                                                            gCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAata                                                                           <  1:1751851/150‑1 (MQ=255)
                                                                             cATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAAATTTAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAatat                                                                          >  2:861167/1‑150 (MQ=255)
                                                                                 aaaaaCTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTa                                                                      >  2:1771418/1‑150 (MQ=255)
                                                                                   aaaCTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAAt                                                                    <  2:1659662/150‑1 (MQ=255)
                                                                                        ctgctgGCATTCACAACTGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAg                                                               <  2:311904/150‑1 (MQ=255)
                                                                                        ctgctgGCATTCACAAATGCGCAGGGGTACAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCCGGTATGGTCATTCCCAATATCGTAATGAAAg                                                               <  2:948295/150‑1 (MQ=255)
                                                                                        ctgatgGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAg                                                               <  2:807660/150‑1 (MQ=255)
                                                                                           ctgGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGaac                                                            >  2:1354558/1‑150 (MQ=255)
                                                                                           ctgGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGaac                                                            <  2:1177773/150‑1 (MQ=255)
                                                                                             ggCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGaactt                                                          >  1:1500464/1‑150 (MQ=255)
                                                                                               cATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGTGGAACAGTGCAGGTATGGTCATTCCCAATATATTAATGAAATaaattaa                                                        >  1:334867/1‑150 (MQ=255)
                                                                                                   cacaAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAacttaactta                                                    >  1:1779433/1‑150 (MQ=255)
                                                                                                    acaAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGTGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTAc                                                   <  2:1431146/150‑1 (MQ=255)
                                                                                                       aaTGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGc                                                >  1:1380729/1‑150 (MQ=255)
                                                                                                         tGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCg                                              >  1:1340255/1‑150 (MQ=255)
                                                                                                           cgcAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCggg                                            >  2:642050/1‑150 (MQ=255)
                                                                                                           cgcAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCggg                                            >  2:394192/1‑150 (MQ=255)
                                                                                                           cgcAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCggg                                            <  2:1136559/150‑1 (MQ=255)
                                                                                                             cAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGc                                          <  2:1482409/150‑1 (MQ=255)
                                                                                                             cAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGc                                          <  2:356162/150‑1 (MQ=255)
                                                                                                             cAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGc                                          <  1:1264779/150‑1 (MQ=255)
                                                                                                              aGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCa                                         <  1:1569267/150‑1 (MQ=255)
                                                                                                                 ggTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCagag                                      <  2:46294/150‑1 (MQ=255)
                                                                                                                     aaaCGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTg                                  <  2:247673/150‑1 (MQ=255)
                                                                                                                     aaaCGATTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTg                                  <  2:812560/150‑1 (MQ=255)
                                                                                                                       aCGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGa                                <  1:392692/150‑1 (MQ=255)
                                                                                                                        cGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAg                               >  2:1721695/1‑150 (MQ=255)
                                                                                                                          tatCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGgc                             <  2:517902/148‑1 (MQ=255)
                                                                                                                              cTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCaa                         >  1:842471/1‑150 (MQ=255)
                                                                                                                               tGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAg                        >  2:794693/1‑150 (MQ=255)
                                                                                                                                 tAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTg                      >  1:838151/1‑150 (MQ=255)
                                                                                                                                 tAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTg                      <  1:1804487/150‑1 (MQ=255)
                                                                                                                                 tAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGAGTGGAGGCGCAAGTg                      >  1:22205/1‑150 (MQ=255)
                                                                                                                                  aGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGa                     >  1:629021/1‑150 (MQ=255)
                                                                                                                                  aGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGa                     <  2:764168/150‑1 (MQ=255)
                                                                                                                                   gCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGTAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAg                    <  1:1433054/150‑1 (MQ=255)
                                                                                                                                   gCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAg                    <  2:1540658/150‑1 (MQ=255)
                                                                                                                                   gCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTAACATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAATTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAg                    <  1:1309441/150‑1 (MQ=255)
                                                                                                                                    cACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTAAAGCCGGGGCACAGAGTGGAGGCGCAAGTGAGc                   >  1:314717/1‑150 (MQ=255)
                                                                                                                                     tccttaaGTTATACAAAGTATAACTTAAGGAGGTGAAGATGCGTATTACCATAAAAAGATGGGGGTACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTAAATCCGGGGCAGAGAGTTGAGGCTCAAGTGACca                  >  1:169998/7‑150 (MQ=255)
                                                                                                                                      ccttaaGTTATACAAAGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGcaa                 >  2:1697982/6‑150 (MQ=255)
                                                                                                                                      ccGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGGATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGGGGAGGCGCAAGTGAGcaa                 <  2:77620/150‑1 (MQ=255)
                                                                                                                                       cttaaGTTATACAAAGTATAACTTAAGGAGTTGCAGATGCGTATTACAATAAAAAGATGGGTTAACAGTGCAGTTATGGTCATTCCCAATATCGTAAAGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGcaac                >  2:37306/5‑150 (MQ=255)
                                                                                                                                       cGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGcaac                <  2:774980/150‑1 (MQ=255)
                                                                                                                                        gTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGcaaca               <  2:656645/150‑1 (MQ=255)
                                                                                                                                         tGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGcaacaa              <  2:682262/150‑1 (MQ=255)
                                                                                                                                           aGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATc            >  2:255410/1‑150 (MQ=255)
                                                                                                                                            gTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCa           >  2:1000070/1‑150 (MQ=255)
                                                                                                                                                     tGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCAACTGATTCt  <  2:1299985/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGACGACTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCAACTGATTCT  >  NC_000913/4448276‑4448574

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