Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,009,731 A→G 5.7% 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 5.7% 361.3 / 8.9 122intergenic (+19/‑90)yedF/yedKputative TusA family sulfurtransferase/DUF159 family protein
Reads supporting (aligned to +/- strand):  ref base A (47/67);  new base G (3/4);  total (51/71)
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.01e-01

Only 100 of 125 total aligned reads displayed.

AAATCCTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTCGCGGAA  >  NC_000913/2009582‑2009880
                                                                                                                                                     |                                                                                                                                                     
aaaTGCTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGAGGGGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGa                                                                                                                                                       <  1:78721/150‑1 (MQ=255)
    ccTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCca                                                                                                                                                   >  1:483923/1‑150 (MQ=255)
     cTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCATAATCAGGTCAATAAGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTGATATAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCcac                                                                                                                                                  >  2:142305/1‑150 (MQ=255)
      tGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCcaca                                                                                                                                                 <  2:393393/150‑1 (MQ=255)
         aaGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCg                                                                                                                                              <  1:1351989/150‑1 (MQ=255)
          aGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCgg                                                                                                                                             <  1:1326726/150‑1 (MQ=255)
              gtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGaat                                                                                                                                         <  2:1002077/150‑2 (MQ=255)
               tgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGAACATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACCGCGGGcaca                                                                                                                                        >  1:1652535/1‑150 (MQ=255)
               ggAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTTGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGaatg                                                                                                                                        <  2:2157941/149‑3 (MQ=255)
                  gCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCg                                                                                                                                     <  2:1333673/150‑1 (MQ=255)
                  gCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGGGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCg                                                                                                                                     <  2:1242164/150‑1 (MQ=255)
                     aCTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTa                                                                                                                                  <  1:25963/150‑1 (MQ=255)
                     aCTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTa                                                                                                                                  <  2:1971537/150‑1 (MQ=255)
                      cTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGATATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTAAGAATTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTAc                                                                                                                                 >  2:2135060/1‑150 (MQ=255)
                       tGTCCGCTGTCGATCAATAATATTCCACTGGATGCGCGTAATCAAGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTAcc                                                                                                                                >  2:2115002/1‑150 (MQ=255)
                        gTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCt                                                                                                                               <  1:2188634/150‑1 (MQ=255)
                           cGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTg                                                                                                                            >  2:962450/1‑150 (MQ=255)
                           cGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACAATTCGTTATTTAATTCATAAGTAATCCCTCATTCCTGCTGTGACCCACAGCGGGAACATCGTTACCTTTg                                                                                                                            <  2:1539871/150‑1 (MQ=255)
                               gTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTAta                                                                                                                        >  1:1871275/1‑149 (MQ=255)
                               gTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATc                                                                                                                        <  1:787288/150‑1 (MQ=255)
                               gTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATc                                                                                                                        <  2:632881/150‑1 (MQ=255)
                               gTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATc                                                                                                                        <  1:1995998/150‑1 (MQ=255)
                                tCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATca                                                                                                                       <  1:1638053/150‑1 (MQ=255)
                                tCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATca                                                                                                                       <  1:1723747/150‑1 (MQ=255)
                                tCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATca                                                                                                                       <  1:62405/150‑1 (MQ=255)
                                tCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATca                                                                                                                       <  1:615782/150‑1 (MQ=255)
                                 cGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcac                                                                                                                      <  2:1339504/150‑1 (MQ=255)
                                 cGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcac                                                                                                                      >  2:2079847/1‑150 (MQ=255)
                                 cGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTAATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcac                                                                                                                      <  1:1873886/150‑1 (MQ=255)
                                  gATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcaca                                                                                                                     >  1:1792272/1‑150 (MQ=255)
                                  gATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcaca                                                                                                                     >  2:807112/1‑150 (MQ=255)
                                       ataattTTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGCATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGATGTAATCCCTCATACCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGAt                                                                                                                <  1:1101059/150‑1 (MQ=255)
                                         aataTTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGAttt                                                                                                              >  2:963648/1‑150 (MQ=255)
                                          atatTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGATACCTTTGTATCACATAGAtttt                                                                                                             >  1:1129149/1‑150 (MQ=255)
                                          atatTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCCTAACCTTAGTAACACAAAGAtgtt                                                                                                             >  2:1844911/1‑150 (MQ=255)
                                           tatTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTa                                                                                                            <  2:1241390/150‑1 (MQ=255)
                                            atTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTTGATATTCAGAAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTAc                                                                                                           <  1:211991/150‑1 (MQ=255)
                                            atTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTAc                                                                                                           <  1:1401835/150‑1 (MQ=255)
                                               ccACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCg                                                                                                        >  2:162228/1‑150 (MQ=255)
                                               ccACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCg                                                                                                        >  1:1494256/1‑150 (MQ=255)
                                                  cTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGccc                                                                                                     >  1:1085321/1‑150 (MQ=255)
                                                   tGGATGCGCTTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGcccc                                                                                                    <  2:2016129/150‑1 (MQ=255)
                                                       tGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACt                                                                                                <  1:2001748/150‑1 (MQ=255)
                                                         cgcgTAATCACGGCTATATGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTcc                                                                                              <  1:1100684/150‑1 (MQ=255)
                                                         cgcgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTcc                                                                                              <  1:740873/150‑1 (MQ=255)
                                                         cgcgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTcc                                                                                              <  1:874041/150‑1 (MQ=255)
                                                         cgcgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTcc                                                                                              <  1:1255313/150‑1 (MQ=255)
                                                           cgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTa                                                                                            >  1:1684076/1‑150 (MQ=255)
                                                                  aCGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCtttt                                                                                     >  2:66475/1‑150 (MQ=255)
                                                                   cGGCTATTCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCTTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTc                                                                                    <  1:1611724/150‑1 (MQ=255)
                                                                   cGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTc                                                                                    <  1:78185/150‑1 (MQ=255)
                                                                    ggCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCa                                                                                   <  1:276162/150‑1 (MQ=255)
                                                                       tataCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTATCCGCCCCCACTCCTATTCTTTTCAGTa                                                                                <  1:543344/150‑1 (MQ=255)
                                                                         cgtcttGCTGAATATCCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGATGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACTTAGATTTTACCCGCCCCCACTCATATTCTTTTCAGTATg                                                                              >  1:1808400/6‑150 (MQ=255)
                                                                            ggTGCTGGATATTCAGCAAGACGGTCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCTCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGAtc                                                                           <  1:2115002/150‑1 (MQ=255)
                                                                            ggTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGAtc                                                                           >  1:35984/1‑150 (MQ=255)
                                                                              tGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGAtctc                                                                         <  2:1100432/150‑1 (MQ=255)
                                                                                   tatatTCAGGAAGACGGGCCGAGCATTCTTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCaggag                                                                    <  2:188055/149‑1 (MQ=255)
                                                                                       ttCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGttt                                                                >  2:579304/1‑150 (MQ=255)
                                                                                          aGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTa                                                             <  2:1170720/150‑1 (MQ=255)
                                                                                           gCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACCGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAt                                                            >  1:1795499/1‑150 (MQ=255)
                                                                                           gCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAATAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAt                                                            >  1:1733724/1‑150 (MQ=255)
                                                                                              aGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtgtg                                                         <  2:483923/150‑1 (MQ=255)
                                                                                                  gggCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGa                                                     <  1:1013311/150‑1 (MQ=255)
                                                                                                    gCCGACCATTCGTTTTTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACg                                                   <  2:1675681/150‑1 (MQ=255)
                                                                                                    gCCGACCATTCGTTATTTAAATCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACg                                                   <  1:420041/150‑1 (MQ=255)
                                                                                                        aCCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCttt                                               >  2:4675/1‑150 (MQ=255)
                                                                                                        aCCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCttt                                               <  2:1766915/150‑1 (MQ=255)
                                                                                                          cATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTAGGATCTCAGGAGGTTTATATGTGTGTCCGCTATGc                                             >  2:393613/1‑150 (MQ=255)
                                                                                                           aatCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGcc                                            <  1:1941014/148‑1 (MQ=255)
                                                                                                           aTTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGcc                                            >  1:10291/1‑150 (MQ=255)
                                                                                                            ttCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGccc                                           <  1:208454/150‑1 (MQ=255)
                                                                                                            ttCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGATTTGccc                                           >  1:268797/1‑150 (MQ=255)
                                                                                                                ttATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATc                                       <  1:66475/150‑1 (MQ=255)
                                                                                                                     tAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTATATGTGTGGACGCTTTGCCCAATCCCaaa                                  >  1:1665206/1‑150 (MQ=255)
                                                                                                                      aaTTCAGAAGTAATCCCTCATTCCCGCTGTGCCCAAAAGCGGGAACATCGTTTCCGTTGTATCACATAGATTTTACCCTCCCCCACGCCTAGTGTTTTCAGTATGATCTCATGAGGTTTCTATGTGTTGACGCTTGGCCCAATCCCaaaa                                 >  1:1566023/1‑149 (MQ=255)
                                                                                                                      aaTTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAAc                                 >  1:15369/1‑150 (MQ=255)
                                                                                                                         tCAGAAGTAATCCCTCATTCCAGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAAcgcg                              >  2:624030/1‑150 (MQ=255)
                                                                                                                            gAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGa                           <  2:1129149/150‑1 (MQ=255)
                                                                                                                            gAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGa                           >  1:1959266/1‑150 (MQ=255)
                                                                                                                               gTAATCCCTCATTCCCGCTGTGACCCACAGCGGGTACATCGTTACCTTTGTATCCCATAGATTTTACCCGCCCCCCCTCCTAATCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGATTTGCCCGATCCCAAACGCGGGAAGa                        >  1:432192/1‑150 (MQ=255)
                                                                                                                               gTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGa                        <  1:655228/150‑1 (MQ=255)
                                                                                                                                 aaTCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTGTATGTGTGGCCGCTTTGCCCAATCCCAAACGCGTGAAGAtt                      >  1:1836718/1‑150 (MQ=255)
                                                                                                                                 aaTCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGAtt                      >  1:942997/1‑150 (MQ=255)
                                                                                                                                  aTCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTa                     >  2:633345/1‑150 (MQ=255)
                                                                                                                                   tCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTAc                    <  2:540445/150‑1 (MQ=255)
                                                                                                                                   tCCCTCATTACCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTAc                    <  2:1648309/150‑1 (MQ=255)
                                                                                                                                     ccTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCt                  <  2:191835/150‑1 (MQ=255)
                                                                                                                                        cATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCATTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTATGCCCAATCCCAAACGCGTGAAGATTACCTTgc               >  2:1208327/1‑150 (MQ=255)
                                                                                                                                           tCCCGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTGCCCGCCCCCCCTCCTATTCTTTTCAGTTTGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCt            >  1:1037597/1‑150 (MQ=255)
                                                                                                                                           tCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCt            >  2:1406864/1‑150 (MQ=255)
                                                                                                                                            cccGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCtt           >  2:618768/1‑150 (MQ=255)
                                                                                                                                            cccGCTGTGGGTCACAGCGGGAACATCGTTACCTATCTATCACATAGATTTTACCCGCCCCCACGCCTAGTCGTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCCATCCCACACGCGTGAAGATAACCTTGCGCtt           >  1:1367448/1‑150 (MQ=255)
                                                                                                                                                 tgtgACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTcgc      <  2:604483/150‑1 (MQ=255)
                                                                                                                                                  gtgACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTcgcg     <  2:1494256/150‑1 (MQ=255)
                                                                                                                                                  gtgACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTcgcg     <  1:914575/150‑1 (MQ=255)
                                                                                                                                                  gtgACCCACAGCGGGAACATCGTTACATTTGTATCACATAGATTTTTCCCGCCCCCACTCCTATTCTTTTCATTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTcgcg     <  2:890412/150‑1 (MQ=255)
                                                                                                                                                    gACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCTGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTGGCGGa   >  2:1244744/1‑150 (MQ=255)
                                                                                                                                                    gACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTCGCGGa   >  1:939603/1‑150 (MQ=255)
                                                                                                                                                     aCCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTATCGCGGaa  >  2:1560018/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAATCCTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTCGCGGAA  >  NC_000913/2009582‑2009880

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: