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

Only 100 of 130 total aligned reads displayed.

GCGAGACGACTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCAAC  >  NC_000913/4448274‑4448567
                                                                                                                                                     |                                                                                                                                                
 cGAGACGACTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCTCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACtt                                                                                                                                                 <  1:52131/150‑1 (MQ=255)
  gagaCGACTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACttt                                                                                                                                                <  2:1623027/150‑1 (MQ=255)
      cgacTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGtat                                                                                                                                            >  1:891969/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)
               ttCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACAAAGTATAACtcacgg                                                                                                                                   <  1:1816475/150‑6 (MQ=255)
               ttCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACAAAGTATAACtcacgg                                                                                                                                   <  2:1281667/150‑6 (MQ=255)
                tCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAagga                                                                                                                                  <  1:450950/150‑1 (MQ=255)
                 cTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAACACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGATTTATACTTTGTATAACTTAaggag                                                                                                                                 >  2:1089607/1‑150 (MQ=255)
                 cTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAagtat                                                                                                                                 >  1:1563850/1‑147 (MQ=255)
                  tGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAaggagg                                                                                                                                <  1:333032/150‑1 (MQ=255)
                   gTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCGTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGt                                                                                                                               >  1:1788856/1‑150 (MQ=255)
                   gTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGt                                                                                                                               >  1:1392713/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)
                         tAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGAGCAGGGGTAAAACGTTTCCTGAAGCCCCGTGAGATATACGTTGTATAACTTAAGGAGGTGCAGAt                                                                                                                         >  1:285896/1‑150 (MQ=255)
                          aaTTCCCTCCCCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATg                                                                                                                        >  1:1625655/1‑150 (MQ=255)
                           aTTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCATGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGc                                                                                                                       >  1:764168/1‑150 (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:1838604/1‑150 (MQ=255)
                                  cACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTCCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGAATAACTAAAGGAGGTGCAGATGCGTATTAc                                                                                                                >  1:1186536/1‑150 (MQ=255)
                                   aCCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTAcc                                                                                                               <  1:779226/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)
                                                  ctcCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAAAGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGAtggtgg                                                                                                >  1:948145/1‑150 (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)
                                                           gCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATTGGGGAACAGTGCa                                                                                       >  1:1696684/1‑150 (MQ=255)
                                                            cTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAg                                                                                      <  1:1567945/150‑1 (MQ=255)
                                                             tttCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAgg                                                                                     <  2:457890/150‑1 (MQ=255)
                                                                ccATATCTATTAACGAATAAAAAACTCTGCTGGCATTCACAAAAGCGCAGGGGTAAAACGTTTCGTGTAGCGCCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGAGTATTACAATAAAAAGATTGTGGAACTGTGCAGGTAt                                                                                  >  2:1811894/1‑150 (MQ=255)
                                                                     tCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCa                                                                             >  2:1423948/1‑150 (MQ=255)
                                                                     tCTATTAACGCATAAAAAACTCTCCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCa                                                                             <  1:220014/150‑1 (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:1177773/150‑1 (MQ=255)
                                                                                               ggCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGaactt                                                   >  1:1500464/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)
                                                                                                             cgcAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATCCTTTGTTTAACTTAAGGAGGGGCAGATGCGTATTAGGATAAAAAGATGGGGGGACAGTGCAGGTATGGTCATTCCCAATATCGTATTGAGAGATCTTACCTTACAGCCggg                                     >  2:1096702/1‑150 (MQ=255)
                                                                                                             cgcAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCggg                                     >  2:642050/1‑150 (MQ=255)
                                                                                                             cgcAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCggg                                     >  2:394192/1‑150 (MQ=255)
                                                                                                               cAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGc                                   <  2:1482409/150‑1 (MQ=255)
                                                                                                               cAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGc                                   <  1:1264779/150‑1 (MQ=255)
                                                                                                               cAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGc                                   <  2:356162/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)
                                                                                                                          cGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAg                        >  2:1721695/1‑150 (MQ=255)
                                                                                                                            tatCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGgc                      <  2:517902/148‑1 (MQ=255)
                                                                                                                               ccTGTAGCAACGTGAGTTATACTTTTTAAAACTTAAGGAGGGGCAGATGCGTAATACCATAAAAAGATGGGGGAACAGTGCAGGTATGGGCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGGGGAGGCGCa                   <  1:1748646/150‑1 (MQ=255)
                                                                                                                                cTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCaa                  >  1:842471/1‑150 (MQ=255)
                                                                                                                                 tGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAg                 >  2:794693/1‑150 (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)
                                                                                                                                       tccttaaGTTATACAAAGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAATGATGGGGGAACAGTGAAGGTTTGGTCATTCCCATTATCGTAATGAAAGAAGTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGca           >  1:340131/7‑150 (MQ=255)
                                                                                                                                       tccttaaGTTATACAAAGTATAACTTAAGGAGGTGAAGATGCGTATTACCATAAAAAGATGGGGGTACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTAAATCCGGGGCAGAGAGTTGAGGCTCAAGTGACca           >  1:169998/7‑150 (MQ=255)
                                                                                                                                       aCCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGca           >  1:521109/1‑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)
                                                                                                                                                     |                                                                                                                                                
GCGAGACGACTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCAAC  >  NC_000913/4448274‑4448567

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: