Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,593,647 C→G 10.5% Q66H (CAG→CAC PP_1399 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,593,6470CG10.5% 333.5 / 4.1 128Q66H (CAG→CACPP_1399hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (63/44);  new base G (2/12);  total (66/62)
Fisher's exact test for biased strand distribution p-value = 3.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

CGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTCGGCGAAGGCGCC  >  NC_002947/1593499‑1593796
                                                                                                                                                    |                                                                                                                                                     
 ggTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTg                                                                                                                                                     <  1:809067/150‑1 (MQ=34)
 ggTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTg                                                                                                                                                     <  1:803371/150‑1 (MQ=34)
 ggTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTg                                                                                                                                                     <  2:3281906/150‑1 (MQ=34)
  gTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTgg                                                                                                                                                    <  1:2624716/150‑1 (MQ=34)
      ggTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGTGGCCACCGGTGAGGTCGAAGCCGTCTTCGGCGAAGCGGTGgccgc                                                                                                                                                >  2:3104114/1‑150 (MQ=255)
        tCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGCccgccg                                                                                                                                              >  2:2229489/1‑150 (MQ=37)
          gAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCCGCGAAGCGGTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGt                                                                                                                                            >  2:1404289/1‑150 (MQ=255)
          gAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGt                                                                                                                                            >  1:1885499/1‑150 (MQ=37)
           aagtcttCTTCAGCGAAGCTTTGGGCAGGGGTGCTGTCTATGCCGTATTAAGTGATGCGTTGGGCTTTGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCTGCGAAGCGATGGCCGCCGGtt                                                                                                                                           <  2:2434799/144‑2 (MQ=2)
                 tCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGTTGGGCACCGGTGGGGTCGAAGGCGTGTTCTGCGTAGCGTTGGCCAGCGGTGGGGTCGAAGCCGTCTTCGGGGAAGCGCTGGCCGCCGGTGCGGTCg                                                                                                                                     >  2:2735637/1‑150 (MQ=11)
                   ttCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGGTGGCCGCCGGTGCGGTCGaa                                                                                                                                   <  2:3414542/150‑1 (MQ=21)
                   ttCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGaa                                                                                                                                   >  2:3920919/1‑150 (MQ=37)
                   ttCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTATTCGGCGAAGCGGTGGCCGCCGGTGCGGTCGaa                                                                                                                                   >  1:2403924/1‑150 (MQ=21)
                    tCAGCTAAGCGTTGGCCACCGTTGCGGTCGAATCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCTAAGCCGTATTCAGCGAAGCGTTGGCCACCGGTGCGGTCTAAGCCGTATTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAg                                                                                                                                  <  1:3233966/150‑1 (MQ=11)
                      aGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGcc                                                                                                                                <  2:3207443/150‑1 (MQ=37)
                      aGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCACCGGTGCGGTCGAAGcc                                                                                                                                <  1:50693/150‑1 (MQ=17)
                        cGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGt                                                                                                                              >  1:2330873/1‑150 (MQ=37)
                        cGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGt                                                                                                                              >  1:465830/1‑150 (MQ=37)
                            gCGTTGGGCCCCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTc                                                                                                                          <  2:2392616/150‑1 (MQ=255)
                              gTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCgg                                                                                                                        >  1:1940290/1‑150 (MQ=37)
                              cttGGCCTCGGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGGTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCgg                                                                                                                        <  2:3940321/149‑1 (MQ=255)
                               ttGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTCGCCGCCGGTGCGGTCGAAGCCGTCTTCgcc                                                                                                                       >  1:4078063/1‑148 (MQ=37)
                                tGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCTTCTTCGGCGa                                                                                                                     <  1:2442231/150‑1 (MQ=255)
                                  gCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGaa                                                                                                                    <  2:1693964/150‑1 (MQ=37)
                                  gCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGaa                                                                                                                    <  2:452495/150‑1 (MQ=37)
                                    cACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGACGTCGGCGTAGc                                                                                                                  >  2:665692/1‑150 (MQ=255)
                                     aCCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCg                                                                                                                 <  2:2578623/150‑1 (MQ=37)
                                     aCCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCg                                                                                                                 >  1:3455234/1‑150 (MQ=37)
                                       cGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGtt                                                                                                               <  2:1751821/150‑1 (MQ=38)
                                          tGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGc                                                                                                            >  1:1675110/1‑150 (MQ=39)
                                              gTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGTCGCGTTGGCCAcc                                                                                                        >  2:1185548/1‑150 (MQ=255)
                                               tCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCg                                                                                                       >  1:950441/1‑150 (MQ=39)
                                               tCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCg                                                                                                       >  1:950390/1‑150 (MQ=39)
                                               tCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCg                                                                                                       >  1:1086014/1‑150 (MQ=39)
                                               tCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCg                                                                                                       >  1:1681102/1‑150 (MQ=39)
                                                 gAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCAGCGGTTAGGTCGAAGCCTACTTCGGCGTAGACCTTGCCGACGGAGCGGTCGAAGCCGTCTTCGGCGAATCCGAGGCCACCGGt                                                                                                     >  1:4468898/1‑150 (MQ=14)
                                                  aaGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTg                                                                                                    >  1:2158512/1‑150 (MQ=39)
                                                    gCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCg                                                                                                  >  1:719767/1‑150 (MQ=39)
                                                     ccGTCTTCCGCGAAGCGTTGGCCAGCGGTGCGGTCTAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCTAAGCCGTCTTCGGCGAAGCGCTGGGCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCgg                                                                                                 <  2:2343888/150‑1 (MQ=14)
                                                     ccGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCgg                                                                                                 >  1:4005204/1‑150 (MQ=39)
                                                     ccGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCgg                                                                                                 >  2:1954093/1‑150 (MQ=39)
                                                     ccGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCgg                                                                                                 >  1:3561419/1‑150 (MQ=39)
                                                     ccGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCgg                                                                                                 <  1:2447159/150‑1 (MQ=39)
                                                     ccGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCgg                                                                                                 >  1:4401426/1‑150 (MQ=39)
                                                        tCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCg                                                                                              >  1:3722813/1‑150 (MQ=39)
                                                         cTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGa                                                                                             <  2:828031/150‑1 (MQ=39)
                                                          ttCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGaa                                                                                            >  1:723766/1‑150 (MQ=39)
                                                          ttCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGaa                                                                                            <  2:315207/150‑1 (MQ=39)
                                                           tCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAg                                                                                           <  2:3061464/150‑1 (MQ=39)
                                                             aGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGcc                                                                                         >  1:3468939/1‑150 (MQ=39)
                                                             aGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGcc                                                                                         >  1:3928424/1‑150 (MQ=39)
                                                               cGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCCCCGGTGCGGTCGAAGCCGt                                                                                       >  2:3216109/1‑150 (MQ=255)
                                                               cGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGt                                                                                       >  2:123137/1‑150 (MQ=39)
                                                                 aaGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGACCACCGGTGCGGTCGAAGCCGTCt                                                                                     >  2:4272543/1‑150 (MQ=255)
                                                                      ttGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGc                                                                                >  2:2731497/1‑150 (MQ=255)
                                                                       tGGCCAACGGAGCGGTCGAAGCCGTATTCACCGAAGCGTTGGCCAACGGTGCGGTCGAAGCCGTCTTCCGCGAAGCGATGGCCCCCGGTGCGGTCGAAGCCGACTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCg                                                                               <  2:113587/150‑1 (MQ=11)
                                                                          ccACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAg                                                                            >  2:1389977/1‑150 (MQ=255)
                                                                          ccACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGGTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAg                                                                            >  2:2932514/1‑150 (MQ=255)
                                                                           cACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGc                                                                           >  1:2137183/1‑150 (MQ=255)
                                                                                gTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGc                                                                      <  1:1954093/150‑1 (MQ=255)
                                                                                   cGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGAc                                                                   >  1:4396412/1‑150 (MQ=255)
                                                                                      tCGAAGCATTCTGCAGCGAAGCCTTGGCCAGCGGTGCTGTCGAAGCCGTATTCGGCGAAGCGGTGGGCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTgg                                                                <  2:3939936/150‑1 (MQ=14)
                                                                                      ccGACGCCGTCTTCAGCGAAGCGTTGGCCCCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGGTGGCCGCCGGTGCGGTCGAAGCCGTCTTCTGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTgg                                                                <  2:1145472/149‑1 (MQ=255)
                                                                                        taaGCCGTCTTCCGCGAAGCGGTGGGCCCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGGTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGt                                                              <  2:708663/149‑1 (MQ=17)
                                                                                        gAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGt                                                              >  1:2542450/1‑150 (MQ=255)
                                                                                          atccttCTGCAGCGCAGCGTTGGGCACGGGTGCGGTCTAAGCCGTATTCGGCGAAGCGGTGGGCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGttt                                                            <  2:4517829/145‑1 (MQ=17)
                                                                                           gCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTg                                                           <  1:2115780/150‑1 (MQ=255)
                                                                                             cggcTTCACCGAAGCGTTGGCCACCGGTGCGGTCGAAGCGGTCTTCGGCGTAGGTCTGGGCGCCGGTTCGGTAGAAGCCGTCTTCGGCGAAGCGTTGGCCTCCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGaa                                                         <  1:3641434/147‑1 (MQ=17)
                                                                                              gTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGCTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGt                                                              <  1:810752/144‑1 (MQ=255)
                                                                                              gTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGCTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGt                                                              >  2:810752/1‑144 (MQ=255)
                                                                                                 tcccgCGAAGCGGCGGCCCCGGGAGCTCTCGAAGCGGTCTCCGGAGAAGCGGTGGGCGCGGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGc                                                     <  2:2226667/146‑1 (MQ=17)
                                                                                                  tCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCg                                                    >  2:3930646/1‑150 (MQ=255)
                                                                                                   cAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCgg                                                   <  1:1498415/150‑1 (MQ=255)
                                                                                                        aaGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGGTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAg                                              <  2:3551513/150‑1 (MQ=255)
                                                                                                          gCGTTGGCCACCGGTGCGGCCGAAGCCGTCTTCTGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGgt                                            <  1:3368055/150‑1 (MQ=255)
                                                                                                                 ccACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCg                                     >  2:4098436/1‑150 (MQ=255)
                                                                                                                    ccGGTGCGGTCGAAGCCGTCTTCTGCGAAGCGTTGGCCTCCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTcg                                  <  1:1338078/150‑1 (MQ=255)
                                                                                                                     cGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGGACTGCGGCAAAGGTGTGAGCGTcgg                                 >  1:4056695/1‑150 (MQ=255)
                                                                                                                            gTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTcggtccggtc                          <  2:3561419/150‑1 (MQ=255)
                                                                                                                             tCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCg                         >  1:3775941/1‑150 (MQ=255)
                                                                                                                              cGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGa                        >  2:741701/1‑150 (MQ=255)
                                                                                                                               gAAGCCGTCTTCTGCGAAGCGCTGGCCTCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGaa                       <  1:3104114/150‑1 (MQ=255)
                                                                                                                               gAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGaa                       <  1:2735637/150‑1 (MQ=255)
                                                                                                                               gAAGCCGTCTCAGGCGAAGCGGTGGGCGCCGGTGCGGTCGAAGCCGTCTTCTGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGaa                       <  2:4489370/150‑1 (MQ=25)
                                                                                                                                aaGCCGTCTTCTGCGAAGCGCTGGCCTCCGGTGCGGTCGAAGCCGTCTTCTGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAg                      <  2:1809208/150‑1 (MQ=255)
                                                                                                                                  gCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGcc                    >  2:1995604/1‑150 (MQ=255)
                                                                                                                                  gCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGcc                    >  2:1995457/1‑150 (MQ=255)
                                                                                                                                  gCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGcc                    >  1:3401075/1‑150 (MQ=255)
                                                                                                                                         taggCGATGCGGTGGGCGTCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTCg             <  2:1317157/148‑1 (MQ=255)
                                                                                                                                          cGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTCgg            <  1:2731497/150‑1 (MQ=255)
                                                                                                                                           ggCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTCGGc           <  2:2330873/150‑1 (MQ=255)
                                                                                                                                              gcaGCGGTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGTTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTCGGCGaa        <  1:3233230/148‑1 (MQ=255)
                                                                                                                                              gAAGCGTTGGCCCCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTCGGCGaa        <  1:1037504/150‑1 (MQ=255)
                                                                                                                                              gAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTCGGCGaa        >  1:3933508/1‑150 (MQ=255)
                                                                                                                                              gAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTCGGCGaa        >  2:1130212/1‑150 (MQ=255)
                                                                                                                                              gAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTCGGCGaa        <  1:4337668/150‑1 (MQ=255)
                                                                                                                                              gAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTCGGCGaa        >  2:887455/1‑150 (MQ=255)
                                                                                                                                                 gcgcTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTCGGCGAAGgc     >  1:1159802/1‑150 (MQ=255)
                                                                                                                                                 gcgaTGGCCTCCGGTGCGGTCGAAGCCGTCTTCCGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTCGGCGAAGgc     <  1:2483944/150‑1 (MQ=255)
                                                                                                                                                    cTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTCGGCGAAGGCGcc  >  1:607465/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTCGGCGAAGGCGCC  >  NC_002947/1593499‑1593796

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