Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,707,781 G→A 6.8% R491C (CGC→TGC)  cysI ← sulfite reductase hemoprotein subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,707,7810GA 6.8% 295.7 / 7.1 101R491C (CGC→TGC) cysIsulfite reductase hemoprotein subunit beta
Reads supporting (aligned to +/- strand):  ref base G (43/51);  new base A (2/5);  total (45/56)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

GGGTCGATGCCGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGCAGCCGGAGATGTTCAGGTCGATTTCGCCGATGTCGAACAG  >  NC_002947/2707633‑2707925
                                                                                                                                                    |                                                                                                                                                
gggTCGATGCCGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGt                                                                                                                                                 >  2:86347/1‑150 (MQ=255)
  gTCGATGCCGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTgc                                                                                                                                               >  2:1585233/1‑150 (MQ=255)
  gTCGATGCCGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTgc                                                                                                                                               >  2:619998/1‑150 (MQ=255)
     gATGCCGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTgcggc                                                                                                                                            >  2:1896951/1‑150 (MQ=255)
     gATGCCGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGAGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTgcggc                                                                                                                                            <  1:1232594/150‑1 (MQ=255)
         ccGATACGCTGGTAGGTGTCGATGAAACGCTTTTGCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCACgcggcgcaa                                                                                                                                        <  2:1245699/150‑4 (MQ=255)
         ccGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCACgcggcgcaa                                                                                                                                        <  2:2034831/150‑4 (MQ=255)
         ccGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCACgcggcgcaa                                                                                                                                        <  2:1565474/150‑4 (MQ=255)
          cGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCACgcggcgcaa                                                                                                                                        >  2:237062/1‑146 (MQ=255)
          cGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCACgcggcgcaa                                                                                                                                        <  1:237062/149‑4 (MQ=255)
           gATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACAAGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTgcc                                                                                                                                      >  1:100403/1‑150 (MQ=255)
                gCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCa                                                                                                                                 <  2:1658238/150‑1 (MQ=255)
                  tGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGc                                                                                                                               <  2:1164414/150‑1 (MQ=255)
                  tGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGc                                                                                                                               >  1:1164414/1‑150 (MQ=255)
                  tGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGc                                                                                                                               <  1:50752/150‑1 (MQ=255)
                  tGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTAGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGc                                                                                                                               <  1:1450266/150‑1 (MQ=255)
                         tgtCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCt                                                                                                                        >  2:1435646/1‑150 (MQ=255)
                          gtCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTg                                                                                                                       >  2:2384006/1‑150 (MQ=255)
                          gtCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTg                                                                                                                       <  2:1335132/150‑1 (MQ=255)
                           tCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTgg                                                                                                                      >  2:1924769/1‑150 (MQ=255)
                             gATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTa                                                                                                                    <  1:86347/150‑1 (MQ=255)
                             gATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTa                                                                                                                    >  1:164673/1‑150 (MQ=255)
                             gATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTAATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTa                                                                                                                    >  2:1682725/1‑150 (MQ=255)
                              aTGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAg                                                                                                                   >  2:1127742/1‑150 (MQ=255)
                                gAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGaa                                                                                                                 <  1:2384006/150‑1 (MQ=255)
                                 aaaCGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAAc                                                                                                                >  1:1455194/1‑150 (MQ=255)
                                   aCGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCACgcggcgcaagc                                                                                                                                      >  2:2321644/1‑121 (MQ=255)
                                   aCGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCACgcggcgcaagc                                                                                                                                      <  1:2321644/126‑6 (MQ=255)
                                       tctTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCa                                                                                                          <  2:572384/150‑1 (MQ=255)
                                            cTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCctt                                                                                                     >  2:462714/1‑150 (MQ=255)
                                            cTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCctt                                                                                                     <  1:1566070/150‑1 (MQ=255)
                                             tCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCTATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCcttc                                                                                                    <  2:1257916/150‑1 (MQ=255)
                                             tCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCcttc                                                                                                    <  2:1558825/150‑1 (MQ=255)
                                             tCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCcttc                                                                                                    <  2:709035/150‑1 (MQ=255)
                                             tCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCcttc                                                                                                    <  2:2250749/150‑1 (MQ=255)
                                               ggTGCGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCcttctt                                                                                                  <  1:1602205/150‑1 (MQ=255)
                                               ggTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGTATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCcttctt                                                                                                  <  2:1775141/150‑1 (MQ=255)
                                               ggTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCcttctt                                                                                                  <  2:1455194/150‑1 (MQ=255)
                                                    gTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAg                                                                                                                                >  2:1788413/1‑115 (MQ=255)
                                                    gTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAg                                                                                                                                <  1:1788413/115‑1 (MQ=255)
                                                          cGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGa                                                                                       <  1:472852/150‑1 (MQ=255)
                                                           gACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAg                                                                                      >  2:560354/1‑150 (MQ=255)
                                                             cGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGa                                                                                    <  1:1352660/150‑1 (MQ=255)
                                                               tACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTCAGCCGCGATTGGTGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGCCCACGCCGAGGATg                                                                                  >  2:1036907/1‑150 (MQ=255)
                                                                  acGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGGTCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCg                                                                               >  2:84892/1‑150 (MQ=255)
                                                                  acGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCg                                                                               >  1:1912776/1‑150 (MQ=255)
                                                                  acGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCg                                                                               >  1:1374222/1‑150 (MQ=255)
                                                                    ggCGATCAGCTTCTCGATCACATCGGCCCTGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGCACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGAt                                                                             <  2:144946/150‑1 (MQ=255)
                                                                      cGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGAtgt                                                                           <  1:2202394/150‑1 (MQ=255)
                                                                         tCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGCTGTGGc                                                                        <  2:1912776/150‑1 (MQ=255)
                                                                         tCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGAGGCGTTGCCGCCCa                                                                                                                                 <  1:1809222/93‑1 (MQ=255)
                                                                         tCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGAGGCGTTGCCGCCCa                                                                                                                                 >  2:1809222/1‑93 (MQ=255)
                                                                                 tcGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACgtgg                                                                <  2:761095/150‑1 (MQ=255)
                                                                                 tcGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACgtgg                                                                >  2:2343785/1‑150 (MQ=255)
                                                                                   gATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACgtggtg                                                              <  1:619998/150‑1 (MQ=255)
                                                                                     tCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGa                                                                                       <  2:309703/123‑1 (MQ=255)
                                                                                     tCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGa                                                                                       >  1:309703/1‑123 (MQ=255)
                                                                                     tCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACgtggtggt                                                            >  1:221606/1‑150 (MQ=255)
                                                                                      cacaTCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACgtggtggtg                                                           <  1:462714/150‑1 (MQ=255)
                                                                                      cacaTCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACgtggtggtg                                                           >  2:629808/1‑150 (MQ=255)
                                                                                       acaTCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGtggtggtgg                                                          <  2:221606/150‑1 (MQ=255)
                                                                                       acaTCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGtggtggtgg                                                          <  2:164673/150‑1 (MQ=255)
                                                                                        caTCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGacac                                                                                                                          <  1:2119646/85‑1 (MQ=255)
                                                                                        caTCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGacac                                                                                                                          >  2:2119646/1‑85 (MQ=255)
                                                                                        caTCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTgg                                                                                                                      <  2:590975/89‑1 (MQ=255)
                                                                                        caTCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTgg                                                                                                                      >  1:590975/1‑89 (MQ=255)
                                                                                        caTCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGc                                                         >  2:670525/1‑150 (MQ=255)
                                                                                         aTCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGcc                                                        >  1:2134405/1‑150 (MQ=255)
                                                                                          tCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGCACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCg                                                       >  1:178205/1‑150 (MQ=255)
                                                                                                  gTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGtt                                               <  1:2343785/150‑1 (MQ=255)
                                                                                                   tCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGTGACACCTGGTAGAACTCCTCATCATTCTTCTCCACGCCGATGATTCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTc                                              <  2:867288/150‑1 (MQ=255)
                                                                                                    cATCCTGCGCGAAGGACGGGCCGAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCa                                             >  2:1792282/1‑150 (MQ=255)
                                                                                                       ccTGCGCGAAGGACGGGCCGAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGc                                          <  1:1792282/150‑1 (MQ=255)
                                                                                                       ccTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGc                                          <  1:1924769/150‑1 (MQ=255)
                                                                                                        cTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGCa                                         >  2:96141/1‑150 (MQ=255)
                                                                                                              gAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGCAGCCGga                                   <  2:100403/150‑1 (MQ=255)
                                                                                                                aGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGCAGCCGgaga                                 <  1:629808/150‑1 (MQ=255)
                                                                                                                     gggCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGCAGCCGGAGATGTTc                            <  2:1374222/150‑1 (MQ=255)
                                                                                                                       gCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGCAGCCGGAGATGTTCAg                          >  1:1013891/1‑150 (MQ=255)
                                                                                                                       gCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGCAGCCGGAGATGTTCAg                          <  1:1585233/150‑1 (MQ=255)
                                                                                                                         cAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCAt                                            >  2:2028911/1‑130 (MQ=255)
                                                                                                                         cAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCAt                                            <  1:2028911/130‑1 (MQ=255)
                                                                                                                          aaGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGCAGCCGGAGATGTTCAGGTc                       >  1:2336373/1‑150 (MQ=255)
                                                                                                                            ggATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGCAGCCGGAGATGTTCAGGTCGa                     <  1:1435646/150‑1 (MQ=255)
                                                                                                                            ggATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGCAGCCGGAGATGTTCAGGTCGa                     <  1:1896951/150‑1 (MQ=255)
                                                                                                                                   gCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGCAGCCGGAGATGTTCAGGTCGATTTCGcc              >  1:1064393/1‑150 (MQ=255)
                                                                                                                                   gCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGCAGCCGGAGATGTTCAGGTCGATTTCGcc              <  2:650357/150‑1 (MQ=255)
                                                                                                                                     ccAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGccc                                                                      >  1:1712437/1‑92 (MQ=255)
                                                                                                                                     ccAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGccc                                                                      <  2:1712437/92‑1 (MQ=255)
                                                                                                                                     ccAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATg                                           >  2:54496/1‑119 (MQ=255)
                                                                                                                                     ccAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATg                                           <  1:54496/119‑1 (MQ=255)
                                                                                                                                     ccAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGCAGCCGGAGATGTTCAGGTCGATTTCGCCGa            >  2:2414501/1‑150 (MQ=255)
                                                                                                                                     ccAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCAAGCCGAGGATGCCGACGTGGCTGACGTGGTGGGGGCCGCAGGCGTTCATGCAGCCGGAGATGTTCAGGTCGATTTCGCCGa            >  1:2114725/1‑150 (MQ=255)
                                                                                                                                      cAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGCAGCCGGAGATGTTCAGGTCGATTTCGCCGAt           >  2:214397/1‑150 (MQ=255)
                                                                                                                                        ggCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGCAGCCGGAGATGTTCAGGTCGATTTCGCCGATGt         <  2:101862/150‑1 (MQ=255)
                                                                                                                                          cTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATg                                                                                  >  1:1258652/1‑75 (MQ=255)
                                                                                                                                          cTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATg                                                                                  <  2:1258652/75‑1 (MQ=255)
                                                                                                                                              cgccgcGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGCAGCCGGAGATGTTCAGGTCGATTTCGCCGATGTCGAACa   <  1:962913/150‑1 (MQ=255)
                                                                                                                                              cgccgcGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGCAGCCGGAGATGTTCAGGTCGATTTCGCCGATGTCGAACa   <  2:1013891/150‑1 (MQ=255)
                                                                                                                                               gccgcGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGCAGCCGGAGATGTTCAGGTCGATTTCGCCGATGTCGAACAg  <  1:962184/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GGGTCGATGCCGATACGCTGGTAGGTGTCGATGAAACGCTCTTCCTCGGTACGTTGCTCGACGTACACGGCGATCAGCTTCTCGATCACATCGGCCATGTCATCCTGCGCGAAGGACGGGCCAAGGATCTTGCCCAGGCTTGCGCCGCGTGCGGCGTTGCCGCCCAGCGACACCTGGTAGAACTCCTCACCCTTCTTGTCCACGCCGAGGATGCCGATGTGGCCCACGTGGTGGTGGCCGCAGGCGTTCATGCAGCCGGAGATGTTCAGGTCGATTTCGCCGATGTCGAACAG  >  NC_002947/2707633‑2707925

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: