Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 980,649 Δ1 bp 5.5% coding (88/522 nt) hscB → DnaJ‑like molecular chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947980,6490C. 5.5% 406.1 / 17.7 127coding (88/522 nt)hscBDnaJ‑like molecular chaperone
Reads supporting (aligned to +/- strand):  ref base C (64/55);  new base . (4/3);  total (69/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

Only 100 of 128 total aligned reads displayed.

AACAACCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCG  >  NC_002947/980500‑980796
                                                                                                                                                       |                                                                                                                                                   
aacaacCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGCTTGACGCTGGCCGGGGGTACTGCAAGTCATTACGCAGTGTTAGACCTCCAGCCAAGCTTGTGTCTGGATCTCGACACCCTAGCCACTAGCTATCGCCAGc                                                                                                                                                     >  2:515020/1‑150 (MQ=255)
     ccccAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCGCTCGCGATAGCGAGTGGcca                                                                                                                                               <  2:164465/150‑2 (MQ=255)
        cAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCAGCTCGCGATAGCGAGTGGCCagc                                                                                                                                             <  2:940278/150‑1 (MQ=255)
         aaCGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCcgcg                                                                                                                                            >  1:1561658/1‑150 (MQ=255)
             tGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGt                                                                                                                                        <  1:662857/150‑1 (MQ=255)
                gccgtgtgtgtGGCTGCGGCGAACGCTTCAACGTTTGAGGCTGGCTATCGTGACTCCTTGTCATTCCGCATTGTTTGTCCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCa                                                                                                                                     <  1:965974/147‑1 (MQ=255)
                  ggTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATc                                                                                                                                   <  2:2949745/150‑1 (MQ=255)
                  ggTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATc                                                                                                                                   <  1:423324/150‑1 (MQ=255)
                    tGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCt                                                                                                                                 >  1:2913067/1‑150 (MQ=255)
                           ggCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGaa                                                                                                                                          >  1:1044192/1‑134 (MQ=255)
                           ggCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGaa                                                                                                                                          <  2:1044192/134‑1 (MQ=255)
                             cTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCTTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCttt                                                                                                                        >  2:2598190/1‑150 (MQ=255)
                             cTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCttt                                                                                                                        <  1:2797201/150‑1 (MQ=255)
                             cTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCttt                                                                                                                        <  2:115329/150‑1 (MQ=255)
                             cTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCttt                                                                                                                        <  2:2758896/150‑1 (MQ=255)
                                 ggcgAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCg                                                                                                                    <  1:1763631/150‑1 (MQ=255)
                                    gAAAGCTTCAACGTTTGAGGCTGGCTGTGAATACTACTCGTAATTTTCGCAGCGATGACCTCCAGCCAGTGTGCCTTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCTCGAAGGCCATCCTGATCGCTTTGCCGACg                                                                                                                 >  1:2066184/1‑150 (MQ=255)
                                     aaaGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGc                                                                                                                >  2:1742900/1‑150 (MQ=255)
                                          ttCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCg                                                                                                           >  1:574846/1‑150 (MQ=255)
                                          ttCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCg                                                                                                           >  2:2581759/1‑150 (MQ=255)
                                           tCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGa                                                                                                          >  2:2713813/1‑150 (MQ=255)
                                           tCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATGACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGa                                                                                                          >  1:1021239/1‑150 (MQ=255)
                                             aaCGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGc                                                                                                        <  1:2079654/150‑1 (MQ=255)
                                               cGTTTGAGGCTGGCTGTGTGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGt                                                                                                      <  2:2066184/150‑1 (MQ=255)
                                                 tttGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCt                                                                                                                                 <  2:2458551/121‑1 (MQ=255)
                                                  ttGAGGCTGGCTGTGGGTACTCCTTGTTATTGCGCATTGTGTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCTAGCTGGCCCGCGAAGTCTATCCTTATCGCTTTTCCGTCGCCTCCGAGCGTGag                                                                                                   <  2:2930494/150‑1 (MQ=255)
                                                   tGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACATCCAGCCAAGCTTCCGCCTGGATCTCGACAAGCTGGCCACTCTCAATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGagc                                                                                                  >  2:2651318/1‑150 (MQ=255)
                                                         tggctgTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGt                                                                                                                                        <  1:2222956/106‑1 (MQ=255)
                                                         tggctgTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGt                                                                                                                                        >  2:2222956/1‑106 (MQ=255)
                                                              gtgGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGaa                                                                                                                                          <  2:2801834/99‑1 (MQ=255)
                                                              gtgGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGaa                                                                                                                                          >  1:2801834/1‑99 (MQ=255)
                                                               tgGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTCCCGACGCCTCCGAGCGTGa                                                                                                    >  1:1185050/1‑136 (MQ=255)
                                                               tgGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTCCCGACGCCTCCGAGCGTGa                                                                                                    <  2:1185050/136‑1 (MQ=255)
                                                                gggTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCCCTCGCTATCGCGAGCTGCCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCt                                                                                     <  2:1668646/150‑1 (MQ=255)
                                                                gggTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCt                                                                                     <  2:1561658/150‑1 (MQ=255)
                                                                gggTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCt                                                                                     <  1:2713813/150‑1 (MQ=255)
                                                                gggTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCt                                                                                     >  2:2726630/1‑150 (MQ=255)
                                                                 ggTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTg                                                                                    >  2:459880/1‑150 (MQ=255)
                                                                  gTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGt                                                                                                                                        <  1:247982/97‑1 (MQ=255)
                                                                  gTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGt                                                                                                                                        >  2:247982/1‑97 (MQ=255)
                                                                    aCTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTg                                                                                                                                                   >  2:2342815/1‑84 (MQ=255)
                                                                    aCTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTg                                                                                                                                                   <  1:2342815/84‑1 (MQ=255)
                                                                      tCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGaaaa                                                                               >  2:2871865/1‑150 (MQ=255)
                                                                        cTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGc                                                                                                                      >  2:1435137/1‑109 (MQ=255)
                                                                        cTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGc                                                                                                                      <  1:1435137/109‑1 (MQ=255)
                                                                        cTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAAt                                                                             >  1:1048168/1‑150 (MQ=255)
                                                                           gTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCcg                                                                          >  2:258087/1‑150 (MQ=255)
                                                                                ttACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTGGCCAAGCTGGCCAGTCGGTATCGCGAGCTGACCTTCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAACTCCGCGGc                                                                      >  1:1148318/1‑150 (MQ=255)
                                                                                ttACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGc                                                                      >  1:340673/1‑150 (MQ=255)
                                                                                  aCGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCttt                                                                                                                        >  2:304434/1‑98 (MQ=255)
                                                                                  aCGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCttt                                                                                                                        <  1:304434/98‑1 (MQ=255)
                                                                                      aTTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGt                                                                                                                                        >  1:256492/1‑78 (MQ=255)
                                                                                      aTTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCttt                                                                                                                        >  1:531280/1‑94 (MQ=255)
                                                                                      aTTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCttt                                                                                                                        <  2:531280/94‑1 (MQ=255)
                                                                                       ttGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAac                                                               <  1:1742900/150‑1 (MQ=255)
                                                                                        tGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAacg                                                              >  1:739802/1‑150 (MQ=255)
                                                                                         gTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAacga                                                             >  1:1472213/1‑150 (MQ=255)
                                                                                          tttGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAacgac                                                            >  2:393541/1‑150 (MQ=255)
                                                                                          tttGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAacgac                                                            >  1:16079/1‑150 (MQ=255)
                                                                                            tGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGc                                                          >  1:1996307/1‑150 (MQ=255)
                                                                                            tGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGc                                                          <  2:1996307/150‑1 (MQ=255)
                                                                                                cTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGc                                                                                                                      <  2:1464466/86‑1 (MQ=255)
                                                                                                cTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGc                                                                                                                      >  1:1464466/1‑86 (MQ=255)
                                                                                                cTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTAc                                                      <  1:2871865/150‑1 (MQ=255)
                                                                                                  ccAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCa                                                    <  1:1898412/150‑1 (MQ=255)
                                                                                                   cAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGc                                                          <  2:1283005/143‑1 (MQ=255)
                                                                                                   cAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGc                                                          >  1:1283005/1‑143 (MQ=255)
                                                                                                          gCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTgc                                            <  2:340673/150‑1 (MQ=255)
                                                                                                           cTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTa                                                                                            <  2:2935362/101‑1 (MQ=255)
                                                                                                            ttCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAac                                                               >  2:62917/1‑129 (MQ=255)
                                                                                                            ttCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAac                                                               <  1:62917/129‑1 (MQ=255)
                                                                                                             tCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTCCCAGACGCTGCGCa                                         <  1:1920434/150‑1 (MQ=255)
                                                                                                                gTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTg                                      >  1:769027/1‑150 (MQ=255)
                                                                                                                  cTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTgcg                                    >  2:788733/1‑150 (MQ=255)
                                                                                                                       tctcGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCcgcg                               >  2:1419687/1‑150 (MQ=255)
                                                                                                                           gACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTcc                                                                                                            <  1:2220177/69‑1 (MQ=255)
                                                                                                                           gACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTcc                                                                                                            >  2:2220177/1‑69 (MQ=255)
                                                                                                                           gACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTAtt                                                                                          >  1:2435407/1‑87 (MQ=255)
                                                                                                                           gACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTAtt                                                                                          <  2:2435407/87‑1 (MQ=255)
                                                                                                                           gACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCcg                                                                          <  1:243035/103‑1 (MQ=255)
                                                                                                                           gACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgc                           >  2:2284976/1‑150 (MQ=255)
                                                                                                                            aCAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcg                          >  1:620025/1‑150 (MQ=255)
                                                                                                                                    gccagcTCGCGATAGCGAGTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcgc                         <  1:699041/138‑1 (MQ=255)
                                                                                                                                    gccagcTCGCGATAGCGAGTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcgc                         >  2:699041/6‑143 (MQ=255)
                                                                                                                                    gccagcTCGCGATAGCGAGTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTAGCGCGCGCTAc                  >  2:1638990/6‑150 (MQ=255)
                                                                                                                                    gCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTAc                  <  2:1105623/150‑1 (MQ=255)
                                                                                                                                      cagcTCGCGATAGCGAGTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACct                >  1:1068263/4‑150 (MQ=255)
                                                                                                                                      cagcTCGCGATAGCGAGTGGCCCGCGAAGTACATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCCGTGCGCCGCGTCGCGCGCGCTACct                >  1:600175/4‑150 (MQ=255)
                                                                                                                                       aCTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctg               >  1:591973/1‑150 (MQ=255)
                                                                                                                                         cTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCATCAAAACTATGTCAACCAGGCCGAGCGTGAGCAGCAGGTTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGAAGCTGCGCACTGCTCCGCGTCGCGCGCGCTACctgc              >  1:1302703/1‑150 (MQ=255)
                                                                                                                                          tCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCa                                                    <  1:944012/111‑1 (MQ=255)
                                                                                                                                           tgCTATCGCGCGCTGGTCCTCGAAATCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgctg            <  2:2859710/149‑1 (MQ=255)
                                                                                                                                           cGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgctg            <  1:459880/150‑1 (MQ=255)
                                                                                                                                             cTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGAAACCCCCGATCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACCACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTgga          >  1:2859710/1‑149 (MQ=255)
                                                                                                                                              tATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGt                                                                                                      >  1:2787721/1‑57 (MQ=255)
                                                                                                                                              tATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGcc         >  2:2084714/1‑150 (MQ=255)
                                                                                                                                                tCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCAt       >  2:2606869/1‑150 (MQ=255)
                                                                                                                                                  gcgAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCa     >  2:427164/1‑150 (MQ=255)
                                                                                                                                                     aGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCg  <  1:2598190/150‑1 (MQ=255)
                                                                                                                                                     aGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCg  <  1:2651318/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
AACAACCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCATTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCG  >  NC_002947/980500‑980796

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