Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,654,144 T→A 5.4% noncoding (290/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,1440TA 5.4% 326.9 / 6.9 92noncoding (290/304 nt)ryfAnovel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base T (43/44);  new base A (2/3);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.88e-01

GCACAATGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCA  >  NC_000913/2653996‑2654292
                                                                                                                                                    |                                                                                                                                                    
gCACAATGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGTTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCttt                                                                                                                                                      <  3:602569/149‑1 (MQ=255)
gCACAATGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCttt                                                                                                                                                      >  2:613116/1‑149 (MQ=255)
  acaATGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTgg                                                                                                                                                    >  4:142368/1‑149 (MQ=255)
    aaTGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggc                                                                                                                                                  <  4:333993/149‑1 (MQ=255)
     aTGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggcg                                                                                                                                                 <  4:447414/149‑1 (MQ=255)
      tGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggcgg                                                                                                                                                <  3:381737/149‑1 (MQ=255)
       gAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggcggc                                                                                                                                               >  3:792314/1‑149 (MQ=255)
       gAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTCggcggc                                                                                                                                               <  4:616275/149‑1 (MQ=255)
         aGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCat                                                                                                                                             <  4:844438/149‑1 (MQ=255)
           aTACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCatat                                                                                                                                           >  4:491970/1‑149 (MQ=255)
           aTACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCattc                                                                                                                                           <  2:666844/149‑3 (MQ=255)
              cGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCattcggg                                                                                                                                        <  4:194037/149‑6 (MQ=255)
               ggTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTa                                                                                                                                       <  3:854767/149‑1 (MQ=255)
               ggTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTa                                                                                                                                       <  2:380997/149‑1 (MQ=255)
               ggTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTa                                                                                                                                       >  3:43834/1‑149 (MQ=255)
                 tGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGa                                                                                                                                     <  4:872588/149‑1 (MQ=255)
                    ttttGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTa                                                                                                                                  >  3:52059/1‑149 (MQ=255)
                    ttttGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTa                                                                                                                                  >  3:156105/1‑149 (MQ=255)
                    ttttGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTa                                                                                                                                  >  2:846525/1‑149 (MQ=255)
                    ttttGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTa                                                                                                                                  >  2:855021/1‑149 (MQ=255)
                    ttttGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTa                                                                                                                                  >  4:156481/1‑149 (MQ=255)
                    ttttGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTa                                                                                                                                  >  1:379667/1‑149 (MQ=255)
                      ttGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATc                                                                                                                                <  3:142368/149‑1 (MQ=255)
                      ttGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATc                                                                                                                                <  4:179560/149‑1 (MQ=255)
                      ttGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATc                                                                                                                                <  4:381870/149‑1 (MQ=255)
                           tCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGAtt                                                                                                                           >  2:202516/1‑149 (MQ=255)
                                 ttATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGttt                                                                                                                     <  4:488778/149‑1 (MQ=255)
                                 ttATTGTTTCTGGTGCGCTGTTAACCCAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATCTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCAATGGCCACCTTTGGGCGGCATCATTTAGATTATTCGGCTCCGTTTa                                                                                                                    >  1:447132/1‑149 (MQ=255)
                                   aTTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTGTAGATTATCCGATTCTGTTTaa                                                                                                                   >  2:746406/1‑149 (MQ=255)
                                    ttGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTaaa                                                                                                                  <  2:726059/149‑1 (MQ=255)
                                    ttGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTaaa                                                                                                                  <  4:52059/149‑1 (MQ=255)
                                    ttGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTaaa                                                                                                                  <  1:226391/149‑1 (MQ=255)
                                         tCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACCAGTTACTGGTGGTTAACAGTCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTTGTTTATCCGATTCTGTTTAAAGTCAc                                                                                                             >  4:849855/1‑149 (MQ=255)
                                         tCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCAc                                                                                                             <  3:27535/149‑1 (MQ=255)
                                               gcgcTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCaaaa                                                                                                       <  4:521556/149‑1 (MQ=255)
                                                 gcTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCaaaaaa                                                                                                     >  3:453927/1‑149 (MQ=255)
                                                  cTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAac                                                                                                    >  1:697263/1‑149 (MQ=255)
                                                  cTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAac                                                                                                    <  3:775139/149‑1 (MQ=255)
                                                   tGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAacc                                                                                                   >  2:506917/1‑149 (MQ=255)
                                                    gTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAacca                                                                                                  <  3:710078/149‑1 (MQ=255)
                                                     ttAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAaccac                                                                                                 >  2:369804/1‑149 (MQ=255)
                                                     ttAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAaccac                                                                                                 <  4:43834/149‑1 (MQ=255)
                                                              gTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGAc                                                                                        >  4:431019/1‑149 (MQ=255)
                                                               tAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACg                                                                                       >  1:116318/1‑149 (MQ=255)
                                                               tAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACg                                                                                       >  3:504409/1‑149 (MQ=255)
                                                                  ataataACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTc                                                                                    <  3:877683/149‑1 (MQ=255)
                                                                   taataaCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTg                                                                                                                                                     <  2:62222/83‑1 (MQ=255)
                                                                   taataaCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTg                                                                                                                                                     >  1:62222/1‑83 (MQ=255)
                                                                    actaaCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCAt                                                                                  <  2:177624/147‑1 (MQ=255)
                                                                     ataaCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATa                                                                                 <  1:846525/149‑1 (MQ=255)
                                                                     ataaCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATa                                                                                 <  1:420744/149‑1 (MQ=255)
                                                                       aaCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGa                                                                               >  1:597598/1‑149 (MQ=255)
                                                                        aCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGaa                                                                              >  1:334043/1‑149 (MQ=255)
                                                                        aCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGaa                                                                              >  1:493518/1‑149 (MQ=255)
                                                                          cGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATg                                                                            <  2:524223/149‑1 (MQ=255)
                                                                            gAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGgc                                                                          >  1:138680/1‑149 (MQ=255)
                                                                                 tctcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAat                                                                     >  2:214914/1‑149 (MQ=255)
                                                                                 tctcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAat                                                                     >  4:859755/1‑149 (MQ=255)
                                                                                   tcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAatat                                                                   >  1:147349/1‑149 (MQ=255)
                                                                                       gCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTc                                                               <  3:591532/149‑1 (MQ=255)
                                                                                       gCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTc                                                               <  1:379856/149‑1 (MQ=255)
                                                                                       gCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTc                                                               <  1:202516/149‑1 (MQ=255)
                                                                                        cGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCg                                                              <  1:214914/149‑1 (MQ=255)
                                                                                           cgaTTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCCTTGGGCGGCATATTTTAGATTCTCCGATTCTGTTCAAAGTCTCGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGta                                                           <  3:849855/147‑1 (MQ=255)
                                                                                              tttACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGtaata                                                        >  2:292250/1‑149 (MQ=255)
                                                                                              tttACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGtaata                                                        <  1:855021/149‑1 (MQ=255)
                                                                                               ttACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGtaataa                                                       >  2:841344/1‑149 (MQ=255)
                                                                                                 aCTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAg                                                     <  1:613116/149‑1 (MQ=255)
                                                                                                     gtggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCg                                                              <  4:721730/136‑1 (MQ=255)
                                                                                                     gtggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCg                                                              >  3:721730/1‑136 (MQ=255)
                                                                                                        gtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTaa                                              <  2:697263/149‑1 (MQ=255)
                                                                                                           aacaacCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTc                                           >  3:136881/1‑149 (MQ=255)
                                                                                                              aacCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCtt                                        >  4:194119/1‑149 (MQ=255)
                                                                                                                 cTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAac                                                                                                    >  2:724981/1‑86 (MQ=255)
                                                                                                                 cTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAac                                                                                                    <  1:724981/86‑1 (MQ=255)
                                                                                                                  ttCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTc                                    >  3:434710/1‑149 (MQ=255)
                                                                                                                  ttCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTc                                    >  2:625428/1‑149 (MQ=255)
                                                                                                                  ttCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTc                                    <  2:597598/149‑1 (MQ=255)
                                                                                                                  ttCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTc                                    >  2:182369/1‑149 (MQ=255)
                                                                                                                  ttCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTCTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTc                                    >  2:796582/1‑149 (MQ=255)
                                                                                                                       agcagcAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGGCGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCtt                               >  3:641869/1‑149 (MQ=255)
                                                                                                                         cagcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAAATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCtttt                             <  1:475382/149‑1 (MQ=255)
                                                                                                                            caAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTc                          <  3:491970/149‑1 (MQ=255)
                                                                                                                             aaGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCa                         <  2:493518/149‑1 (MQ=255)
                                                                                                                              aGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCAc                        <  1:30116/149‑1 (MQ=255)
                                                                                                                                   taaaataTGCCGCCCAAAGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGAc                   >  3:863829/7‑149 (MQ=255)
                                                                                                                                    aaaataTGCCGCCCAAAGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACa                  >  3:845699/6‑149 (MQ=255)
                                                                                                                                      cGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACaaa                >  4:514603/1‑149 (MQ=255)
                                                                                                                                         aTGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAatat             <  3:156481/149‑1 (MQ=255)
                                                                                                                                            ccgccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGa          <  3:514603/149‑1 (MQ=255)
                                                                                                                                             cgccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGAt         <  1:20326/149‑1 (MQ=255)
                                                                                                                                             cgccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGAt         <  1:369804/149‑1 (MQ=255)
                                                                                                                                              gccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGAtt        >  1:524105/1‑149 (MQ=255)
                                                                                                                                               cccTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTg       >  2:534684/1‑149 (MQ=255)
                                                                                                                                                  tttGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCGATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGtt    >  4:144319/1‑149 (MQ=255)
                                                                                                                                                    tGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCa  <  1:116798/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCACAATGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTTCA  >  NC_000913/2653996‑2654292

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