Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 729,376 G→T 6.1% intergenic (+36/‑207) 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,3760GT 6.1% 757.8 / 27.5 212intergenic (+36/‑207)ybfA/rhsCDUF2517 family protein/Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base G (114/84);  new base T (2/11);  total (116/96)
Fisher's exact test for biased strand distribution p-value = 3.55e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

Only 100 of 215 total aligned reads displayed.

CTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATA  >  NC_000913/729229‑729524
                                                                                                                                                   |                                                                                                                                                    
cTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGaa                                                                                                                                                    <  2:2970503/150‑1 (MQ=255)
  ggAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAAc                                                                                                                                                  >  1:3244069/1‑150 (MQ=255)
    aaaGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGc                                                                                                                                                >  1:3533333/1‑150 (MQ=255)
     aaGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCt                                                                                                                                               <  1:2551776/150‑1 (MQ=255)
        aCCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGc                                                                                                                                            >  1:928743/1‑150 (MQ=255)
          cgcgCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCgg                                                                                                                                          <  2:110972/150‑1 (MQ=255)
           gcgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTTATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGt                                                                                                                                         >  1:2063405/1‑150 (MQ=255)
           gcgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGt                                                                                                                                         >  1:2407272/1‑150 (MQ=255)
           gcgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGt                                                                                                                                         >  1:4442635/1‑150 (MQ=255)
            cgcCCGCTTTTACAGCTACCTGCATCGCTTCTGGTCGAAAACCAGCGATTAACCGGTGTGGATGGATCAGGCCGAATAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGtt                                                                                                                                        <  2:4006212/150‑1 (MQ=255)
            cgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGtt                                                                                                                                        <  1:2984837/150‑1 (MQ=255)
            cgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGtt                                                                                                                                        <  2:4006173/150‑1 (MQ=255)
            cgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtt                                                                                                                                        <  2:2715828/150‑1 (MQ=255)
             gcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGttt                                                                                                                                       <  2:1624837/150‑1 (MQ=255)
             gcCCGCTATTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGAGTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGttt                                                                                                                                       <  2:3243793/150‑1 (MQ=255)
                cGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCTATAAACCGGTGTGGTTGGATCAGGCCGAAAAGGCATCTGGTGTTTTTTATTGTTTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtttttt                                                                                                                                    <  1:780975/150‑1 (MQ=255)
                cGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCTATAAACCGGTGTGGATGTATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTTAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtttttt                                                                                                                                    <  1:646235/150‑1 (MQ=255)
                cGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTTTTGATTGCTTAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtttttt                                                                                                                                    <  1:2991755/150‑1 (MQ=255)
                cGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtttttt                                                                                                                                    <  2:2289791/150‑1 (MQ=255)
                      tACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTAATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCa                                                                                                                              <  2:3706853/150‑1 (MQ=255)
                       aCAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAg                                                                                                                             >  2:2994077/1‑150 (MQ=255)
                        cAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGCGTGGTTGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGATTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAgg                                                                                                                            <  1:3582711/150‑1 (MQ=255)
                           cTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAg                                                                                                                         >  2:1573462/1‑150 (MQ=255)
                           cTACCGGCATCGCGTCTGGTCGAAACCCAGCGATAAACCGGGGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAg                                                                                                                         <  2:4121155/150‑1 (MQ=255)
                               cTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGt                                                                                                                     >  1:2632882/1‑150 (MQ=255)
                                 gCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTc                                                                                                                   >  2:2470207/1‑150 (MQ=255)
                                    tCGCGTCTTGTAGAAAACCAGCTATAAACCGTTGTTGATTTATCTGGCCGAAATGTAAAATTGTGATTTTTTTTGTTTGCTGATTTTTTAATACTTTGAAAACCGCAAGCGGAAACGGTGGCGGTTTTTTGTATCAGGAAGATGTTCaga                                                                                                                <  1:1489335/150‑1 (MQ=255)
                                    tCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCaga                                                                                                                <  1:2798264/150‑1 (MQ=255)
                                    tCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCaga                                                                                                                >  1:4031924/1‑150 (MQ=255)
                                            ggTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTATATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATAtt                                                                                                        <  2:3798556/150‑1 (MQ=255)
                                             gTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATAttt                                                                                                       >  2:1664330/1‑150 (MQ=255)
                                              tCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTa                                                                                                      >  2:1498935/1‑150 (MQ=255)
                                               cGAAAACCAGCGATTAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAg                                                                                                     <  1:3506714/150‑1 (MQ=255)
                                               cGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGTTTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAg                                                                                                     <  1:2144665/150‑1 (MQ=255)
                                               cGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAg                                                                                                     <  1:1592762/150‑1 (MQ=255)
                                                gAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGa                                                                                                    >  2:3454565/1‑150 (MQ=255)
                                                     ccAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGAAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCg                                                                                               <  1:3646894/150‑1 (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:4478180/1‑150 (MQ=255)
                                                       aGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGtt                                                                                             >  1:3758634/1‑150 (MQ=255)
                                                       aGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGACAGCGGAAACGCTGGCGCTTTTTTGTATCATGAAGATGTTCAGAGGAATATTTAGAAATAGtt                                                                                             >  1:2343625/1‑150 (MQ=255)
                                                        gCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTta                                                                                            <  2:2036272/150‑1 (MQ=255)
                                                          gATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTtata                                                                                          >  1:3038295/1‑150 (MQ=255)
                                                            tAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATAcc                                                                                        >  2:3075855/1‑150 (MQ=255)
                                                            tAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATAcc                                                                                        >  2:1985655/1‑150 (MQ=255)
                                                                 cGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCg                                                                                   >  2:697028/1‑150 (MQ=255)
                                                                     gtgGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAAc                                                                               >  2:4153734/1‑150 (MQ=255)
                                                                     gtgGAAGGATCAGGCCGAAAAGGCAACTGGTGATATTTAATGATTGCTGAGTTATGAATACATAGAAAACCGCCAGCGGAAACCCTGGCGCAATATAGAATCACGAAGATGTTCACAGGAATATATAGACATCGTTATAACAATCGAAAc                                                                               >  2:3658268/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)
                                                                           ggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatatat                                                                         >  1:2916325/1‑150 (MQ=255)
                                                                               cAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCAt                                                                     <  1:3356582/150‑1 (MQ=255)
                                                                               cAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCAt                                                                     <  1:4544900/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:545920/1‑150 (MQ=255)
                                                                                   ccGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATAACAATCGAAACATATATTAATGaaa                                                                 >  1:2218695/1‑150 (MQ=255)
                                                                                    cGAAAACCCAACTCGAGATATTTAAAGATAGCTCAGAAATGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGATTTATCTAACAGGAAGATGTTCACACGAATATTTAGACATCGATAAACCAATCGAAACATAAAATCATCaaaa                                                                >  2:1474213/1‑149 (MQ=255)
                                                                                          ggCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGACACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAatatatat                                                          >  1:1437339/1‑150 (MQ=255)
                                                                                             aaCTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATaaa                                                       >  2:280125/1‑150 (MQ=255)
                                                                                             aaCTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAAAAAATATTCATGAAATATATATaaa                                                       >  1:1301455/1‑150 (MQ=255)
                                                                                                    gATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTc                                                <  2:141038/150‑1 (MQ=255)
                                                                                                    gATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTATTTTGTATCAGGACGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTc                                                <  2:518289/150‑1 (MQ=255)
                                                                                                    gATTTAGATTGATTGCTGAGTTTTGAATCCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTc                                                <  2:4205533/150‑1 (MQ=255)
                                                                                                     aTTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCACCGGAAACCCAGGCGGTTTTTTGTATCAGGAAGAGGTTCACAGGAATATTTAGAAATCGTTATACCAAACGAAACATATAATCATGAAAAAAATATAAAAAttttac                                               >  1:197638/1‑148 (MQ=255)
                                                                                                       ttttATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTa                                             >  2:2131899/1‑150 (MQ=255)
                                                                                                        tttATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTaa                                            >  2:4258617/1‑150 (MQ=255)
                                                                                                        tttATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAAAAGTTCAGAGGAATAATTAGACAACCTTATACCAATCGAAACAAAAATTCAAGAAAAATATATAAATATTTTCCTaa                                            >  1:3452336/1‑150 (MQ=255)
                                                                                                        tttATTGATTGCTGAGTTATGAATACATTGAAAACCGCCAGCCGAAACGCTGGCGGTTTTTTGTATCAGGAAGAAGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAACTATATATAAATATTTTCCTaa                                            >  2:2541618/1‑150 (MQ=255)
                                                                                                         ttATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTCTCCTAAt                                           <  1:2432383/150‑1 (MQ=255)
                                                                                                           attgattgCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATAAGGAATATGTTCAGTGGAATATTTAGACAACGGTATACCAATCGAAACAAATATTAATGAAATATATATAAATATTTTAATAAttt                                         >  1:444073/1‑149 (MQ=255)
                                                                                                              gattgCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATAAGGAAGATGTTCAGAGGAATATTTAGACAACGTTATACAAAACCAAACAAAAATTCAATAAATAAATAAAAAAATTTTACTAAATTTTc                                      >  1:471452/1‑150 (MQ=255)
                                                                                                                ttgCTGAGTTTTGAATACATTGAAAACCTCCAGCGGAAACGCTGGCTGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCtt                                    <  2:2876482/150‑1 (MQ=255)
                                                                                                                 tgCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATAAGGAAGTTGTTAAGAGGAATATTTAGACATCGTTATACCAATCGAAAAATATATTCATGAAATATATATAAAAATTTTCCTAATTGTTCTTa                                   >  1:2651553/1‑150 (MQ=255)
                                                                                                                  gCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTGATTGTTCTTAt                                  <  2:2073745/150‑1 (MQ=255)
                                                                                                                  gCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTGATTGTTCTTAt                                  <  2:569900/150‑1 (MQ=255)
                                                                                                                  gCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTAt                                  <  2:3554132/150‑1 (MQ=255)
                                                                                                                    tGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTTTATGAGGAAGATGTTCAGTGGAATATTTAGACAAAGGTATACCAAACGAAACATATATTCAAGAAATAAAAATTAATATTTTAAAAATTGTTATTATCt                                >  1:1599205/1‑150 (MQ=255)
                                                                                                                    tGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCt                                >  2:3610547/1‑150 (MQ=255)
                                                                                                                     gAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTg                               >  2:4240575/1‑150 (MQ=255)
                                                                                                                     gAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTg                               >  2:1856688/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)
                                                                                                                                       aaaCCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAg             <  2:3662508/150‑1 (MQ=255)
                                                                                                                                            cccAGCGCAAACGCTGGCGGATATTTGTATCAGGAAGATGTTCAGAGCAATATTTAGACAACGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCttt        <  2:1336019/149‑1 (MQ=255)
                                                                                                                                               aGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCtt     >  2:3990558/1‑150 (MQ=255)
                                                                                                                                                gCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTta    >  2:189853/1‑150 (MQ=255)
                                                                                                                                                 cGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTtat   <  2:4450713/150‑1 (MQ=255)
                                                                                                                                                  ggAAACGCTGGCGGTTTTTTGTATCAGGAAGAGGTTCAGAGGTCTATTTAGTCTTCGTTATACCAATCGAAACATATATTCATGAATTATATATAAATATTTTCCTATTTGTTCTTATCGGACAGATATCTCACTTAAGGCTTTCTtata  <  2:3422200/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATA  >  NC_000913/729229‑729524

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: