Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,654,142 T→A 5.0% noncoding (288/304 nt) ryfA → novel sRNA, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,654,1420TA 5.0% 335.6 / 6.6 101noncoding (288/304 nt)ryfAnovel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base T (44/51);  new base A (2/3);  total (46/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

Only 100 of 102 total aligned reads displayed.

GAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTT  >  NC_000913/2654003‑2654290
                                                                                                                                           |                                                                                                                                                     
gAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggcgg                                                                                                                                               >  2:863888/1‑148 (MQ=255)
gAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggcgg                                                                                                                                               <  2:413016/148‑1 (MQ=255)
   gATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCat                                                                                                                                            <  1:443511/148‑1 (MQ=255)
    aTACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCata                                                                                                                                           >  1:954770/1‑148 (MQ=255)
    aTACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAATGGGCGGCatt                                                                                                                                           <  1:717155/148‑2 (MQ=255)
    aTACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCattc                                                                                                                                          <  4:835790/149‑3 (MQ=255)
       cGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAGGCCCGAATGCCGCCCTTTGGGCGGCATAttt                                                                                                                                        >  2:343294/1‑148 (MQ=255)
        ggTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATAtttt                                                                                                                                       <  1:618013/148‑1 (MQ=255)
        ggTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATAtttt                                                                                                                                       <  1:521151/148‑1 (MQ=255)
        ggTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATAtttt                                                                                                                                       <  1:95246/148‑1 (MQ=255)
        ggTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATAtttt                                                                                                                                       <  1:217009/148‑1 (MQ=255)
        ggTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTa                                                                                                                                      <  4:827707/149‑1 (MQ=255)
         gTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTa                                                                                                                                      >  1:541652/1‑148 (MQ=255)
             ttttGTATCGTACTTATTGTTTCTGGTGCGCTGTTATCCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGGATGCCGCCCTTTGGGCGGCATATTTTAGATTa                                                                                                                                 >  3:803567/1‑149 (MQ=255)
             ttttGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTa                                                                                                                                 >  4:821692/1‑149 (MQ=255)
              tttGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTa                                                                                                                                 >  2:479471/1‑148 (MQ=255)
               ttGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTAt                                                                                                                                <  1:870539/148‑1 (MQ=255)
                        aCTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTg                                                                                                                       <  2:185723/148‑1 (MQ=255)
                        aCTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGt                                                                                                                      <  3:821692/149‑1 (MQ=255)
                         cTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGt                                                                                                                      <  1:264627/148‑1 (MQ=255)
                          ttATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTCGGGCGGCATATTTTAGATTATCCGATTCTGttt                                                                                                                    >  4:840111/1‑149 (MQ=255)
                           tATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGttt                                                                                                                    >  1:839457/1‑148 (MQ=255)
                             ttGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAGGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTaa                                                                                                                  <  1:343294/148‑1 (MQ=255)
                             ttGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTaaa                                                                                                                 <  3:225541/149‑1 (MQ=255)
                                tttCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGt                                                                                                               <  1:83099/148‑1 (MQ=255)
                                  tCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCa                                                                                                             <  1:863888/148‑1 (MQ=255)
                                   cTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACg                                                                                                           <  3:301536/149‑1 (MQ=255)
                                     ggTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGc                                                                                                          >  1:19536/1‑148 (MQ=255)
                                           ctcttAACCGAGGTAAATAATAACCGGAGTCTTTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTGTAGATTATCCGATTCTGTTTAAAGTCAGGCAAAATac                                                                                                  >  4:609397/4‑149 (MQ=255)
                                           cTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCaaaaaa                                                                                                   <  1:586732/148‑1 (MQ=255)
                                            tGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCAc                                                                                                            >  3:408318/1‑139 (MQ=255)
                                            tGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCAc                                                                                                            <  4:408318/139‑1 (MQ=255)
                                             gTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAacca                                                                                                <  3:401202/149‑1 (MQ=255)
                                             gTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAacc                                                                                                 >  2:501725/1‑148 (MQ=255)
                                                 aCCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCAcccc                                                                                            >  3:521821/1‑149 (MQ=255)
                                                        tAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATc                                                                                                                               <  4:499205/108‑1 (MQ=255)
                                                        tAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATc                                                                                                                               >  3:499205/1‑108 (MQ=255)
                                                          aataatAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGt                                                                                    <  2:541652/148‑1 (MQ=255)
                                                           ataataACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTc                                                                                  <  3:817194/149‑1 (MQ=255)
                                                           ataataACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTc                                                                                  >  3:475694/1‑149 (MQ=255)
                                                            taataaCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCa                                                                                 >  3:275052/1‑149 (MQ=255)
                                                             aataaCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCAt                                                                                >  3:239645/1‑149 (MQ=255)
                                                              ataaCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCAt                                                                                >  2:852470/1‑148 (MQ=255)
                                                                 aCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGaa                                                                            <  4:759107/149‑1 (MQ=255)
                                                                 aCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGa                                                                             >  2:421149/1‑148 (MQ=255)
                                                                 aCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGa                                                                             <  1:777044/148‑1 (MQ=255)
                                                                  ccGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGaa                                                                            <  1:852470/148‑1 (MQ=255)
                                                                   ctgagaCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGTGGCATATTTTAGATTAGCCGATTCTTTTTAACGTCACGAAACAATCCACCCCAGCGACGTTCATAGact                                                                           <  2:391675/146‑3 (MQ=255)
                                                                   cGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAAt                                                                           <  2:408313/148‑1 (MQ=255)
                                                                    ggAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCCGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATg                                                                          <  2:142129/148‑1 (MQ=255)
                                                                     gAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGgc                                                                        >  3:10009/1‑149 (MQ=255)
                                                                          tctcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAat                                                                   >  3:741079/1‑149 (MQ=255)
                                                                           ctcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAat                                                                   >  1:303568/1‑148 (MQ=255)
                                                                                gCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTc                                                             <  4:182687/149‑1 (MQ=255)
                                                                                  gACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCgg                                                           <  3:174666/149‑1 (MQ=255)
                                                                                  gACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCg                                                            <  1:226625/148‑1 (MQ=255)
                                                                                        ttACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGCTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGtaata                                                      >  1:247681/1‑148 (MQ=255)
                                                                                          aCTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAg                                                   <  4:149676/149‑1 (MQ=255)
                                                                                           cTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAg                                                   <  2:303568/148‑1 (MQ=255)
                                                                                            tggtggTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTg                                                 >  4:265690/1‑149 (MQ=255)
                                                                                             ggtggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCGGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGt                                                  <  2:698507/147‑1 (MQ=255)
                                                                                             ggtggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGt                                                  >  2:631725/1‑147 (MQ=255)
                                                                                             ggtggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGt                                                  <  1:631725/147‑1 (MQ=255)
                                                                                             ggtggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGt                                                  >  1:698507/1‑147 (MQ=255)
                                                                                             ggtggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTgg                                                <  4:452366/149‑1 (MQ=255)
                                                                                              gtggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTc                                                             >  4:711923/1‑135 (MQ=255)
                                                                                              gtggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTc                                                             <  3:711923/135‑1 (MQ=255)
                                                                                              gtggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGc                                               >  4:632071/1‑149 (MQ=255)
                                                                                                 gtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTaa                                            >  4:905745/1‑149 (MQ=255)
                                                                                                   tAACGACCTTGAGAGCAGCAAGTAAGCCCGAATGTAACATCTTGGACGGCAGATTTTAGATTATCCGATTCTATTTAAAGTCACGCAAGCAACCACCCCGGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAAtt                                          <  4:502900/149‑1 (MQ=255)
                                                                                                   tAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAAt                                           <  2:839457/148‑1 (MQ=255)
                                                                                                      caacCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTcc                                        <  1:501725/148‑1 (MQ=255)
                                                                                                      caacCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTcc                                        >  1:520346/1‑148 (MQ=255)
                                                                                                       aacCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCt                                       <  2:855507/148‑1 (MQ=255)
                                                                                                           ttCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGt                                                  >  4:357615/1‑133 (MQ=255)
                                                                                                           ttCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGt                                                  <  3:357615/133‑1 (MQ=255)
                                                                                                            tCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTcc                                        <  1:358942/142‑1 (MQ=255)
                                                                                                            tCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTcc                                        >  2:358942/1‑142 (MQ=255)
                                                                                                            tCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTc                                  >  2:885482/1‑148 (MQ=255)
                                                                                                             cAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGtt                                                                                   <  3:534845/98‑1 (MQ=255)
                                                                                                             cAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGtt                                                                                   >  4:534845/1‑98 (MQ=255)
                                                                                                              agagCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTcccc                               >  4:176138/1‑149 (MQ=255)
                                                                                                                agcagcAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCtt                             <  3:265690/149‑1 (MQ=255)
                                                                                                                  cagcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCtttt                           <  4:521821/149‑1 (MQ=255)
                                                                                                                        gTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCAcc                     >  4:615287/1‑149 (MQ=255)
                                                                                                                         tAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTCTCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCCTTTTTCACCt                    <  4:622178/149‑1 (MQ=255)
                                                                                                                         tAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCt                    >  3:496533/1‑149 (MQ=255)
                                                                                                                             cccGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGAc                 >  2:451168/1‑148 (MQ=255)
                                                                                                                                 taTGCCGCCCAAAGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAata            >  4:836451/2‑149 (MQ=255)
                                                                                                                                 taTGCCGCCCAAAGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAata            >  4:861404/2‑149 (MQ=255)
                                                                                                                                      cgccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGa        >  1:414144/1‑148 (MQ=255)
                                                                                                                                      cgccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGAt       <  3:176138/149‑1 (MQ=255)
                                                                                                                                       gccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGa        <  2:653395/148‑1 (MQ=255)
                                                                                                                                       gccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAGAACCACCCCAGCGACGTTCATAGAATGGCGCAATGTGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGAtt      <  3:384110/149‑1 (MQ=255)
                                                                                                                                       gccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGAtt      >  4:474680/1‑149 (MQ=255)
                                                                                                                                       gccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGAt       >  1:498508/1‑148 (MQ=255)
                                                                                                                                        cccTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTg     <  4:514502/149‑1 (MQ=255)
                                                                                                                                        cccTTTGGGCAGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGAtt      <  2:1501/148‑1 (MQ=255)
                                                                                                                                         ccTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTgg    >  4:216041/1‑149 (MQ=255)
                                                                                                                                           tttGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGtt  >  3:305097/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                     
GAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTT  >  NC_000913/2654003‑2654290

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: