Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,372 T→C 5.4% intergenic (+33/‑483) sseA → / → ryfA 3‑mercaptopyruvate sulfurtransferase/novel sRNA, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,653,3720TC 5.4% 343.5 / 6.5 111intergenic (+33/‑483)sseA/ryfA3‑mercaptopyruvate sulfurtransferase/novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base T (52/53);  new base C (4/2);  total (56/55)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 112 total aligned reads displayed.

CAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCG  >  NC_000913/2653224‑2653521
                                                                                                                                                    |                                                                                                                                                     
cAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTg                                                                                                                                                      >  1:1430022/1‑150 (MQ=255)
 aGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTgg                                                                                                                                                     <  1:801949/150‑1 (MQ=255)
 aGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTgg                                                                                                                                                     <  1:431634/150‑1 (MQ=255)
   ccGTGGTTTTGTTAGCACTCGCGACGCTGGAAGTGCCAAACGTGAAACTGTACGACCGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCg                                                                                                                                                   <  1:1694544/150‑1 (MQ=255)
        gTTTTCTTAGCACTCGCGACGCTGGACGTGCCAAACGTGAAACTGTACGACGGCGCATGGATTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGtttt                                                                                                                                              <  1:1149509/150‑1 (MQ=255)
         ttttGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGttttt                                                                                                                                             <  2:1793470/150‑1 (MQ=255)
              ttAGCACTCGCGACGCTGGATGTGCCAAACGGGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTg                                                                                                                                        <  1:388032/150‑1 (MQ=255)
                aGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCt                                                                                                                                      >  2:1864525/1‑150 (MQ=255)
                 gCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTa                                                                                                                                     >  1:1680385/1‑150 (MQ=255)
                 gCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTa                                                                                                                                     >  1:1678175/1‑150 (MQ=255)
                    cTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGGTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTc                                                                                                                                  <  2:1324299/150‑1 (MQ=255)
                    cTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGAGAATGGGGCGCGCGGGCAGATTTACCGGATGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTATTTTATTGCTAGtg                                                                                                                                  >  1:1345736/1‑149 (MQ=255)
                      cgcgACGCTGGATGTGCCAAACGTGAAACTGTACGCCGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTg                                                                                                                                >  2:1549434/1‑150 (MQ=255)
                      cgcgACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTg                                                                                                                                >  1:868691/1‑150 (MQ=255)
                      cgcgACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAAGGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTg                                                                                                                                <  2:1319237/150‑1 (MQ=255)
                       gcgACGCTGGATGTGCCATACGTGAAACTGTACGACTTCGCATGGAGTGAATGGGGCGCGCGGGCAGAATTACCGGGTGAGCCTGTGAAATAAGTATTTTACAGGCAATAAAAACCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTgg                                                                                                                               <  1:201349/150‑1 (MQ=255)
                       gcgACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTgg                                                                                                                               <  2:1762248/150‑1 (MQ=255)
                         gACGCTGGAGGTGCGAAACGTGAAACTGTACGACGGCGGATGGTGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGtt                                                                                                                             <  2:1019399/150‑1 (MQ=255)
                           cGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAAGTTGGCGGGTTTTTATTGCTAGTCTGGTTcg                                                                                                                           >  2:1835380/1‑150 (MQ=255)
                            gCTGGAGGTGCCAAACGTGAAACTGTACGACGGCGCAAGGAGTGACTGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTGTACAGGCCATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTcgc                                                                                                                          <  1:213234/150‑1 (MQ=255)
                             cTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTcgcg                                                                                                                         >  1:224761/1‑150 (MQ=255)
                               ggATGTGCCAACCGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGc                                                                                                                       >  1:963564/1‑150 (MQ=255)
                                    tgCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTc                                                                                                                  >  2:480269/1‑150 (MQ=255)
                                     gCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTcc                                                                                                                 >  1:522742/1‑150 (MQ=255)
                                      ccAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGCGCCAGTGAAATAAGTATTATACAGGCAATAAAAAACCGCCGAAATTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCca                                                                                                                >  2:1764370/1‑150 (MQ=255)
                                       cAAACGTGAAACTGTACGACGGCGCATGGTGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcag                                                                                                               >  2:1706360/1‑150 (MQ=255)
                                                aaCTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGAc                                                                                                      >  1:1621743/1‑150 (MQ=255)
                                                aaCTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGAc                                                                                                      <  1:1166593/150‑1 (MQ=255)
                                                    gTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgt                                                                                                  >  1:1513064/1‑150 (MQ=255)
                                                    gTACGACGGCGCATGGAGGGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgt                                                                                                  <  2:313204/150‑1 (MQ=255)
                                                       cgacgGCTCATGTAGTGAATGGGGCGCGCGGGCAGATTTACCGGTAGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGGTCGCGGCCTTTCCAGCAGGTTGACTTGGGtt                                                                                               <  2:1634023/150‑1 (MQ=39)
                                                         acgTCGAATGGAGTGATTGGGGCGCTGGGGGAGATTTACCGGTTGATCCAGTGAAATAAGTATTTTACAGGCAATATTAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGAGTCCTTTCCAGCAGGTGGACTTGTGTTAc                                                                                             <  2:841190/150‑1 (MQ=255)
                                                           ggCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAATAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACAt                                                                                           <  1:1056147/150‑1 (MQ=255)
                                                               cATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGc                                                                                       >  1:1501401/1‑150 (MQ=255)
                                                                      tGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAg                                                                                <  1:176692/150‑1 (MQ=255)
                                                                      tGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAg                                                                                <  2:1486868/150‑1 (MQ=255)
                                                                       gAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTTGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAgg                                                                               >  2:1144467/1‑150 (MQ=255)
                                                                           gggtCGCGCGGTCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACATGCAATAAAAAACCGCCGAATTTGGCTGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCt                                                                           <  2:1612842/150‑1 (MQ=255)
                                                                             ggCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTc                                                                         >  2:782881/1‑150 (MQ=255)
                                                                             ggCGCGCGGGCAGATTTACCGGGTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTGCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTc                                                                         <  1:428706/150‑1 (MQ=255)
                                                                               cgcgcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTcac                                                                       <  2:1577961/150‑1 (MQ=255)
                                                                               cgcgcgGGCAGATTTAACTGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAATCCGCAGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTcac                                                                       <  2:69340/150‑1 (MQ=255)
                                                                                  gcgGGCAGATCTACCGGGTGAGCCAGTGAAATAATGATTTTCCACGCAATAAAAAACCGCCGAACTTGGCGGTTTTTTATTGCTAGTCTGGTGCGCGGCCTTGCCAGCAGGTTGACTTGTGTTACCTGAGCAACGCAGGTGCTTCACAGc                                                                    <  1:194298/150‑1 (MQ=255)
                                                                                     ggCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCaaa                                                                 >  2:1735717/1‑150 (MQ=255)
                                                                                     ggCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCaaa                                                                 <  2:1227380/150‑1 (MQ=255)
                                                                                       cAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTTGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAAc                                                               >  1:1741638/1‑150 (MQ=255)
                                                                                       cAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAAc                                                               >  1:892159/1‑150 (MQ=255)
                                                                                       cAGATTTACCGGGTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAAc                                                               <  1:1730494/150‑1 (MQ=255)
                                                                                       cAGATTTACCGGATGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAAc                                                               <  1:1006776/150‑1 (MQ=255)
                                                                                         gATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACaa                                                             >  1:1502560/1‑150 (MQ=255)
                                                                                          aTTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACaaa                                                            >  2:1127273/1‑149 (MQ=255)
                                                                                           tttACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATa                                                           <  2:24985/150‑1 (MQ=255)
                                                                                            ttACCGGTTGATCCATTGAAATAACTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATAc                                                          <  1:201168/150‑1 (MQ=255)
                                                                                            ttACCGGTTGAGGAAGTGAAATAAGTTTTTTAAAGGCAATAAAAAACCGCCGAATTTGGCTGTTTTTTATTGCTAGGCTGGTGCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATAc                                                          <  2:100146/150‑1 (MQ=255)
                                                                                            ttACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATAc                                                          <  2:815857/150‑1 (MQ=255)
                                                                                               ccGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCa                                                       >  2:801522/1‑150 (MQ=255)
                                                                                                cGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCAc                                                      >  2:674888/1‑150 (MQ=255)
                                                                                                 ggTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCAcc                                                     <  1:1828781/150‑1 (MQ=255)
                                                                                                    tGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGt                                                  >  2:1046275/1‑150 (MQ=255)
                                                                                                     gAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTa                                                 >  2:1846342/1‑150 (MQ=255)
                                                                                                       gCCAGTGAGATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAAc                                               >  1:1634182/1‑150 (MQ=255)
                                                                                                        ccAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTAACCAGTAAct                                              >  1:1841554/1‑150 (MQ=255)
                                                                                                          aGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAActct                                            >  2:521899/1‑150 (MQ=255)
                                                                                                           gTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAActctc                                           <  2:194869/150‑1 (MQ=255)
                                                                                                            tGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAActctct                                          >  2:460878/1‑150 (MQ=255)
                                                                                                               aaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCtttt                                       <  1:1299569/150‑1 (MQ=255)
                                                                                                                aTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCttttt                                      >  1:705908/1‑150 (MQ=255)
                                                                                                                 tAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTg                                     <  2:1501401/150‑1 (MQ=255)
                                                                                                                 tAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTg                                     <  1:903025/150‑1 (MQ=255)
                                                                                                                 tAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTg                                     <  2:593813/150‑1 (MQ=255)
                                                                                                                   aGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTGCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTc                                   <  1:196755/150‑1 (MQ=255)
                                                                                                                    gTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGCCGTGTGTTACATGAGCAACGGAGGTGCTTCACAGCAAAACAATAATCCCCAGTAACTCTCGTGTTGTCa                                  >  1:1305417/1‑150 (MQ=38)
                                                                                                                    gTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCa                                  <  2:1634182/150‑1 (MQ=255)
                                                                                                                    gTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCa                                  >  2:1764384/1‑150 (MQ=255)
                                                                                                                      aTTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGAATTGTGGTAGATGAGCAACGCAGGTGCTAGACAGCAAATCAATCCTCAAAATTAACTCTCTTTTTGTCAAg                                >  2:1815474/1‑150 (MQ=37)
                                                                                                                       ttttACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGc                               <  2:21944/150‑1 (MQ=255)
                                                                                                                         agactagCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACACCCCAGTAACTCTCTTTATGTCAAGCaa                             >  2:326545/7‑150 (MQ=37)
                                                                                                                          tACAGGCTATAAAAAACCGCCGAATTTGGAGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCaaa                            <  1:1181074/150‑1 (MQ=255)
                                                                                                                           aCAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCaaaa                           <  1:579817/150‑1 (MQ=255)
                                                                                                                            cAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAag                          <  2:862902/150‑1 (MQ=255)
                                                                                                                              ggCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAagag                        >  1:1837679/1‑150 (MQ=255)
                                                                                                                                 aaTAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTGCTTGACAAAAAGAGAGt                     <  2:1476875/150‑1 (MQ=38)
                                                                                                                                  aTAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTGCTTGACAAAAAGAGAGtt                    <  2:1502560/150‑2 (MQ=38)
                                                                                                                                  aTAAAAAACCGCAGAATTTGGCGGTTTTTTATTGCAAGTCTGGTTCGCGGCCGTTCCAGCAGGTTTACTTGTGTTACATGAGAAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTGCTTGACAAAAAGAGAGtt                    <  2:1219168/150‑2 (MQ=255)
                                                                                                                                    aaaaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAat                  >  1:397987/1‑150 (MQ=255)
                                                                                                                                    aaaaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAat                  >  2:1627190/1‑150 (MQ=255)
                                                                                                                                     aaaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAatt                 >  2:674617/1‑150 (MQ=255)
                                                                                                                                     aaaaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCATAAGTGAGTAatt                 >  2:1488158/1‑150 (MQ=255)
                                                                                                                                     aaaaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGAATTGAGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCAAAAGTAACTCTATTTTTGTCAAGCAAAAGAGATTAatt                 >  2:1269073/1‑150 (MQ=37)
                                                                                                                                       aaccccccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTGGACTTGTGTTACATGCGCAACGCAGGTGATTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAattat               <  2:161201/144‑1 (MQ=38)
                                                                                                                                       aaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAattat               >  2:107420/1‑150 (MQ=255)
                                                                                                                                       aaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAattat               >  1:1202983/1‑150 (MQ=255)
                                                                                                                                          ccgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGAGTTACATTAGCAAGCGAGTAGCTGCACCGCAAGACAATACTCACAATTAACTCTCTTATTGTCAAGCTATAGAGAGTAAGTATTGt            >  1:1692508/1‑150 (MQ=25)
                                                                                                                                          ccgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTAACCATTAGTTACATGAGAAACGCAGGTGCTTCACAGCAAATCAATACTAACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATGAgtgt            >  1:1517311/1‑150 (MQ=37)
                                                                                                                                           cgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGtt           >  1:889627/1‑150 (MQ=255)
                                                                                                                                           cgccgAAGTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCTCAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGtt           >  1:35110/1‑150 (MQ=255)
                                                                                                                                               gAATTTGGCGGTTTTTTATTGCTAGTCTGTTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGtttatt       <  2:1616308/150‑1 (MQ=255)
                                                                                                                                               gAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCTGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGtttatt       <  2:1513064/150‑1 (MQ=255)
                                                                                                                                                aaTATGTCGGTTTTTTATTGCTAGTCTGGCTCGCGGCCTTTCGAGCATGTTGACTAGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCTCCAGTAACTCTCTTTATGGCAAGCAAACGAGAGTAATTATTGtttattt      <  2:1345736/150‑1 (MQ=255)
                                                                                                                                                    tGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCg  <  1:1474958/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCG  >  NC_000913/2653224‑2653521

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: