Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,009,731 A→G 11.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,7310AG11.5% 230.4 / 21.4 96intergenic (+19/‑90)yedF/yedKputative TusA family sulfurtransferase/DUF159 family protein
Reads supporting (aligned to +/- strand):  ref base A (42/42);  new base G (2/9);  total (45/51)
Fisher's exact test for biased strand distribution p-value = 5.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.96e-01

TCCTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTCGCGG  >  NC_000913/2009585‑2009878
                                                                                                                                                  |                                                                                                                                                   
tCCTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGAccc                                                                                                                                                  <  1:684052/150‑1 (MQ=255)
     gAAGTGGTGAGCGACTGTCCGCAGTCGATCTATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCATGAGGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGc                                                                                                                                             >  1:273880/1‑150 (MQ=255)
      atgttgtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTTGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCg                                                                                                                                            <  1:1338829/145‑1 (MQ=255)
       aGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGGAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCgg                                                                                                                                           <  2:742706/150‑1 (MQ=255)
       aGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCgg                                                                                                                                           <  1:892959/150‑1 (MQ=255)
        gtggtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCggg                                                                                                                                          <  2:1582159/150‑1 (MQ=255)
        gtggtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCggg                                                                                                                                          <  1:1652481/150‑1 (MQ=255)
         tggtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGa                                                                                                                                         <  2:664695/150‑1 (MQ=255)
         tggtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGa                                                                                                                                         <  1:551887/150‑1 (MQ=255)
          ggtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGaa                                                                                                                                        >  1:1330259/1‑150 (MQ=255)
           gtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGaat                                                                                                                                       <  2:117314/150‑2 (MQ=255)
              aGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGaatgag                                                                                                                                    <  1:516918/150‑5 (MQ=255)
                cGACTGTCCGCAGTCGATGAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCAGTGACCCACAGCGGGAACATCGt                                                                                                                                  >  1:319785/1‑150 (MQ=255)
                cGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCCCGGCTATCCGGTGCTGGATATTCAGCAAGCCGGGCCGCCCATTCGGTATTTACTTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGATCATCGt                                                                                                                                  >  1:276028/1‑150 (MQ=255)
                cGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTCTACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGaatgaggg                                                                                                                                  <  2:144998/150‑7 (MQ=255)
                cGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGt                                                                                                                                  <  2:1171380/150‑1 (MQ=255)
                       ccGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCttt                                                                                                                           >  1:49894/1‑150 (MQ=255)
                       ccGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCttt                                                                                                                           <  2:412237/150‑1 (MQ=255)
                        cGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCgg                                                                                                                                           <  1:508009/133‑1 (MQ=255)
                        cGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCgg                                                                                                                                           >  2:508009/1‑133 (MQ=255)
                        cGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTg                                                                                                                          <  2:864860/150‑1 (MQ=255)
                        cGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTg                                                                                                                          <  2:532515/150‑1 (MQ=255)
                         gCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGGTATTTAATTCAGAAGTAAGCCCTCATTCCCGCTGTGACCCACAGCCGGAACATCGTTACCTTTGt                                                                                                                         >  2:293962/1‑150 (MQ=255)
                         gCAGACGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGCCGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCCCAGCGGCAACATCGTTACGTTTGt                                                                                                                         >  2:343680/1‑150 (MQ=255)
                          cattCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTa                                                                                                                        <  1:351963/147‑1 (MQ=255)
                          cAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTa                                                                                                                        <  1:1566233/150‑1 (MQ=255)
                           aGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTAt                                                                                                                       >  2:802597/1‑150 (MQ=255)
                               gATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcaca                                                                                                                   >  1:383200/1‑150 (MQ=255)
                               gATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcac                                                                                                                    >  2:1698711/1‑149 (MQ=255)
                               gATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcac                                                                                                                    <  1:1698711/149‑1 (MQ=255)
                                          ttCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTAcc                                                                                                        <  1:1060712/150‑1 (MQ=255)
                                          ttCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATCTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTAcc                                                                                                        <  1:150685/150‑1 (MQ=255)
                                          ttCCACTGGATGCGCGTAATCAAGTCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTAcc                                                                                                        <  2:1367865/150‑1 (MQ=255)
                                             cACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGc                                                                                                     >  2:1124580/1‑150 (MQ=255)
                                              aCTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcaca                                                                                                                   >  1:1155600/1‑135 (MQ=255)
                                              aCTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcaca                                                                                                                   <  2:1155600/135‑1 (MQ=255)
                                                        cgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGCCCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTa                                                                                          >  1:528086/1‑150 (MQ=255)
                                                         gTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTTATCCCTCATACCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTAt                                                                                         <  1:199427/150‑1 (MQ=255)
                                                           aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTc                                                                                       <  2:1677834/150‑1 (MQ=255)
                                                           aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTc                                                                                       <  2:1671478/150‑1 (MQ=255)
                                                            aTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCt                                                                                      <  2:189632/150‑1 (MQ=255)
                                                               aCGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCtttt                                                                                   >  2:522150/1‑150 (MQ=255)
                                                               aCGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCtttt                                                                                   >  1:1060626/1‑150 (MQ=255)
                                                                cGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTc                                                                                  <  1:1552006/150‑1 (MQ=255)
                                                                cGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTc                                                                                  >  1:243588/1‑150 (MQ=255)
                                                                  gCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAg                                                                                <  2:1452060/150‑1 (MQ=255)
                                                                  gCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAg                                                                                >  1:487426/1‑150 (MQ=255)
                                                                   cTATACGGTGCTGGATCTTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGt                                                                               <  2:273880/150‑1 (MQ=255)
                                                                   cTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGt                                                                               <  2:1330259/150‑1 (MQ=255)
                                                                    tataCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTa                                                                              <  2:694764/150‑1 (MQ=255)
                                                                        cGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGAt                                                                          >  2:1426072/1‑150 (MQ=255)
                                                                         ggTGCTGGATATTCAGCAAGACGGGCCGACGTTTCGTTATTTAATTCATAATTAATCCCTCATTCCCGTTTTGACCCACAGCGGGAACTTCGTTAGCTTTGTATCACATAGATTTTACCCGTCCCCACTCCTTTTCTTTTCAGTTTGTtc                                                                         <  1:1364918/150‑1 (MQ=255)
                                                                         ggTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGAtc                                                                         <  2:1658634/150‑1 (MQ=255)
                                                                         ggTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCATCATTCCCGCTGTGACCCACAGCGGGACCATCGTTAACTTTGTATCACATAGATTTTACCCGCCCCCAGACATATTCTTTTCATTATGAtt                                                                         >  1:261080/1‑149 (MQ=255)
                                                                              tGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCagg                                                                    <  2:464357/150‑1 (MQ=255)
                                                                                   atTCAGCAAGACGGGCCGACCAATCGTTATTTAATTCAGACGTAATCCCTCAAGCCCGCGGGGACCCACAGCGTGAACATAGTTACCTTTCTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGtt                                                               <  2:470045/150‑1 (MQ=255)
                                                                                    ttCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGttt                                                              >  1:1713585/1‑150 (MQ=255)
                                                                                      cAGCAAGACGGGCCGACCATTCGTGATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGACCATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCt                                                            >  2:1718724/1‑150 (MQ=255)
                                                                                        gCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAt                                                          >  1:1626185/1‑150 (MQ=255)
                                                                                         cAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtg                                                         <  2:440931/150‑1 (MQ=255)
                                                                                           aGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtgtg                                                       >  1:979735/1‑150 (MQ=255)
                                                                                           aGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtgtg                                                       >  2:130419/1‑150 (MQ=255)
                                                                                           aGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCAGCTGTGACCCACAGCGGGATCATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATCTTCACAGGAGGTTTCTAtgtg                                                       >  1:163315/1‑150 (MQ=255)
                                                                                           aGAAGGGCCGACCATTAGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGATCATCGTTACCTTTGTATCACTTAGATTTTACCCGCCCCCTCTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtgtg                                                       <  1:273672/150‑1 (MQ=255)
                                                                                            gACGGGCGGAACATTCGTTATTTAACTCAGAATTAATCCCTCATTCCCGATGTGACACACATCGGGAACATCGTTACCTTTGTATCAAATATATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtgtgt                                                      <  2:309562/150‑1 (MQ=255)
                                                                                             aCGGGCCGACCATTCGTAATTTAATTCCGAAGTAATCCCTCATTCCCGCGGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtgtgtg                                                     <  2:1369312/150‑1 (MQ=255)
                                                                                              cGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTgg                                                    >  1:1492127/1‑150 (MQ=255)
                                                                                                  ccGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGc                                                >  1:1299062/1‑150 (MQ=255)
                                                                                                    gACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCtt                                              >  2:973670/1‑150 (MQ=255)
                                                                                                       cATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACGCACAGCGGGAACATCGTTAGCTTTGTATCACATAGATTTTACGCGCCCCCTCTCCTATTCTTTTCAGTATTATCTCAGGAGGTTTCTATGTGTGGACGCTTTGc                                           <  1:157229/150‑1 (MQ=255)
                                                                                                       cATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGc                                           >  2:1149035/1‑150 (MQ=255)
                                                                                                         ttCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGccc                                         <  2:1060626/150‑1 (MQ=255)
                                                                                                             ttATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATc                                     <  2:49894/150‑1 (MQ=255)
                                                                                                              tATTTAATTCAGAAGTAATCCCTCAATCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATcc                                    <  2:276028/150‑1 (MQ=255)
                                                                                                               aTTTAATTCAGAAGTAATCCCTCCTTCCCGGAGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATATCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATccc                                   <  2:261080/150‑1 (MQ=255)
                                                                                                               aTTTAATTCAGAAGTAATCCCTCATTCACGCTGTGACCCACACCGGGAACATCTTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTATTTTCAGTATGATCTCAGGACGTTTCTATGTGTGGACGCTATGCCCAATccc                                   >  2:523606/1‑150 (MQ=255)
                                                                                                               aTTTAATTAAGAAGAAATCCCTCATTCCCTATGTGACCCACCGAGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCACCACACCTATTCTTTTCAGTATGATCAAAGGAGGTTTCTATGTGTGGACGATTTGACCAATccc                                   >  1:149080/1‑150 (MQ=255)
                                                                                                                tttAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCa                                  <  1:802597/150‑1 (MQ=255)
                                                                                                                tttAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCa                                  >  2:1523931/1‑150 (MQ=255)
                                                                                                                   aaTTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAAc                               <  1:664455/150‑1 (MQ=255)
                                                                                                                    aTTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAAcg                              >  1:1579555/1‑150 (MQ=255)
                                                                                                                    aTTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAAcg                              >  1:1021664/1‑150 (MQ=255)
                                                                                                                                tCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTAc                  <  2:979735/150‑1 (MQ=255)
                                                                                                                                tCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCATTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTAc                  >  2:1329735/1‑150 (MQ=255)
                                                                                                                                   cTCATTCCCGCTGTGACCCCCAGCGGGAACATCGTTACCTTTGTATGCCAGAGATTTTTCCCGCCCCCACTCCTATTCTTTTCAGTTTGAGCTCAGGAGGGTTCGATGTGTGGACGCTTTGCCCAATCGCAAACGCGTGAAGCTTACCtt               >  1:881838/1‑150 (MQ=255)
                                                                                                                                    tCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTg              <  2:383200/150‑1 (MQ=255)
                                                                                                                                       ttCCCGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCCCTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTgcgc           >  1:860147/1‑150 (MQ=255)
                                                                                                                                        tCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCt          >  2:1128914/1‑150 (MQ=255)
                                                                                                                                         cccGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACCTAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGAGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCtt         >  1:232724/1‑150 (MQ=255)
                                                                                                                                         cccGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCtt         >  2:1167715/1‑150 (MQ=255)
                                                                                                                                          ccGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTtc        >  2:1450199/1‑150 (MQ=255)
                                                                                                                                            gCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTtctc      <  1:1187754/150‑1 (MQ=255)
                                                                                                                                            gCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTtctc      >  2:909958/1‑150 (MQ=255)
                                                                                                                                             cTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTcg     >  2:1111375/1‑150 (MQ=255)
                                                                                                                                             cTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTcg     >  2:1117192/1‑150 (MQ=255)
                                                                                                                                                tgACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTCGCgg  <  1:973670/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCCTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTCGCGG  >  NC_000913/2009585‑2009878

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: