Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,788,444:1 +A 7.8% coding (724/1563 nt) ahpF → alkyl hydroperoxide reductase subunit F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,788,4441.A 7.8% 390.9 / 28.9 103coding (724/1563 nt)ahpFalkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base . (45/50);  new base A (4/4);  total (49/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

Only 100 of 105 total aligned reads displayed.

GCAAGATCGACACCAACGCTGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAAGAGCACGTTAAACAGTACGACGTCGACATCATGAACCTGCAGCGCGGT  >  NC_002947/2788295‑2788587
                                                                                                                                                      |                                                                                                                                               
gCAAGATCGACACCAACGCTGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGccgccg                                                                                                                                                  >  2:622458/1‑150 (MQ=255)
  aaGATCGACACCAACGCTGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAg                                                                                                                                               >  1:844789/1‑150 (MQ=255)
  aaGATCGACACCAACGCTGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTGGCCGCCGaa                                                                                                                                               >  1:2057214/1‑149 (MQ=255)
           acCAACGCTGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTcggc                                                                                                                                      <  1:996154/150‑1 (MQ=255)
           acCAACGCTGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTcggc                                                                                                                                      <  1:622458/150‑1 (MQ=255)
           acCAACGCTGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTcggc                                                                                                                                      >  1:374379/1‑150 (MQ=255)
             cAACGCTGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAAGCGCTcggcgg                                                                                                                                    <  2:1671603/150‑1 (MQ=255)
             cAACGCTGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAAGCGCTcggcgg                                                                                                                                    >  1:1671603/1‑150 (MQ=255)
             cAACGCTGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTcggcgg                                                                                                                                    >  1:1692687/1‑150 (MQ=255)
              aaCGCTGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCggcggc                                                                                                                                   <  1:1376767/150‑1 (MQ=255)
               aCGCTGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGcc                                                                                                                                  >  1:1621882/1‑150 (MQ=255)
                 gCTGGCGCCCGCCAGGGCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCGGCCATCTACGCTGCCCGCAAGGGCAGCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAg                                                                                                                                <  1:1019592/150‑1 (MQ=255)
                 gCTGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAg                                                                                                                                <  2:1519805/150‑1 (MQ=255)
                  cTGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAgg                                                                                                                               <  1:107050/150‑1 (MQ=255)
                   tGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAAGCGCTcggcgg                                                                                                                                    >  1:772885/1‑144 (MQ=255)
                   tGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAAGCGCTcggcgg                                                                                                                                    <  2:772885/144‑1 (MQ=255)
                   tGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGt                                                                                                                              <  2:1273296/150‑1 (MQ=255)
                   tGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTGGCCGCCGAAGCGCTCggcggcaacac                                                                                                                              <  2:2057214/150‑6 (MQ=255)
                   tGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCAGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGt                                                                                                                              <  1:377887/150‑1 (MQ=255)
                      cgcccgccAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCt                                                                                                                           <  1:1972742/150‑1 (MQ=255)
                           gccAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGaca                                                                                                                      <  2:1121670/150‑1 (MQ=255)
                                gCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTg                                                                                                                 >  2:1025559/1‑150 (MQ=255)
                                gCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTg                                                                                                                 <  2:1931901/150‑1 (MQ=255)
                                   gagaAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGcc                                                                                                              <  2:764052/150‑1 (MQ=255)
                                     gaagaTCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCAt                                                                                                            >  2:1657297/1‑150 (MQ=255)
                                     gaagaTCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCAt                                                                                                            <  2:1715910/150‑1 (MQ=255)
                                      aagaTCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATc                                                                                                           <  2:365601/150‑1 (MQ=255)
                                        gagcAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCga                                                                                                         >  1:82210/4‑150 (MQ=255)
                                        gaTCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGAATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCga                                                                                                         <  2:414485/150‑1 (MQ=255)
                                        gaTCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGAATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCga                                                                                                         >  1:414485/1‑150 (MQ=255)
                                        gaTCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGCGCATCCGCACCTGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCga                                                                                                         <  2:498884/150‑1 (MQ=255)
                                         aTCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCgag                                                                                                        <  2:1692687/150‑1 (MQ=255)
                                           cAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGaa                                                                                                      <  1:86245/150‑1 (MQ=255)
                                           cAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGaa                                                                                                      <  2:445170/150‑1 (MQ=255)
                                           cAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGaa                                                                                                      <  1:2312624/150‑1 (MQ=255)
                                           cAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGaa                                                                                                      <  2:182216/150‑1 (MQ=255)
                                               gCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGGGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTc                                                                                                  <  2:374379/150‑1 (MQ=255)
                                               gCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAAGCGCTCggcggcaacac                                                                                                                              >  2:547741/1‑117 (MQ=255)
                                               gCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAAGCGCTCggcggcaacac                                                                                                                              <  1:547741/122‑6 (MQ=255)
                                               gCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTc                                                                                                  >  1:311023/1‑150 (MQ=255)
                                                     gAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCg                                                                                            >  1:2057527/1‑150 (MQ=255)
                                                       aGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCGGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGt                                                                                          >  1:1402503/1‑150 (MQ=255)
                                                          cTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCg                                                                                                          >  2:685419/1‑131 (MQ=255)
                                                          cTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCg                                                                                                          <  1:685419/131‑1 (MQ=255)
                                                          cTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCACCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCa                                                                                       >  1:1170000/1‑150 (MQ=255)
                                                           tttGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAg                                                                                      <  1:1025559/150‑1 (MQ=255)
                                                           tttGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAg                                                                                      >  1:34631/1‑150 (MQ=255)
                                                              gACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGaa                                                                                   <  2:1621882/150‑1 (MQ=255)
                                                                cGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAAc                                                                                 >  1:1724864/1‑150 (MQ=255)
                                                                     tGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCgg                                                                                           >  1:215745/1‑135 (MQ=255)
                                                                     tGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCgg                                                                                           <  2:215745/135‑1 (MQ=255)
                                                                       gtcgtcGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTcg                                                                                                                                        <  2:1850695/88‑1 (MQ=255)
                                                                       gtcgtcGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTcg                                                                                                                                        >  1:1850695/1‑88 (MQ=255)
                                                                       gtcgtcGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAggg                                                                          >  2:1222217/1‑150 (MQ=255)
                                                                          gtcgGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCg                                                                       >  1:763278/1‑150 (MQ=255)
                                                                            cggcggTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGaa                                                                     <  2:82210/150‑1 (MQ=255)
                                                                                ggTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTg                                                                 <  2:1403906/150‑1 (MQ=255)
                                                                                 gTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCACCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTgg                                                                <  2:1170000/150‑1 (MQ=255)
                                                                                  tGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGc                                                               >  1:937845/1‑150 (MQ=255)
                                                                                     ccccGCTGGCGCCGCTGCTGCCATCTAAGCTGCCAGAATAGGAATGCGCCTTTTTACTTCCTGCGAGCTCTTCTGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCAc                                                            >  1:1797929/1‑150 (MQ=255)
                                                                                       ccGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCga                                                                                                         >  1:2207941/1‑103 (MQ=255)
                                                                                       ccGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCga                                                                                                         <  2:2207941/103‑1 (MQ=255)
                                                                                       ccGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACgg                                                          <  2:34631/150‑1 (MQ=255)
                                                                                       ccGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACgg                                                          <  1:417341/150‑1 (MQ=255)
                                                                                       ccGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACgg                                                          >  1:1299051/1‑150 (MQ=255)
                                                                                        cGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAgg                                                                                                                               <  2:1433360/80‑1 (MQ=255)
                                                                                        cGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAgg                                                                                                                               >  1:1433360/1‑80 (MQ=255)
                                                                                        cGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGgc                                                         >  1:1263023/1‑150 (MQ=255)
                                                                                               gccgcTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGaa                                                  >  1:777341/1‑150 (MQ=255)
                                                                                                 cgcTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAaga                                                <  2:311023/150‑1 (MQ=255)
                                                                                                 cgcTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAaga                                                <  1:117781/150‑1 (MQ=255)
                                                                                                      ctgcCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAAGAGCACg                                           <  2:1263023/150‑1 (MQ=255)
                                                                                                        cccATCTTCCCTCCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTTGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAAGAGCACGtt                                         <  1:979054/149‑1 (MQ=255)
                                                                                                         ccATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCga                                                                                                         <  2:1681128/85‑1 (MQ=255)
                                                                                                         ccATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCga                                                                                                         >  1:1681128/1‑85 (MQ=255)
                                                                                                          cATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAAGAGCACGTTaa                                       >  2:1356845/1‑150 (MQ=255)
                                                                                                           aTCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGc                                                                                                               >  2:1080297/1‑77 (MQ=255)
                                                                                                           aTCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGc                                                                                                               <  1:1080297/77‑1 (MQ=255)
                                                                                                           aTCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAaga                                                <  2:28763/140‑1 (MQ=255)
                                                                                                           aTCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAaga                                                >  1:28763/1‑140 (MQ=255)
                                                                                                           aTCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAAGAGCACGTTaaa                                      <  2:1724864/150‑1 (MQ=255)
                                                                                                           aTCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAAGAGCACGTTaaa                                      >  1:1073494/1‑150 (MQ=255)
                                                                                                              tACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTgg                                                                <  1:1760883/121‑1 (MQ=255)
                                                                                                              tACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTgg                                                                >  2:1760883/1‑121 (MQ=255)
                                                                                                                 gCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAAGAGCACGTTAAACAGTac                                >  2:2198261/1‑150 (MQ=255)
                                                                                                                      ccGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAAACGAAGGGCCGAAGCTGGCCAAGGCGCTGGAAGAGCACGTTAAACAGTacgacgg                           >  2:244978/1‑149 (MQ=255)
                                                                                                                      ccGCAAGGGCATACGCACCGGAGTTTCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAAGAGCACGTTAAACAGTACGACGt                           >  2:1959495/1‑150 (MQ=255)
                                                                                                                        gCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGGGTCTGGCCACGGCGCTGGAAGAGCACGTTAAACAGTACGACGTCg                         >  1:840901/1‑150 (MQ=255)
                                                                                                                        gCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAAGAGCACGTTAAACAGTACGACGTc                          <  1:1658527/149‑1 (MQ=255)
                                                                                                                        gCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAAGAGCACGTTAAACAGTACGACGTc                          >  2:1658527/1‑149 (MQ=255)
                                                                                                                        gCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAAGAGCACGTTAAACAGTACGACGTCg                         >  1:1586960/1‑150 (MQ=255)
                                                                                                                        gCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGAAGAGAGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGCAGGGCCGAAGCTGGCCACGGCGCTGGAAGAGCACGTTAAACAGTACGACGTCg                         >  1:1590502/1‑150 (MQ=255)
                                                                                                                          aaGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGc                                                                                        >  2:583097/1‑85 (MQ=255)
                                                                                                                                      aCCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGacac                                                                                                                     >  2:1163529/1‑44 (MQ=255)
                                                                                                                                      aCCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGacac                                                                                                                     <  1:1163529/44‑1 (MQ=255)
                                                                                                                                        cGGTGTTCCCGCCGACCCCTTCGGCGGCCAGGTGCTGGACCCGCGGGCCATCGAGAACTGCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAAGAGCACGTTAAACAGTACGACGTCGACATCATGAACCTGCa         <  2:763278/150‑1 (MQ=255)
                                                                                                                                           tgtTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAAGAGCACGTTAAACAGTACGACGTCGACATCATGAACCTGCAgcg      >  1:1893601/1‑150 (MQ=255)
                                                                                                                                             ttGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAAGAGCACGTTAAACAGTacga                              >  2:1999757/1‑124 (MQ=255)
                                                                                                                                             ttGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAAGAGCACGTTAAACAGTacga                              <  1:1999757/124‑1 (MQ=255)
                                                                                                                                               gccgccGAAGCGCTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAAGAGCACGTTAAACAGTACGACGTCGACATCATGAACCTGCAGCGCGGt  >  2:1185909/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
GCAAGATCGACACCAACGCTGGCGCCCGCCAGGCCGAGAAGATCAACGCCAAGGAAGCCTTTGACGTGCTGGTCGTCGGCGGTGGCCCCGCTGGCGCCGCTGCTGCCATCTACGCTGCCCGCAAGGGCATCCGCACCGGTGTTGCCGCCGAGCGCTTCGGCGGCCAGGTGCTGGACACTCTGGCCATCGAGAACTTCATCTCGGTGCAGGAAACCGAAGGGCCGAAGCTGGCCACGGCGCTGGAAGAGCACGTTAAACAGTACGACGTCGACATCATGAACCTGCAGCGCGGT  >  NC_002947/2788295‑2788587

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