Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,009,731 A→G 5.4% intergenic (+19/‑90) yedF → / → yedK putative sulfurtransferase YedF/putative SOS response‑associated peptidase YedK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,009,7310AG 5.4% 283.5 / 5.0 96intergenic (+19/‑90)yedF/yedKputative sulfurtransferase YedF/putative SOS response‑associated peptidase YedK
Reads supporting (aligned to +/- strand):  ref base A (44/45);  new base G (2/3);  total (47/49)
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.64e-01

AATCCTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTC  >  NC_000913/2009583‑2009874
                                                                                                                                                    |                                                                                                                                               
aaTCCTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGAc                                                                                                                                                <  1:1300505/150‑1 (MQ=255)
      ggAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAg                                                                                                                                          >  1:710904/1‑150 (MQ=255)
        aaGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCCGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCg                                                                                                                                        <  1:309642/150‑1 (MQ=255)
         aGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCgg                                                                                                                                       <  1:1178252/150‑1 (MQ=255)
         aGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCgg                                                                                                                                       <  2:1801916/150‑1 (MQ=255)
          gtggtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCggg                                                                                                                                      <  2:558841/150‑1 (MQ=255)
          gtggtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCggg                                                                                                                                      <  2:1045703/150‑1 (MQ=255)
          gtggtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCggg                                                                                                                                      <  1:1904341/150‑1 (MQ=255)
          gtggtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGa                                                                                                                                     >  1:1696401/1‑150 (MQ=255)
           tggtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGa                                                                                                                                     <  2:840694/150‑1 (MQ=255)
            ggtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGaa                                                                                                                                    <  1:1301145/150‑1 (MQ=255)
            ggtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGaa                                                                                                                                    >  2:511400/1‑150 (MQ=255)
             gtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAAc                                                                                                                                   <  1:978303/150‑1 (MQ=255)
                   gACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGtt                                                                                                                             <  1:511400/150‑1 (MQ=255)
                      tGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTAcc                                                                                                                          <  1:196789/150‑1 (MQ=255)
                       gTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCt                                                                                                                         >  1:1486055/1‑150 (MQ=255)
                       gTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCt                                                                                                                         <  2:842451/150‑1 (MQ=255)
                         ccGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCttt                                                                                                                       >  2:258399/1‑150 (MQ=255)
                         ccGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCttt                                                                                                                       >  2:1133644/1‑150 (MQ=255)
                          cGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTg                                                                                                                      >  1:719814/1‑150 (MQ=255)
                              gTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATc                                                                                                                  <  1:919213/150‑1 (MQ=255)
                               tCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATca                                                                                                                 <  2:1093421/150‑1 (MQ=255)
                               tCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATca                                                                                                                 <  2:1467919/150‑1 (MQ=255)
                                cGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCACCAAGACGGGCCGACCATTCGTTATTTAATTAAGAAGAAAACCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTAGTAtaaa                                                                                                                >  1:1852487/1‑147 (MQ=255)
                                    taataatATTCCACTGGATGCGCTTAATCACGTCTATTCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAg                                                                                                            <  1:1970832/149‑1 (MQ=255)
                                     aataatATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGa                                                                                                           >  2:654825/1‑150 (MQ=255)
                                     aataatATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGa                                                                                                           <  2:1153748/150‑1 (MQ=255)
                                      ataataTTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGAt                                                                                                          >  1:663541/1‑150 (MQ=255)
                                              ccACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCg                                                                                                  >  1:1833294/1‑150 (MQ=255)
                                                aCTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGcc                                                                                                <  1:258399/150‑1 (MQ=255)
                                                   ggATGCGCGTAATCCCGGCTATACGGTGCTGGATATTCAGCAAGCCGGGCCGCCCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGCCCCACAGCGGGAACATCGTTACCTATGTATCACATAGATGTTACCCGccccc                                                                                             >  2:1836624/1‑150 (MQ=255)
                                                    gATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCa                                                                                            <  1:557654/150‑1 (MQ=255)
                                                      tGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACt                                                                                          <  1:1492860/150‑1 (MQ=255)
                                                       gcgcGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTc                                                                                         <  2:503446/150‑1 (MQ=255)
                                                          cgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTa                                                                                      >  1:480972/1‑150 (MQ=255)
                                                          cgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTa                                                                                      >  1:488061/1‑150 (MQ=255)
                                                          cgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTa                                                                                      >  2:1492251/1‑150 (MQ=255)
                                                              aTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCt                                                                                  >  1:742894/1‑150 (MQ=255)
                                                              aTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCt                                                                                  >  2:524157/1‑150 (MQ=255)
                                                              aTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCt                                                                                  >  1:369553/1‑150 (MQ=255)
                                                               tCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCtt                                                                                 <  1:244101/150‑1 (MQ=255)
                                                                cACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCttt                                                                                >  1:1717343/1‑150 (MQ=255)
                                                                  cGGCTATACGGTGCTGGATCTTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCCGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTc                                                                              <  1:685102/150‑1 (MQ=255)
                                                                  cGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTc                                                                              <  2:1717343/150‑1 (MQ=255)
                                                                  cGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTc                                                                              >  2:597114/1‑150 (MQ=255)
                                                                  cGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTc                                                                              <  2:459580/150‑1 (MQ=255)
                                                                   ggCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCa                                                                             <  2:1215805/150‑1 (MQ=255)
                                                                    gCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAg                                                                            >  1:88868/1‑150 (MQ=255)
                                                                     cTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGt                                                                           <  2:1453979/150‑1 (MQ=255)
                                                                       ataCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTAt                                                                         >  2:667781/1‑150 (MQ=255)
                                                                       ataCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTAt                                                                         >  1:637049/1‑150 (MQ=255)
                                                                         aCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGa                                                                       <  1:597114/150‑1 (MQ=255)
                                                                           ggTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGAtc                                                                     <  1:1492251/150‑1 (MQ=255)
                                                                                 ggATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCagga                                                               >  2:344908/1‑150 (MQ=255)
                                                                                  gATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCaggag                                                              <  2:369553/150‑1 (MQ=255)
                                                                                    tatTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGt                                                            >  2:1066639/1‑150 (MQ=255)
                                                                                     atTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGtt                                                           <  2:872826/150‑1 (MQ=255)
                                                                                     atTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGtt                                                           >  2:454449/1‑150 (MQ=255)
                                                                                      ttCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATca                                                                                                                 >  2:1268766/1‑95 (MQ=255)
                                                                                      ttCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATca                                                                                                                 <  1:1268766/95‑1 (MQ=255)
                                                                                          gCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAt                                                      >  1:873062/1‑150 (MQ=255)
                                                                                             aGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCATGAGGTTTCTAtgtg                                                   <  2:710904/150‑1 (MQ=255)
                                                                                             aGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtg                                                     <  2:1139413/148‑1 (MQ=255)
                                                                                             aGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtg                                                     >  1:1139413/1‑148 (MQ=255)
                                                                                                cGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTgg                                                <  2:1686826/150‑1 (MQ=255)
                                                                                                cGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTgg                                                <  2:1486055/150‑1 (MQ=255)
                                                                                                    ccGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGc                                            >  1:1470169/1‑150 (MQ=255)
                                                                                                      gACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCtt                                          <  1:560518/150‑1 (MQ=255)
                                                                                                          aTTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGcc                                      >  2:381289/1‑150 (MQ=255)
                                                                                                           ttCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGATCATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGccc                                     <  1:1836624/150‑1 (MQ=255)
                                                                                                           ttCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGccc                                     <  2:719814/150‑1 (MQ=255)
                                                                                                               ttATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATc                                 >  1:714979/1‑150 (MQ=255)
                                                                                                                tATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATcc                                <  2:88868/150‑1 (MQ=255)
                                                                                                                 aTTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATccc                               <  1:699118/150‑1 (MQ=255)
                                                                                                                 aTTTAAATCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGAGGTAACATCGTTAAATTTGTATCACAAAGATTTTAACCGCCCCCACTACTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCAAATccc                               >  1:1365032/1‑150 (MQ=255)
                                                                                                                  tttAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCa                              <  1:1211205/150‑1 (MQ=255)
                                                                                                                   ttAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATGTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCaa                             >  2:1703144/1‑150 (MQ=255)
                                                                                                                      aTTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATATCAGGAGGTTTCTATGTGTGGACGCTTTGACCAATCCCAAAcg                          >  2:1109432/1‑150 (MQ=255)
                                                                                                                       ttCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAAcgc                         >  2:1105170/1‑150 (MQ=255)
                                                                                                                        tCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAAcgcg                        >  2:240306/1‑150 (MQ=255)
                                                                                                                          aGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTg                      <  1:1066639/150‑1 (MQ=255)
                                                                                                                           gAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGa                     >  2:383460/1‑150 (MQ=255)
                                                                                                                           gAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGa                     >  1:668771/1‑150 (MQ=255)
                                                                                                                              gTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGa                  >  1:843514/1‑150 (MQ=255)
                                                                                                                               tAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTGCCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAgga                 >  1:39146/1‑148 (MQ=255)
                                                                                                                                aaTCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGAAGCTTTGCCCAATCACAAACGCGTGAAGCtt                >  2:1052725/1‑150 (MQ=255)
                                                                                                                                  taccTCATTCCCGGCGGGCCCCACACCGGGTACATCGGTACCTTTGTATCACATACATTTTTACCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTAc              <  2:46160/148‑1 (MQ=255)
                                                                                                                                    ccTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCt            >  1:1953062/1‑150 (MQ=255)
                                                                                                                                       cATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTgc         <  2:668771/150‑1 (MQ=255)
                                                                                                                                         ttCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCTGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTgcgc       <  1:1133644/150‑1 (MQ=255)
                                                                                                                                         ttCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTgcgc       <  2:1027697/150‑1 (MQ=255)
                                                                                                                                          tCCCGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCt      >  2:1496150/1‑150 (MQ=255)
                                                                                                                                           cccTCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCtt     >  2:530751/1‑150 (MQ=255)
                                                                                                                                           cccGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCtt     >  2:803858/1‑150 (MQ=255)
                                                                                                                                           cccGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCtt     >  2:1817386/1‑150 (MQ=255)
                                                                                                                                              gCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTtctc  <  1:722704/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
AATCCTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTC  >  NC_000913/2009583‑2009874

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