Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 729,377 A→T 6.5% intergenic (+37/‑206) ybfA → / → rhsC DUF2517 family protein/Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913729,3770AT 6.5% 759.9 / 27.9 212intergenic (+37/‑206)ybfA/rhsCDUF2517 family protein/Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (113/85);  new base T (3/11);  total (116/96)
Fisher's exact test for biased strand distribution p-value = 1.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 216 total aligned reads displayed.

CTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAA  >  NC_000913/729229‑729526
                                                                                                                                                    |                                                                                                                                                     
cTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGaa                                                                                                                                                      <  2:2970503/150‑1 (MQ=255)
  ggAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAAc                                                                                                                                                    >  1:3244069/1‑150 (MQ=255)
  ggAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAAc                                                                                                                                                    <  1:2305187/150‑1 (MQ=255)
     aaGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCt                                                                                                                                                 <  1:2551776/150‑1 (MQ=255)
       gACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTgg                                                                                                                                               <  2:1360232/150‑1 (MQ=255)
          cgcgCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCgg                                                                                                                                            <  2:110972/150‑1 (MQ=255)
           gcgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTTATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGt                                                                                                                                           >  1:2063405/1‑150 (MQ=255)
           gcgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGt                                                                                                                                           >  1:4442526/1‑150 (MQ=255)
           gcgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGt                                                                                                                                           >  1:2407272/1‑150 (MQ=255)
           gcgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGt                                                                                                                                           <  2:2846766/150‑1 (MQ=255)
            cgcCCGCTTTTACAGCTACCTGCATCGCTTCTGGTCGAAAACCAGCGATTAACCGGTGTGGATGGATCAGGCCGAATAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGtt                                                                                                                                          <  2:4006212/150‑1 (MQ=255)
            cgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGtt                                                                                                                                          <  2:4006173/150‑1 (MQ=255)
            cgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGtt                                                                                                                                          <  1:2985052/150‑1 (MQ=255)
            cgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtt                                                                                                                                          <  2:2715828/150‑1 (MQ=255)
             gcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGttt                                                                                                                                         <  2:1624837/150‑1 (MQ=255)
              cccGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtttt                                                                                                                                        >  1:3433429/1‑150 (MQ=255)
                ggCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCTATAAACCGGTGTGGTTGTATCAGGCCGAAATGGCAACTGGTTATTTTTATTGTTTGCTTATTTTTTAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtttttt                                                                                                                                      <  1:2992160/149‑1 (MQ=255)
                 gCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTg                                                                                                                                     >  2:2776402/1‑150 (MQ=255)
                      tACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTAATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCa                                                                                                                                <  2:3706853/150‑1 (MQ=255)
                       aCAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAg                                                                                                                               >  2:2994077/1‑150 (MQ=255)
                           cTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAg                                                                                                                           >  2:1573462/1‑150 (MQ=255)
                           cTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTATTTTGTATCAGGAAg                                                                                                                           <  2:3164126/150‑1 (MQ=255)
                              ccTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTTTCAGGTAGATg                                                                                                                        >  1:4027961/1‑150 (MQ=255)
                                 gCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTc                                                                                                                     >  2:2470207/1‑150 (MQ=255)
                                    tCGCGTCTTGTAGAAAACCAGCTATAAACCGTTGTTGATTTATCTGGCCGAAATGTAAAATTGTGATTTTTTTTGTTTGCTGATTTTTTAATACTTTGAAAACCGCAAGCGGAAACGGTGGCGGTTTTTTGTATCAGGAAGATGTTCaga                                                                                                                  <  1:1489335/150‑1 (MQ=255)
                                    tCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCTAAAAGGCAACTGGTTATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCaga                                                                                                                  <  1:2798607/150‑1 (MQ=255)
                                    tCGCGTCGGGTCTAAAACCCGCGATAAACCGGTGTGGAGGGATCAGGCCGAATTGGCAACTGGTTATTTTTTTTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCaga                                                                                                                  <  1:4459366/150‑1 (MQ=255)
                                            ggTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTATATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATAtt                                                                                                          <  2:3798556/150‑1 (MQ=255)
                                             gTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATAttt                                                                                                         >  1:3554132/1‑150 (MQ=255)
                                             gTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATAttt                                                                                                         >  1:474258/1‑150 (MQ=255)
                                             gTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATAttt                                                                                                         >  2:1664288/1‑150 (MQ=255)
                                              tCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTa                                                                                                        >  2:1498935/1‑150 (MQ=255)
                                              tCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTa                                                                                                        >  2:1678942/1‑150 (MQ=255)
                                               cGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGTTTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAg                                                                                                       <  1:2144665/150‑1 (MQ=255)
                                               cGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAg                                                                                                       <  1:1592762/150‑1 (MQ=255)
                                                gAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGa                                                                                                      >  2:3454565/1‑150 (MQ=255)
                                                      cAGCGATAAACCGGTGTGGATGGATCAGTCCGAAAAGGCAACTGGTGATATTAATAGATTGCTGAGATTTGAATACATTGAAATCCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCgg                                                                                                >  2:4211940/1‑149 (MQ=255)
                                                      cAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGACTATTTAGACATCGt                                                                                                >  1:1237974/1‑150 (MQ=255)
                                                      cAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGt                                                                                                >  2:1184588/1‑150 (MQ=255)
                                                      cAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATAAGGAAGATGTTCAGAGGAATATTTAGACATCgg                                                                                                >  1:3811380/1‑149 (MQ=255)
                                                       aGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGtt                                                                                               >  1:3758634/1‑150 (MQ=255)
                                                       aGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGACAGCGGAAACGCTGGCGCTTTTTTGTATCATGAAGATGTTCAGAGGAATATTTAGAAATAGtt                                                                                               >  1:2343625/1‑150 (MQ=255)
                                                        gCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTta                                                                                              <  2:2036272/150‑1 (MQ=255)
                                                            tAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATAcc                                                                                          >  2:1985655/1‑150 (MQ=255)
                                                            tAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATAcc                                                                                          >  2:3075855/1‑150 (MQ=255)
                                                              aaCCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCaa                                                                                        >  2:2104308/1‑150 (MQ=255)
                                                              aaCCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCaa                                                                                        >  2:2286835/1‑150 (MQ=255)
                                                                 cGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCg                                                                                     >  2:3654789/1‑150 (MQ=255)
                                                                 cGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCg                                                                                     >  2:697028/1‑150 (MQ=255)
                                                                    tgtgGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaaa                                                                                  <  1:3591275/150‑1 (MQ=255)
                                                                    tgtgGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaaa                                                                                  <  1:650778/150‑1 (MQ=255)
                                                                     gtgGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAAc                                                                                 >  2:4153734/1‑150 (MQ=255)
                                                                      tggatggaTCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGAAATCGGTATACCAAACGAAACa                                                                                >  1:1832080/1‑150 (MQ=255)
                                                                       ggatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACat                                                                               <  2:2760631/150‑1 (MQ=255)
                                                                           ggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatatat                                                                           >  2:2905348/1‑150 (MQ=255)
                                                                               cAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCAt                                                                       <  1:4544900/150‑1 (MQ=255)
                                                                               cAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCAt                                                                       <  1:3356582/150‑1 (MQ=255)
                                                                                 ggCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGa                                                                     <  1:1821797/150‑1 (MQ=255)
                                                                                  gCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGaa                                                                    >  2:854198/1‑150 (MQ=255)
                                                                                  gCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGaa                                                                    >  1:1921511/1‑150 (MQ=255)
                                                                                   ccGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGaaa                                                                   <  2:1969947/150‑1 (MQ=255)
                                                                                   ccGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGaaa                                                                   <  1:2960434/150‑1 (MQ=255)
                                                                                    cGAAAAGGCAACAGGTGATTTTTATTGATTCCTGAGTTTTGAAAACAAAGAAAACCGCCAGCGGAAACGCTGGCGCTTATATGTATCAGGAAGATGTTCAGACCAAAAAATAGACATCGTTATACCAATCGAAACATATATTCATGaaaa                                                                  >  2:1474179/1‑149 (MQ=255)
                                                                                       aaaGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAatata                                                               <  1:3489132/150‑1 (MQ=255)
                                                                                          ggCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGACACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAatatatat                                                            >  1:1437339/1‑150 (MQ=255)
                                                                                             aaCTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAAAAAATATTCATGAAATATATATaaa                                                         >  1:1301455/1‑150 (MQ=255)
                                                                                                tGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACCTCGTTATACCAATCGAAACATATATTAATGAAATATATATAAatat                                                      >  1:602586/1‑150 (MQ=255)
                                                                                                    gATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTATTTTGTATCAGGACGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTc                                                  <  2:518289/150‑1 (MQ=255)
                                                                                                    gATTTAGATTGATTGCTGAGTTTTGAATCCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTc                                                  <  2:4205533/150‑1 (MQ=255)
                                                                                                     aTTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCACCGGAAACCCAGGCGGTTTTTTGTATCAGGAAGAGGTTCACAGGAATATTTAGAAATCGTTATACCAAACGAAACATATAATCATGAAAAAAATATAAAAAttttac                                                 >  1:197638/1‑148 (MQ=255)
                                                                                                        tttATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTTTATAAGGAAGATGTTAAGAGGAAAATTTAGAAAAAGTTATAACAAAAGAAAAAAATAATCATTAAATATATATAAATATTTTACTaa                                              >  1:4058934/1‑150 (MQ=255)
                                                                                                        tttATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTaa                                              >  2:4258617/1‑150 (MQ=255)
                                                                                                        tttATTGATTGCTGAGTTATGAATACATTGAAAACCGCCAGCCGAAACGCTGGCGGTTTTTTGTATCAGGAAGAAGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAACTATATATAAATATTTTCCTaa                                              >  2:2541618/1‑150 (MQ=255)
                                                                                                         ttATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTCTCCTAAt                                             <  1:2432383/150‑1 (MQ=255)
                                                                                                          tATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAAtt                                            <  2:2063405/150‑1 (MQ=255)
                                                                                                              gattgCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTc                                        <  2:4031924/150‑1 (MQ=255)
                                                                                                              gattgCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATAAGGAAGATGTTCAGAGGAATATTTAGACAACGTTATACAAAACCAAACAAAAATTCAATAAATAAATAAAAAAATTTTACTAAATTTTc                                        >  1:471452/1‑150 (MQ=255)
                                                                                                                 tgCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGTCGGTTTTTTGTATCAGGAAGATGTTAAGAGGAATATTTAGACATCTTTATAACAATCCAAAAAAATATACATGAAAAATATATAAATATATTACTAATTGttatta                                     >  1:2651574/1‑150 (MQ=255)
                                                                                                                 tgCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATAAGGAAGTTGTTAAGAGGAATATTTAGACATCGTTATACCAATCGAAAAATATATTCATGAAATATATATAAAAATTTTCCTAATTGTTCTTa                                     >  1:2651553/1‑150 (MQ=255)
                                                                                                                  gCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTGATTGTTCTTAt                                    <  2:569900/150‑1 (MQ=255)
                                                                                                                  gCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTAt                                    <  2:3554132/150‑1 (MQ=255)
                                                                                                                   cTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTAACAGGAAGAGGTTCAGAGGAAAAGTTAGACATCCTGATACCAAACCAAACAAATATTCATGAAATAAATATAAATATTTTACAAATTGTTCTTATc                                   >  1:3515137/1‑150 (MQ=255)
                                                                                                                    tGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTTTATGAGGAAGATGTTCAGTGGAATATTTAGACAAAGGTATACCAAACGAAACATATATTCAAGAAATAAAAATTAATATTTTAAAAATTGTTATTATCt                                  >  1:1599205/1‑150 (MQ=255)
                                                                                                                    tGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCt                                  >  2:3610547/1‑150 (MQ=255)
                                                                                                                     gAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTg                                 >  2:4240575/1‑150 (MQ=255)
                                                                                                                       gTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGAc                               <  2:4519397/150‑1 (MQ=255)
                                                                                                                        ttttGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTTTATAAGGAAGATGTTAAGAGGAAAATTTATACAACGTTGTACCAATAGAAACAAATAATAAAGAAAAATATATAAATATTTTACTAATTTTtattatattaaa                              >  1:155666/1‑144 (MQ=255)
                                                                                                                         tttGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAg                             >  2:3245084/1‑150 (MQ=255)
                                                                                                                          ttGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGa                            <  2:3125781/150‑1 (MQ=255)
                                                                                                                             aaTACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGatat                         <  2:432180/150‑1 (MQ=255)
                                                                                                                               tACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATAtct                       >  2:122540/1‑150 (MQ=255)
                                                                                                                                   ttGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACt                   >  2:1579864/1‑150 (MQ=255)
                                                                                                                                    tGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACtt                  >  2:1314495/1‑150 (MQ=255)
                                                                                                                                      aaaaCCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTaa                <  1:792711/150‑1 (MQ=255)
                                                                                                                                       aaaCCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAg               <  2:3662441/150‑1 (MQ=255)
                                                                                                                                               aGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCtt       >  2:3990558/1‑150 (MQ=255)
                                                                                                                                                gCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTta      >  2:1892164/1‑150 (MQ=255)
                                                                                                                                                 cGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTtat     <  2:4450713/150‑1 (MQ=255)
                                                                                                                                                  ggAAACGCTGGCGGTTTTTTGTATCAGGAAGAGGTTCAGAGGTCTATTTAGTCTTCGTTATACCAATCGAAACATATATTCATGAATTATATATAAATATTTTCCTATTTGTTCTTATCGGACAGATATCTCACTTAAGGCTTTCTtata    <  2:3422200/150‑1 (MQ=255)
                                                                                                                                                    aaaCGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTATACAACGTTATACCAATCGAAACAAATATTCATGAAAAATATATAAAAATTTTACAAATTGTTATTATCTGACAGATATATCACTTAATGCTTTATTATaaa  >  1:2630234/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAA  >  NC_000913/729229‑729526

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