Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 980,643 C→A 5.9% R28S (CGC→AGC)  hscB → DnaJ‑like molecular chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947980,6430CA 5.9% 438.7 / 10.2 135R28S (CGC→AGC) hscBDnaJ‑like molecular chaperone
Reads supporting (aligned to +/- strand):  ref base C (68/59);  new base A (4/4);  total (72/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.56e-01

Only 100 of 139 total aligned reads displayed.

TCAACAACCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCAT  >  NC_002947/980498‑980791
                                                                                                                                                   |                                                                                                                                                    
tCAACAACCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGa                                                                                                                                                  <  1:16243/150‑1 (MQ=255)
 caacaaCCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAg                                                                                                                                                 <  1:2434092/150‑1 (MQ=255)
  aacaacCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGCTTGACGCTGGCCGGGGGTACTGCAAGTCA‑TTACGCAGTGTTAGACCTCCAGCCAAGCTTGTGTCTGGATCTCGACACCCTAGCCA‑CTAGCTATCGCCAGc                                                                                                                                                >  2:515020/1‑150 (MQ=255)
    caacCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCAGCTCGCGATAGCGAgtg                                                                                                                                               <  1:2885726/150‑3 (MQ=255)
       ccccAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCG‑CTCGCGATAGCGAG‑TGGcca                                                                                                                                          <  2:164465/150‑2 (MQ=255)
          cAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCAGCTCGCGATAGCGAG‑TGGCCagc                                                                                                                                        <  2:940278/150‑1 (MQ=255)
                  gccgtgtgtgtGGCTGCGGCGAACGCTTCAACGTTTGAGGCTGGCTATCGTGACTCCTTGTCA‑TTCCGCATTGTTTGTCCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCa                                                                                                                                <  1:965974/147‑1 (MQ=255)
                   cgGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGAATGGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCAt                                                                                                                               <  1:2830587/150‑1 (MQ=255)
                    ggTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATc                                                                                                                              <  2:2949745/150‑1 (MQ=255)
                    ggTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATc                                                                                                                              <  1:423324/150‑1 (MQ=255)
                    ggTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATc                                                                                                                              <  1:1894307/150‑1 (MQ=255)
                            tgGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGa                                                                                                                                                  <  2:2023989/122‑1 (MQ=255)
                            tgGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGa                                                                                                                                                  >  1:2023989/1‑122 (MQ=255)
                             ggCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGaa                                                                                                                                     >  1:1044192/1‑134 (MQ=255)
                             ggCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGaa                                                                                                                                     <  2:1044192/134‑1 (MQ=255)
                               cTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCTTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCttt                                                                                                                   >  2:2598190/1‑150 (MQ=255)
                               cTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCttt                                                                                                                   <  2:2758896/150‑1 (MQ=255)
                               cTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCttt                                                                                                                   <  1:2797201/150‑1 (MQ=255)
                               cTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCttt                                                                                                                   <  1:515020/150‑1 (MQ=255)
                               cTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCttt                                                                                                                   <  2:115329/150‑1 (MQ=255)
                                   ggcgAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCg                                                                                                               <  1:1763631/150‑1 (MQ=255)
                                     ggAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGTCTTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGAc                                                                                                             <  1:812722/149‑1 (MQ=255)
                                      gAAAGCTTCAACGTTTGAGGCTGGCTGTGAATACTACTCGTAATTTTCGCAGCG‑ATGACCTCCAGCCAGTGTGCCTTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCTCGAAGGCCATCCTGATCGCTTTGCCGACg                                                                                                            >  1:2066184/1‑150 (MQ=255)
                                       aaaGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGc                                                                                                           >  2:1742900/1‑150 (MQ=255)
                                            ttCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCg                                                                                                      >  2:2581759/1‑150 (MQ=255)
                                            ttCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCg                                                                                                      >  1:574846/1‑150 (MQ=255)
                                             tCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGa                                                                                                     >  2:2713813/1‑150 (MQ=255)
                                             tCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TGACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGa                                                                                                     >  1:1021239/1‑150 (MQ=255)
                                               aaCGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGc                                                                                                   <  1:2079654/150‑1 (MQ=255)
                                                   tttGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCt                                                                                                                            >  1:2458551/1‑121 (MQ=255)
                                                     tGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACATCCAGCCAAGCTTCCGCCTGGATCTCGACAAGCTGGCCA‑CTCTCAATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGagc                                                                                             >  2:2651318/1‑150 (MQ=255)
                                                             gctgTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTAtt                                                                                     <  2:2913067/150‑1 (MQ=255)
                                                              ctgTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTg                                                                                    >  1:1548502/1‑150 (MQ=255)
                                                                gtgGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGaa                                                                                                                                     >  1:2801834/1‑99 (MQ=255)
                                                                 tgGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTCCCGACGCCTCCGAGCGTGa                                                                                               <  2:1185050/136‑1 (MQ=255)
                                                                 tgGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTCCCGACGCCTCCGAGCGTGa                                                                                               >  1:1185050/1‑136 (MQ=255)
                                                                  gggTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCC‑CTCGCTATCGCGAGCTGCCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCt                                                                                <  2:1668646/150‑1 (MQ=255)
                                                                  gggTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCt                                                                                <  2:1561658/150‑1 (MQ=255)
                                                                  gggTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCt                                                                                <  1:2713813/150‑1 (MQ=255)
                                                                   ggTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTg                                                                               >  2:459880/1‑150 (MQ=255)
                                                                    gTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGt                                                                                                                                   <  1:247982/97‑1 (MQ=255)
                                                                    gTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGt                                                                                                                                   >  2:247982/1‑97 (MQ=255)
                                                                      aCTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTg                                                                                                                                              >  2:2342815/1‑84 (MQ=255)
                                                                      aCTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTg                                                                                                                                              <  1:2342815/84‑1 (MQ=255)
                                                                        tCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATc                                                                                                                                                      <  2:2453999/74‑1 (MQ=255)
                                                                        tCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATc                                                                                                                                                      >  1:926466/1‑74 (MQ=255)
                                                                        tCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATc                                                                                                                                                      <  2:926466/74‑1 (MQ=255)
                                                                        tCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAg                                                                                                                                                 >  2:394964/1‑79 (MQ=255)
                                                                        tCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGaaaa                                                                          >  2:2871865/1‑150 (MQ=255)
                                                                          cTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGTTCGCGTCGTCAACGCGGATCTCGTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAAt                                                                        >  1:2664332/1‑150 (MQ=255)
                                                                          cTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAAt                                                                        >  1:1048168/1‑150 (MQ=255)
                                                                                  ttACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTGGCCAAGCTGGCCA‑GTCGGTATCGCGAGCTGACCTTCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAACTCCGCGGc                                                                 >  1:1148318/1‑150 (MQ=255)
                                                                                  ttACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGc                                                                 >  1:340673/1‑150 (MQ=255)
                                                                                  ttACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGATGGGCCGCGAAGTCCATCCTGATCGCTTTTCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGc                                                                 >  1:690279/1‑150 (MQ=255)
                                                                                    aCGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCttt                                                                                                                   <  1:304434/98‑1 (MQ=255)
                                                                                        aTTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGt                                                                                                                                   <  2:256492/78‑1 (MQ=255)
                                                                                        aTTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGt                                                                                                                                   >  1:256492/1‑78 (MQ=255)
                                                                                           gTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAacga                                                        >  1:1472213/1‑150 (MQ=255)
                                                                                            tttGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAacgac                                                       >  1:16079/1‑150 (MQ=255)
                                                                                            tttGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAacgac                                                       >  2:393541/1‑150 (MQ=255)
                                                                                              tGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGc                                                     <  2:1996307/150‑1 (MQ=255)
                                                                                              tGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGc                                                     >  1:1996307/1‑150 (MQ=255)
                                                                                                  cTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGc                                                                                                                 <  2:1464466/86‑1 (MQ=255)
                                                                                                  cTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGc                                                                                                                 >  1:1464466/1‑86 (MQ=255)
                                                                                                  cTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTAc                                                 <  1:2871865/150‑1 (MQ=255)
                                                                                                    ccAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCa                                               <  1:1898412/150‑1 (MQ=255)
                                                                                                       gCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGAc                                            >  1:1105623/1‑150 (MQ=255)
                                                                                                            gCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTgc                                       <  2:340673/150‑1 (MQ=255)
                                                                                                             cTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTa                                                                                       <  2:2935362/101‑1 (MQ=255)
                                                                                                             cTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTa                                                                                       >  1:2935362/1‑101 (MQ=255)
                                                                                                              ttCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAac                                                          <  1:62917/129‑1 (MQ=255)
                                                                                                              ttCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAac                                                          >  2:62917/1‑129 (MQ=255)
                                                                                                                 cGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAg                                                                                                                                                 >  1:2369665/1‑39 (MQ=255)
                                                                                                                  gTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTg                                 >  1:769027/1‑150 (MQ=255)
                                                                                                                     tGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTgcgc                              <  2:16079/150‑1 (MQ=255)
                                                                                                                        aTCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGc                                                                                                   <  2:1430799/78‑1 (MQ=255)
                                                                                                                        aTCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGc                                                                                                   >  1:1430799/1‑78 (MQ=255)
                                                                                                                         tctcGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCcgcg                          >  2:1419687/1‑150 (MQ=255)
                                                                                                                         tctcGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCcgcg                          >  2:629327/1‑150 (MQ=255)
                                                                                                                             gACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTcc                                                                                                       <  1:2220177/69‑1 (MQ=255)
                                                                                                                             gACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTcc                                                                                                       >  2:2220177/1‑69 (MQ=255)
                                                                                                                             gACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTAtt                                                                                     >  1:2435407/1‑87 (MQ=255)
                                                                                                                             gACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTAtt                                                                                     <  2:2435407/87‑1 (MQ=255)
                                                                                                                             gACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCcg                                                                     >  2:243035/1‑103 (MQ=255)
                                                                                                                             gACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgc                      >  2:2284976/1‑150 (MQ=255)
                                                                                                                              aCAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcg                     >  1:620025/1‑150 (MQ=255)
                                                                                                                                      gccagcTCGCGATAGCGAG‑TGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcgc                    >  2:699041/6‑143 (MQ=255)
                                                                                                                                      gccagcTCGCGATAGCGAG‑TGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcgc                    <  1:699041/138‑1 (MQ=255)
                                                                                                                                        cagcTCGCGATAGCGAG‑TGGCCCGCGAAGTACATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCCGTGCGCCGCGTCGCGCGCGCTACct           >  1:600175/4‑150 (MQ=255)
                                                                                                                                         a‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctg          >  1:591973/1‑150 (MQ=255)
                                                                                                                                           cTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCATCAAAACTATGTCAACCAGGCCGAGCGTGAGCAGCAGGTTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGAAGCTGCGCACTGCTCCGCGTCGCGCGCGCTACctgc         >  1:1302703/1‑150 (MQ=255)
                                                                                                                                             tgCTATCGCGCGCTGGTCCTCGAAATCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgctg       <  2:2859710/149‑1 (MQ=255)
                                                                                                                                             cGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgctg       <  1:459880/150‑1 (MQ=255)
                                                                                                                                              gCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTgg      >  1:231093/1‑150 (MQ=255)
                                                                                                                                               caaTCGCGAGCTGGCCCGCGAAGTCCATCGTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGc     <  2:1148318/148‑1 (MQ=255)
                                                                                                                                                tATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGt                                                                                                 <  2:2787721/57‑1 (MQ=255)
                                                                                                                                                tATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGcc    >  2:2084853/1‑150 (MQ=255)
                                                                                                                                                tATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGcc    >  2:2084714/1‑150 (MQ=255)
                                                                                                                                                  tCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCAt  >  2:2643643/1‑150 (MQ=255)
                                                                                                                                                  tCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCAt  >  2:2606869/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCAACAACCCCAACGTGCGCGGTGAGTGTGGCTGCGGCGAAAGCTTCAACGTTTGAGGCTGGCTGTGGGTACTCCTTGTCA‑TTACGCATTGTTTGACCTCCAGCCAAGCTTCCGTCTGGATCTCGACAAGCTGGCCA‑CTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCAT  >  NC_002947/980498‑980791

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