Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,009,731 A→G 7.5% intergenic (+19/‑90) yedF → / → yedK putative TusA family sulfurtransferase/DUF159 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,009,7310AG 7.5% 525.7 / 19.7 167intergenic (+19/‑90)yedF/yedKputative TusA family sulfurtransferase/DUF159 family protein
Reads supporting (aligned to +/- strand):  ref base A (80/74);  new base G (6/7);  total (86/81)
Fisher's exact test for biased strand distribution p-value = 7.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 168 total aligned reads displayed.

AATCCTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTCGC  >  NC_000913/2009583‑2009876
                                                                                                                                                    |                                                                                                                                                   
 aTCCTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGAc                                                                                                                                                    >  4:1496043/1‑149 (MQ=255)
  tCCTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGAcc                                                                                                                                                   >  4:634093/1‑149 (MQ=255)
    cTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCca                                                                                                                                                 >  4:1649572/1‑149 (MQ=255)
    cTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCca                                                                                                                                                 >  1:1251516/1‑149 (MQ=255)
     tGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCcac                                                                                                                                                <  4:1539013/149‑1 (MQ=255)
      ggAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCcaca                                                                                                                                               >  4:282816/1‑149 (MQ=255)
          gtggtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGACGTAATCCCTCGTTCCCGCTGTGGGTCACAGCgg                                                                                                                                           <  4:40269/149‑1 (MQ=255)
          gtggtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCgg                                                                                                                                           <  2:1522990/149‑1 (MQ=255)
          gtggtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCgg                                                                                                                                           <  3:1649572/149‑1 (MQ=255)
             gtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGaa                                                                                                                                       <  4:1319617/149‑1 (MQ=255)
                aGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGaatga                                                                                                                                    <  4:606261/149‑4 (MQ=255)
                aGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACAt                                                                                                                                    <  2:1070944/149‑1 (MQ=255)
                 gCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATc                                                                                                                                   <  3:634093/149‑1 (MQ=255)
                 gCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGGAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATc                                                                                                                                   >  4:128900/1‑149 (MQ=255)
                  cGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCg                                                                                                                                  >  1:1414759/1‑149 (MQ=255)
                   gACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGt                                                                                                                                 <  3:1496043/149‑1 (MQ=255)
                          cGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCttt                                                                                                                          <  3:926843/149‑1 (MQ=255)
                          cGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCttt                                                                                                                          <  2:1150878/149‑1 (MQ=255)
                          cGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCttt                                                                                                                          <  4:1087685/149‑1 (MQ=255)
                           gCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTg                                                                                                                         >  3:1039959/1‑149 (MQ=255)
                            cAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGt                                                                                                                        <  2:775585/149‑1 (MQ=255)
                               tCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATc                                                                                                                     <  4:783325/149‑1 (MQ=255)
                                 gATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcac                                                                                                                   >  1:1143970/1‑149 (MQ=255)
                                 gATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcac                                                                                                                   >  2:138727/1‑149 (MQ=255)
                                     aataatATTCCACTGGATGCGCGTAATCACGGCTATACGTTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTAATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAg                                                                                                               >  3:902359/1‑149 (MQ=255)
                                     aataatATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcaca                                                                                                                  >  4:364022/1‑146 (MQ=255)
                                     aataatATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcaca                                                                                                                  <  3:364022/146‑1 (MQ=255)
                                     aataatATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAg                                                                                                               >  3:1612770/1‑149 (MQ=255)
                                     aataatATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAg                                                                                                               >  3:837772/1‑149 (MQ=255)
                                     aataatATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAg                                                                                                               >  3:1408036/1‑149 (MQ=255)
                                      ataataTTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGa                                                                                                              <  2:1087711/149‑1 (MQ=255)
                                          tatTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGAtttt                                                                                                          >  4:1123621/1‑149 (MQ=255)
                                           atTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTa                                                                                                         <  1:652478/149‑1 (MQ=255)
                                           atTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTa                                                                                                         <  4:212974/149‑1 (MQ=255)
                                            ttCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTAc                                                                                                        >  3:409636/1‑149 (MQ=255)
                                            gtCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTCATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTAc                                                                                                        >  4:1093275/2‑149 (MQ=255)
                                               cACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATc                                                                                                                                   >  2:1372924/1‑119 (MQ=255)
                                               cACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCg                                                                                                    >  3:566822/1‑149 (MQ=255)
                                               cACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCg                                                                                                    >  1:420053/1‑149 (MQ=255)
                                                         gcgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTcc                                                                                          <  1:872960/149‑1 (MQ=255)
                                                           gTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTa                                                                                        >  1:400246/1‑149 (MQ=255)
                                                           gTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTa                                                                                        >  1:979773/1‑149 (MQ=255)
                                                             aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGAtttt                                                                                                          <  3:123967/130‑1 (MQ=255)
                                                               tCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCt                                                                                    >  1:46775/1‑149 (MQ=255)
                                                                 aCGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCttt                                                                                  >  2:377955/1‑149 (MQ=255)
                                                                 aCGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCttt                                                                                  >  1:1302625/1‑149 (MQ=255)
                                                                 aCGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCttt                                                                                  >  1:656954/1‑149 (MQ=255)
                                                                    gCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCa                                                                               <  3:1584886/149‑1 (MQ=255)
                                                                         aCGGTGCTTGATTTTCAACAAAGCGGGCCGAACAATCGTTAATTTAATCCGAAGTAATCCCTTATTCCCGCTGTGGCCCCCAGCGGGAAACACCTTTACTTTGTATCACCTTGATTTTACCCGGCCCCCACTCCTATTCTTTTCAGTAt                                                                           <  2:1143970/149‑1 (MQ=255)
                                                                               cTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCa                                                                    <  2:1495468/149‑1 (MQ=255)
                                                                                tGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCag                                                                   >  1:133893/1‑149 (MQ=255)
                                                                                 ggATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCagg                                                                  <  1:807359/149‑1 (MQ=255)
                                                                                 ggATATTCAGCAAGACGCGCCGACCATTCGTTATTTAATTCAGAAGCAATCCCTCATTCCCGCTGTGACCCACAGCGGGACCGTCGTTACCTTTGTATCACATAGATTTTACCCGCCCACACTCCTATTCTTTTCAGTATGATCTCagg                                                                  <  4:1602142/149‑1 (MQ=255)
                                                                                  gATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTAt                                                                           <  3:248684/139‑1 (MQ=255)
                                                                                  gATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCagga                                                                 >  2:1643862/1‑149 (MQ=255)
                                                                                     atTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGt                                                              <  1:1351312/149‑1 (MQ=255)
                                                                                      ttCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGtt                                                             <  2:426237/149‑1 (MQ=255)
                                                                                       tCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTa                                                                                                         >  2:237308/1‑105 (MQ=255)
                                                                                       tCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGttt                                                            <  4:254853/149‑1 (MQ=255)
                                                                                         aGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCt                                                          >  4:1344978/1‑149 (MQ=255)
                                                                                          gCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGTGGTTTCTa                                                         >  1:621981/1‑149 (MQ=255)
                                                                                             aGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtgt                                                      <  4:409636/149‑1 (MQ=255)
                                                                                              gACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCagga                                                                 <  1:361201/137‑1 (MQ=255)
                                                                                              gACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtgtg                                                     >  3:895294/1‑149 (MQ=255)
                                                                                               aCGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCGGTATGATCTCAGGAGGTTTCTAtgtgt                                                    <  3:1088720/149‑1 (MQ=255)
                                                                                               aCGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtgtgt                                                    >  4:1122484/1‑149 (MQ=255)
                                                                                               aCGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtgtgt                                                    <  2:905188/149‑1 (MQ=255)
                                                                                                       aCCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCtt                                            <  3:1325969/149‑1 (MQ=255)
                                                                                                        ccATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGAATGATCTCAGGAGGTTTCTATGTGTGGACGCttt                                           <  3:655212/149‑1 (MQ=255)
                                                                                                          aTTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGc                                         >  1:1035489/1‑149 (MQ=255)
                                                                                                          aTTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGc                                         <  4:858296/149‑1 (MQ=255)
                                                                                                           ttCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtgt                                                      >  3:240794/1‑135 (MQ=255)
                                                                                                           ttCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtgt                                                      <  4:240794/135‑1 (MQ=255)
                                                                                                           ttCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGcc                                        <  2:567991/149‑1 (MQ=255)
                                                                                                            tCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGccc                                       >  1:176897/1‑149 (MQ=255)
                                                                                                                  tttAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCtt                                                                                   >  3:503414/1‑99 (MQ=255)
                                                                                                                  tttAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCtt                                                                                   <  4:503414/99‑1 (MQ=255)
                                                                                                                  tttAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATccc                                 <  2:654735/149‑1 (MQ=255)
                                                                                                                     aaTTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCaaa                              <  1:1133136/149‑1 (MQ=255)
                                                                                                                          aGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCTAATCCCAAACGCGt                         <  3:1413004/149‑1 (MQ=255)
                                                                                                                          aGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGt                         >  1:555718/1‑149 (MQ=255)
                                                                                                                          aGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGt                         >  1:909742/1‑149 (MQ=255)
                                                                                                                             aGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGaa                      <  1:138727/149‑1 (MQ=255)
                                                                                                                               tAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGt                         <  2:422167/144‑1 (MQ=255)
                                                                                                                                 aTCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGAtt                  >  3:833713/1‑149 (MQ=255)
                                                                                                                                  tCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCaa                                     >  2:646656/1‑129 (MQ=255)
                                                                                                                                  tCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTa                 >  4:67468/1‑149 (MQ=255)
                                                                                                                                  tCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTa                 >  1:1169199/1‑149 (MQ=255)
                                                                                                                                    ccTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGaa                      <  3:640504/142‑1 (MQ=255)
                                                                                                                                     cTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCt              >  4:1226439/1‑149 (MQ=255)
                                                                                                                                      tCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCtt             <  4:759728/149‑1 (MQ=255)
                                                                                                                                       cATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTg            <  4:645504/149‑1 (MQ=255)
                                                                                                                                        gttCCCGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTgc           >  2:1518522/2‑149 (MQ=255)
                                                                                                                                        aTTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTgc           <  2:941375/149‑1 (MQ=255)
                                                                                                                                         ttCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGt                         >  1:1077345/1‑134 (MQ=255)
                                                                                                                                         ttCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGt                         <  2:1077345/134‑1 (MQ=255)
                                                                                                                                           cccGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCt        >  2:344669/1‑149 (MQ=255)
                                                                                                                                            ccGCTGTGGGTCACAGCGGGAATATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCtt       >  1:646737/1‑149 (MQ=255)
                                                                                                                                            ccGCTGTGGGTCACAGCGGGAACATCGCTACCTTTGTATCACATAGCTTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGc                                              <  2:1229763/110‑1 (MQ=255)
                                                                                                                                            ccGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCtt       >  3:1139976/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AATCCTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTCGC  >  NC_000913/2009583‑2009876

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