Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,522,839 T→A 6.2% intergenic (+25/+15) PP_4856 → / ← yhjG Dps family ferritin/protein YhjG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,522,8390TA 6.2% 296.8 / 8.0 97intergenic (+25/+15)PP_4856/yhjGDps family ferritin/protein YhjG
Reads supporting (aligned to +/- strand):  ref base T (52/39);  new base A (4/2);  total (56/41)
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.06e-01

Only 100 of 101 total aligned reads displayed.

GCTACCGCGAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCAC  >  NC_002947/5522690‑5522988
                                                                                                                                                     |                                                                                                                                                     
gCTACCGCGAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCt                                                                                                                                                       <  1:578676/150‑1 (MQ=255)
          caTCATTTAATACATCGGGAATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGccc                                                                                                                                             <  2:870340/149‑1 (MQ=255)
             cATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCtt                                                                                                                                          <  2:1787717/150‑1 (MQ=255)
               ttCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGc                                                                                                                                        >  2:2029437/1‑150 (MQ=255)
               ttCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCAAAAGCGGCCCCtt                                                                                                                                          <  2:179067/148‑1 (MQ=255)
               ttCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCAAAAGCGGCCCCtt                                                                                                                                          >  1:179067/1‑148 (MQ=255)
               ttCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCAAAAGCGGCCCCTTga                                                                                                                                        >  2:415390/1‑149 (MQ=255)
               ttCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGATTTTGCGGCCCCTTGc                                                                                                                                        >  1:199281/1‑150 (MQ=255)
                  aaTACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTAt                                                                                                                                     >  1:926007/1‑150 (MQ=255)
                   aTACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTAtt                                                                                                                                    >  1:2136248/1‑150 (MQ=255)
                    tACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTGttt                                                                                                                                   <  2:1200622/150‑1 (MQ=255)
                      cATCGTGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTcc                                                                                                                                 <  2:1850956/150‑1 (MQ=255)
                        tCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGt                                                                                                                               <  2:1423189/150‑1 (MQ=255)
                        tCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGt                                                                                                                               >  1:1898225/1‑150 (MQ=255)
                               tAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCa                                                                                                                        <  2:199075/150‑1 (MQ=255)
                                      ccGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGt                                                                                                                 >  2:53640/1‑150 (MQ=255)
                                      ccGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGt                                                                                                                 >  1:543822/1‑150 (MQ=255)
                                           tACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCTTGCTTTCACCGCCGTCGgcg                                                                                                            >  1:1299090/1‑150 (MQ=255)
                                             cccGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGcc                                                                                                          <  1:175913/150‑1 (MQ=255)
                                             cccGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGcc                                                                                                          >  1:1941206/1‑150 (MQ=255)
                                             cccGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGcc                                                                                                          >  1:1971506/1‑150 (MQ=255)
                                                  cgtATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTtgcc                                                                                                     <  1:604097/150‑1 (MQ=255)
                                                   gtATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTtgcct                                                                                                    >  2:1143626/1‑150 (MQ=255)
                                                      tCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTtgcctgcc                                                                                                 <  2:543822/150‑1 (MQ=255)
                                                      tCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTtgcctgcc                                                                                                 <  2:1898225/150‑1 (MQ=255)
                                                       cTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTTCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTgcctgcct                                                                                                <  2:337102/150‑1 (MQ=255)
                                                       cTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTgcctgcct                                                                                                <  2:313930/150‑1 (MQ=255)
                                                              gACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATctg                                                                                         <  2:188538/150‑1 (MQ=255)
                                                              gACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACAACCTAAAGGCGCCGCTTTTGCGGCCCCTTGCTATTTCCGTTCTTTCACCTCCGTCGGCGCCTTGCCTGCCTTCATctg                                                                                         <  2:882649/150‑1 (MQ=255)
                                                                cATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCCTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATctgct                                                                                       <  1:53640/150‑1 (MQ=255)
                                                                 aTCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATctgctg                                                                                      >  1:1219929/1‑150 (MQ=255)
                                                                 aTCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATccgctg                                                                                      <  2:1239032/150‑1 (MQ=255)
                                                                  tCCTGGCGCAGGAAGAGGAGCACGCGGNCGATATGTCCGANCTNTTGNAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCtgctgc                                                                                     >  1:2128870/1‑150 (MQ=255)
                                                                  tCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCtgctgc                                                                                     >  1:2367892/1‑150 (MQ=255)
                                                                    cTGGCGGAGGAAGAGGAGCACGCGGACGATATGTCCGACGTGTGGCAAGGCCTTTGATTTCTGACATCCTCTAGGGACCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCTGTCGGCACCATGCCTGCTTTCATATGCTGcaa                                                                                   >  2:780221/1‑150 (MQ=255)
                                                                    cTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATctgctg                                                                                      >  1:959179/1‑147 (MQ=255)
                                                                     tGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGcaac                                                                                  >  2:2005222/1‑150 (MQ=255)
                                                                         gcAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGAGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACg                                                                              >  2:1179878/1‑150 (MQ=255)
                                                                         gcAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCNTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTTCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTCCAACAACg                                                                              >  1:862393/1‑150 (MQ=255)
                                                                                 aggagCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTtgcctgc                                                                                                  >  2:1099941/1‑122 (MQ=255)
                                                                                 aggagCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTtgcctgc                                                                                                  <  1:1099941/122‑1 (MQ=255)
                                                                                 aggagCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCAc                                                                      >  1:2183938/1‑150 (MQ=255)
                                                                                  ggagCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACt                                                                     >  1:1351135/1‑150 (MQ=255)
                                                                                    agCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGGCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTtgcct                                                                                                    >  1:1143510/1‑117 (MQ=255)
                                                                                    agCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTtgcct                                                                                                    <  2:1143510/117‑1 (MQ=255)
                                                                                    agCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTgg                                                                   >  2:376779/1‑150 (MQ=255)
                                                                                     gCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGt                                                                  >  1:215201/1‑150 (MQ=255)
                                                                                        cgcgGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGcaacaa                                                                                <  2:2075952/133‑1 (MQ=255)
                                                                                        cgcgGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGcaacaa                                                                                >  1:2075952/1‑133 (MQ=255)
                                                                                         gcgGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCgg                                                              <  1:9462/150‑1 (MQ=255)
                                                                                           ggACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGt                                                            >  1:2346042/1‑150 (MQ=255)
                                                                                           ccaCGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGt                                                            >  1:84497/3‑150 (MQ=255)
                                                                                              cGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCAt                                                         <  1:997046/150‑1 (MQ=255)
                                                                                                 tatGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTgcctgcct                                                                                                >  2:220693/1‑108 (MQ=255)
                                                                                                 tatGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTgcctgcct                                                                                                <  1:220693/108‑1 (MQ=255)
                                                                                                 tatGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATccc                                                      <  2:2367892/150‑1 (MQ=255)
                                                                                                     tCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCGGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCgct                                                  >  1:738866/1‑150 (MQ=255)
                                                                                                      ccGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCgctg                                                 <  1:376779/150‑1 (MQ=255)
                                                                                                       cGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCAAAAGCGGCCCCTTgagg                                                                                                                                      <  1:1293442/64‑4 (MQ=18)
                                                                                                       cGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCAAAAGCGGCCCCTTgagg                                                                                                                                      >  2:1293442/1‑61 (MQ=18)
                                                                                                         aCCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGctgctg                                              <  2:1971506/150‑1 (MQ=255)
                                                                                                           cTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTggg                                            <  2:2346042/150‑1 (MQ=255)
                                                                                                             gTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGgcg                                          >  1:2329738/1‑150 (MQ=255)
                                                                                                                gCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCATTAATCTGGTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCaa                                       >  2:64550/1‑150 (MQ=255)
                                                                                                                 cAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAAt                                      >  2:120635/1‑150 (MQ=255)
                                                                                                                  aaGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATc                                     >  2:2231337/1‑150 (MQ=255)
                                                                                                                     gCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAgc                                  <  1:2005222/150‑1 (MQ=255)
                                                                                                                       cTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAgcgc                                <  2:2183938/150‑1 (MQ=255)
                                                                                                                        tttGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGt                                                            >  1:497319/1‑121 (MQ=255)
                                                                                                                        tttGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGt                                                            <  2:497319/121‑1 (MQ=255)
                                                                                                                        tttGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTgg                                             <  2:2309117/136‑1 (MQ=255)
                                                                                                                        tttGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTgg                                             >  1:2309117/1‑136 (MQ=255)
                                                                                                                        tttGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGcc                               >  1:1124606/1‑150 (MQ=255)
                                                                                                                         ttGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGgcg                                          <  1:2063755/138‑1 (MQ=255)
                                                                                                                         ttGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGgcg                                          >  2:2063755/1‑138 (MQ=255)
                                                                                                                             tttCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCACAAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCaa                          >  1:1836153/1‑150 (MQ=255)
                                                                                                                                cTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCggg                                                             >  1:2264996/1‑112 (MQ=255)
                                                                                                                                cTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCggg                                                             <  2:2264996/112‑1 (MQ=255)
                                                                                                                                 tGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAAcccc                      <  2:1351135/150‑1 (MQ=255)
                                                                                                                                  gACACCCTCAAGGGGCCGCTTTTGCGGCCCGTGGCTATTTCCGTCCCTTCACCCCCCCCGGCGCCTTGCCTGCGTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCGTCCCCGCTGCTGGGCGCAATCAGCCGCCGCAACCccg                     >  2:1541621/1‑150 (MQ=255)
                                                                                                                                        cTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCccgccgccg               >  1:797541/1‑150 (MQ=255)
                                                                                                                                         tCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCgg              >  1:2323170/1‑150 (MQ=255)
                                                                                                                                         tCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCgg              >  1:2366373/1‑150 (MQ=255)
                                                                                                                                          cAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCAccgccgccg                                                                                                               >  2:1780837/1‑52 (MQ=255)
                                                                                                                                          cAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCg                                                                                                               <  1:1780837/52‑1 (MQ=255)
                                                                                                                                           aaGGGGCCGCAAAAGCGGCCCCTTGCTAGTTGCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCGGCAACAACGGTGCGTCATCGATCGTGGCATCCCCGCTGCTGGGCGACATCATGGACGGCATCGGCGCCGCCGGcc            >  1:1488632/1‑150 (MQ=255)
                                                                                                                                             ggggCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCt          >  2:314954/1‑150 (MQ=255)
                                                                                                                                             ggggCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCt          >  2:1321780/1‑150 (MQ=255)
                                                                                                                                             ggggCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCt          <  2:926007/150‑1 (MQ=255)
                                                                                                                                               ggCCGCTTTTGCGGCCCCTTGCTATTTCCTTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCACTGc        <  2:1299090/150‑1 (MQ=255)
                                                                                                                                               ggCCGATTTTGAGGCCCATCAGTCTGTACGTGTTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGc        <  2:2323170/150‑1 (MQ=255)
                                                                                                                                                gCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGcc       >  2:1765673/1‑150 (MQ=255)
                                                                                                                                                gCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGcc       >  2:2335377/1‑150 (MQ=255)
                                                                                                                                                 ccGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGcca      >  1:801302/1‑150 (MQ=255)
                                                                                                                                                 ccGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGcca      >  1:1350541/1‑150 (MQ=255)
                                                                                                                                                   gCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCa                              >  1:1540554/1‑124 (MQ=255)
                                                                                                                                                   gCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCa                              <  2:1540554/124‑1 (MQ=255)
                                                                                                                                                     ttttGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCcaccac  >  1:562723/1‑150 (MQ=255)
                                                                                                                                                     ttttGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCcaccac  <  1:1541621/150‑1 (MQ=255)
                                                                                                                                                     ttttGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCcaccac  >  1:1492751/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTACCGCGAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCAC  >  NC_002947/5522690‑5522988

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