Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,527,788 T→A 5.4% 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,527,7880TA 5.4% 930.7 / 29.2 279intergenic (+25/+15)PP_4856/yhjGDps family ferritin/protein YhjG
Reads supporting (aligned to +/- strand):  ref base T (146/115);  new base A (6/9);  total (154/125)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.80e-01

Only 100 of 294 total aligned reads displayed.

CTACCGCGAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCAC  >  NC_002947/5527640‑5527937
                                                                                                                                                    |                                                                                                                                                     
cTACCGCGAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCtt                                                                                                                                                      >  2:1371293/1‑150 (MQ=255)
   ccGCGAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTg                                                                                                                                                   <  2:5560512/150‑1 (MQ=255)
    cgcgAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGc                                                                                                                                                  <  1:3725508/150‑1 (MQ=255)
    cgcgAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGc                                                                                                                                                  <  2:2320505/150‑1 (MQ=255)
         gaTCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGccc                                                                                                                                             >  2:1271684/1‑150 (MQ=255)
          aTCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGcccc                                                                                                                                            <  1:4468533/150‑1 (MQ=255)
             aTTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCtttt                                                                                                                                                    <  2:2811974/139‑1 (MQ=255)
             aTTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCtttt                                                                                                                                                    >  1:2811974/1‑139 (MQ=255)
              ttCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGc                                                                                                                                        >  1:5072875/1‑150 (MQ=255)
              ttCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCAAAAGCGGCCCCTTga                                                                                                                                        <  2:2780072/150‑2 (MQ=255)
              ttCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCAAAAGCGGCCCCTTga                                                                                                                                        <  1:2510559/150‑2 (MQ=255)
               tCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCAAAAGCGGCCCCTTgag                                                                                                                                       <  1:6484500/150‑3 (MQ=255)
                 aaTACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTAt                                                                                                                                     <  2:3975500/150‑1 (MQ=255)
                 aaTACATCGGGGATAAAGACCCGACTACCCGTCGTATCGTCGAAGCCATCCTGGCGCAGGAAGAGGCGCACGCGGACGAGATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCGTGTGCGGCCCCTGGCTAt                                                                                                                                     >  1:3404228/1‑150 (MQ=255)
                       tCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGt                                                                                                                               <  2:4053883/150‑1 (MQ=255)
                          gggATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCt                                                                                                                            >  1:5132296/1‑150 (MQ=255)
                           ggATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCtt                                                                                                                           >  1:933023/1‑150 (MQ=255)
                                     ccGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGt                                                                                                                 >  1:835979/1‑150 (MQ=255)
                                      cgacgaCCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCGATTTCCGTGCTTTCAccgcccag                                                                                                                <  2:2190257/150‑4 (MQ=255)
                                            cccGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGcc                                                                                                          >  2:687014/1‑150 (MQ=255)
                                              cgtcgtATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCtt                                                                                                        <  2:2904496/150‑1 (MQ=255)
                                                    aTCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTtgcctgc                                                                                                  <  1:972750/150‑1 (MQ=255)
                                                     tCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTtgcctgcc                                                                                                 <  2:4405337/150‑1 (MQ=255)
                                                      cTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTgcctgcct                                                                                                <  1:198915/150‑1 (MQ=255)
                                                       ttCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCtt                                                                                               >  1:2423216/1‑150 (MQ=255)
                                                         cGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCa                                                                                             <  2:5601945/150‑1 (MQ=255)
                                                         cGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCa                                                                                             <  2:376658/150‑1 (MQ=255)
                                                          gAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCAt                                                                                            >  2:1348992/1‑150 (MQ=255)
                                                          gAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCAt                                                                                            <  1:4207227/150‑1 (MQ=255)
                                                           aaGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGgcgc                                                                                                           >  1:611727/1‑134 (MQ=255)
                                                            aGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATct                                                                                          >  2:3497719/1‑150 (MQ=255)
                                                             caCAGCCCGGCGCCGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATctg                                                                                         <  1:7310227/149‑1 (MQ=255)
                                                               cATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATctgct                                                                                       <  2:6323449/150‑1 (MQ=255)
                                                                aTCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATctgctg                                                                                      <  1:4572166/150‑1 (MQ=255)
                                                                 tCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCtgctgc                                                                                     <  1:4030675/150‑1 (MQ=255)
                                                                    tGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGcaac                                                                                  <  2:3167373/150‑1 (MQ=255)
                                                                      gcgcAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATATGCTGcaacaa                                                                                >  1:702263/1‑150 (MQ=255)
                                                                       cgcAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCaacaac                                                                               <  1:1857235/150‑1 (MQ=255)
                                                                        gcAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACg                                                                              <  2:372822/150‑1 (MQ=255)
                                                                        gcAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACg                                                                              >  2:359656/1‑150 (MQ=255)
                                                                          aGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGgt                                                                            <  1:2874741/150‑1 (MQ=255)
                                                                          aGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGgt                                                                            >  2:3425117/1‑150 (MQ=255)
                                                                            gAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGgtgt                                                                          <  2:589681/150‑1 (MQ=255)
                                                                            gAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGgtgt                                                                          <  2:660930/150‑1 (MQ=255)
                                                                            gAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGgtgt                                                                          >  1:1062681/1‑150 (MQ=255)
                                                                            gAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGgtgt                                                                          >  1:295000/1‑150 (MQ=255)
                                                                              agagGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGc                                                                        >  2:6182646/1‑150 (MQ=255)
                                                                               gaggagCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCa                                                                       <  1:6429049/150‑1 (MQ=255)
                                                                               gaggagCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCa                                                                       <  1:6171401/150‑1 (MQ=255)
                                                                                  gagCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTg                                                                    <  2:4892254/150‑1 (MQ=255)
                                                                                        gcgGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCgg                                                              >  1:1109364/1‑150 (MQ=255)
                                                                                         cgGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCggg                                                             <  2:1542964/150‑1 (MQ=255)
                                                                                          ggACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGt                                                            <  2:2056706/150‑1 (MQ=255)
                                                                                              gATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATc                                                        >  2:2437856/1‑150 (MQ=255)
                                                                                              gATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATc                                                        >  1:7155020/1‑150 (MQ=255)
                                                                                                tatGTCCGACATGTTGCAAGGCCTATGATTTCAGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATccc                                                      >  2:7226983/1‑150 (MQ=255)
                                                                                                 atGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATcccc                                                     <  1:1027148/150‑1 (MQ=255)
                                                                                                  tGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCg                                                    >  1:6554144/1‑150 (MQ=255)
                                                                                                   gTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCtttt                                                                                                                                                    <  2:5109310/53‑1 (MQ=255)
                                                                                                   gTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCgc                                                   >  1:2265513/1‑150 (MQ=255)
                                                                                                        aCCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGctgctg                                              <  2:4253193/150‑1 (MQ=255)
                                                                                                                 aaGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATc                                     >  2:4521805/1‑150 (MQ=255)
                                                                                                                   ggCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAg                                   >  1:1631965/1‑150 (MQ=255)
                                                                                                                    gCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAgc                                  <  2:2245246/150‑1 (MQ=255)
                                                                                                                    gCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAgc                                  <  2:900020/150‑1 (MQ=255)
                                                                                                                       tttGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCCATCAGCGcc                               >  1:5266369/1‑150 (MQ=255)
                                                                                                                       tttGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGcc                               >  2:4639057/1‑150 (MQ=255)
                                                                                                                       tttGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGcc                               >  1:6682235/1‑150 (MQ=255)
                                                                                                                       tttGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGcc                               >  2:6412443/1‑150 (MQ=255)
                                                                                                                       tttGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGcc                               >  1:7009727/1‑150 (MQ=255)
                                                                                                                       tttGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGcc                               >  1:5428149/1‑150 (MQ=255)
                                                                                                                       tttGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGcc                               >  1:3436801/1‑150 (MQ=255)
                                                                                                                       tttGATTTCTGACACACTCAAGGGGCCGCTTTTTCGGCCCCTTGCTATTTCCGTTCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATACCCGCAGCTGGGCGCAATCAGCGcc                               >  2:6255132/1‑150 (MQ=255)
                                                                                                                        ttGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGTGCGCAATCAGCGCCa                              <  1:4349624/150‑1 (MQ=255)
                                                                                                                        ttGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCa                              >  2:7272326/1‑150 (MQ=255)
                                                                                                                        ttGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCa                              <  2:3884252/150‑1 (MQ=255)
                                                                                                                        ttGATTTCTGACACCCTCAAGGGGCCGCCTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCa                              >  1:4281739/1‑150 (MQ=255)
                                                                                                                           aTTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCa                           >  2:6982176/1‑150 (MQ=255)
                                                                                                                            tttCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCCGCGCCAGCac                          >  1:2093785/1‑149 (MQ=255)
                                                                                                                            tttCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCaa                          >  1:3692514/1‑150 (MQ=255)
                                                                                                                            tttCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCaa                          >  2:5795810/1‑150 (MQ=255)
                                                                                                                             ttCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAAc                         <  2:1737420/150‑1 (MQ=255)
                                                                                                                             ttCGGACCCCCTCAAGGGGCCGCGTGTGCGGCCCCGTGCTAGGTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCCTCTGCTGCAACACCGGTGTGCCCTGGGTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAAc                         <  2:3404228/150‑1 (MQ=255)
                                                                                                                                tGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAAcccc                      >  1:3159837/1‑150 (MQ=255)
                                                                                                                                  acacCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTc                                                                                              <  1:1319010/76‑1 (MQ=255)
                                                                                                                                  acacCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCATGCCTGCCTTCATCTGCTGCAACAACGCTGTGCACAGGTTCGGGTCAACCCCGCAGCAGGGCGCAATCACCGCCCGCAAaccccc                    >  1:2639032/1‑148 (MQ=255)
                                                                                                                                  acacCCTCAAGGGGCCGCTTTTGCGGCCCCTTGATATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCccgc                    <  1:6344896/150‑1 (MQ=255)
                                                                                                                                    acCCTCAAGGGGCCGCTTTTGCGGCCCCGTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACAGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACccccccc                  >  1:3094687/1‑149 (MQ=255)
                                                                                                                                       cTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCccgccgccg               >  2:1968321/1‑150 (MQ=255)
                                                                                                                                       cTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCccgccgccg               >  1:1528050/1‑150 (MQ=255)
                                                                                                                                       cTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCccgccgccg               >  2:4784667/1‑150 (MQ=255)
                                                                                                                                          aaGGGGCCGCAAAAGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCa                              >  2:6081973/1‑132 (MQ=255)
                                                                                                                                          aaGGGGCCGCAAAAGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCa                              <  1:6081973/132‑1 (MQ=255)
                                                                                                                                           aGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGccc           >  1:4843310/1‑150 (MQ=255)
                                                                                                                                            ggggCCGCTTTTGCGGCCCCTCGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCt          >  2:831984/1‑150 (MQ=255)
                                                                                                                                               gCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCAGcc       >  1:3150366/1‑150 (MQ=255)
                                                                                                                                                 cGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGccac     >  2:5471222/1‑150 (MQ=255)
                                                                                                                                                    ttttGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCcaccac  >  1:4675389/1‑150 (MQ=255)
                                                                                                                                                    ttttGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCcaccac  <  1:561254/150‑1 (MQ=255)
                                                                                                                                                    ttttGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCcaccac  >  1:3834955/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTACCGCGAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTGCCACCAC  >  NC_002947/5527640‑5527937

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