Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,799,989 G→T 8.7% R94M (AGG→ATG) ‡ PP_3360 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,799,9890GT 8.7% 290.1 / 10.9 92R94M (AGG→ATG) ‡PP_3360membrane protein
Reads supporting (aligned to +/- strand):  ref base G (43/41);  new base T (4/4);  total (47/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

CGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTG  >  NC_002947/3799841‑3800136
                                                                                                                                                    |                                                                                                                                                   
cgcCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAgg                                                                                                                                                    <  2:1462743/150‑1 (MQ=255)
 gcCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAggg                                                                                                                                                   <  1:1357326/150‑1 (MQ=255)
  ccGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCACCTTACTGTTGTTGCTAGCCTCCTTTCCGGCCTGCTGTAACATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGt                                                                                                                                                  <  1:657421/150‑1 (MQ=255)
              ctgttCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCtgtg                                                                                                                                      <  1:2139051/146‑1 (MQ=255)
               tGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTcag                                                                                                                                               >  1:2054649/1‑140 (MQ=255)
               tGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTcag                                                                                                                                               <  2:2054649/140‑1 (MQ=255)
                ggTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGt                                                                                                                                    <  1:435398/150‑1 (MQ=255)
                 gTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGtt                                                                                                                                   >  1:1202613/1‑150 (MQ=255)
                   ccGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGAAGCGTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCg                                                                                                                                 <  1:1232316/150‑1 (MQ=255)
                         acGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTcg                                                                                                                           >  1:289252/1‑150 (MQ=255)
                           gAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCCTGGTcagcca                                                                                                                                            <  2:1967725/131‑3 (MQ=255)
                           gAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCCTGGTcagcca                                                                                                                                            >  1:1967725/1‑129 (MQ=255)
                           gAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTcgcg                                                                                                                         >  1:1651260/1‑150 (MQ=255)
                           gAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCCGCCCGAGCGCTGGCTGACCAGGGTCAGCATCCTGTGGTTCGTGGTcgcg                                                                                                                         >  2:1590061/1‑150 (MQ=255)
                              ggCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGcc                                                                                                                      <  1:2116945/150‑1 (MQ=255)
                                 ccACTGAGATATAGAATGTCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTg                                                                                                                   <  2:516916/150‑1 (MQ=255)
                                  cACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGc                                                                                                                  <  1:580206/150‑1 (MQ=255)
                                   aCTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCa                                                                                                                 >  2:1015246/1‑150 (MQ=255)
                                       aGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGAtt                                                                                                             <  1:152914/150‑1 (MQ=255)
                                         ccAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTg                                                                                                                               <  1:1466378/130‑1 (MQ=255)
                                         ccAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTg                                                                                                                               >  2:1466378/1‑130 (MQ=255)
                                          cAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCCTGGTcagcca                                                                                                                                            >  2:496605/1‑114 (MQ=255)
                                          cAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCCTGGTcagcca                                                                                                                                            <  1:496605/116‑3 (MQ=255)
                                          cAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGc                                                                                                          <  2:1688499/150‑1 (MQ=255)
                                               aaTGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCtgtg                                                                                                     <  2:2213461/150‑1 (MQ=255)
                                               aaTGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCtgtg                                                                                                     <  2:1990647/150‑1 (MQ=255)
                                                aTGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTgg                                                                                                    >  2:300573/1‑150 (MQ=255)
                                                 tGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTggg                                                                                                   >  2:2024806/1‑150 (MQ=255)
                                                     cAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCaa                                                                                               >  1:578851/1‑150 (MQ=255)
                                                     cAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCaa                                                                                               >  1:2358715/1‑150 (MQ=255)
                                                     cAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCaa                                                                                               >  1:1588191/1‑150 (MQ=255)
                                                      aaGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAAc                                                                                              >  2:376987/1‑150 (MQ=255)
                                                        ggTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGCGCGCTGGCTGACTTGGGGGAGCAGCCTGTGGTGCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACcg                                                                                            >  1:285466/1‑150 (MQ=255)
                                                        ggTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACcg                                                                                            >  1:25787/1‑150 (MQ=255)
                                                           gCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCgcgc                                                                                         >  2:1524994/1‑150 (MQ=255)
                                                            cTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCCTGGTcagcca                                                                                                                                            <  1:1580378/98‑3 (MQ=255)
                                                            cTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCCTGGTcagcca                                                                                                                                            >  2:1580378/1‑96 (MQ=255)
                                                            cTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGcc                                                                                        >  1:456757/1‑150 (MQ=255)
                                                               tCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATc                                                                                     >  1:16307/1‑150 (MQ=255)
                                                               tCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATc                                                                                     <  2:2358715/150‑1 (MQ=255)
                                                                  ggCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCggc                                                                                  >  2:2207463/1‑150 (MQ=255)
                                                                   gCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCggcc                                                                                 >  2:203663/1‑150 (MQ=255)
                                                                      ccAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCggcctgg                                                                              <  2:25787/150‑1 (MQ=255)
                                                                       cAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTcagcag                                                                                                                                            >  2:432960/1‑87 (MQ=255)
                                                                       cAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTcagcag                                                                                                                                            <  1:432960/87‑1 (MQ=255)
                                                                        aGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCggcctgggc                                                                            <  1:1472281/150‑1 (MQ=255)
                                                                          cAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTAAGCAGCGTGTGGTTCGTGGTCGCGGCCTTGCAGCTTGGCCTGTGGGCCAACCGGGCCATCggcctgggcct                                                                          <  2:1017603/150‑1 (MQ=255)
                                                                           aaCTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGc                                                                                                          <  2:536659/117‑1 (MQ=255)
                                                                           aaCTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGc                                                                                                          >  1:536659/1‑117 (MQ=255)
                                                                           aaCTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCggcctgggccgg                                                                         >  2:1116188/1‑148 (MQ=255)
                                                                            aCTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCCGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGa                                                                        <  1:1590061/150‑1 (MQ=255)
                                                                             cTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTgg                                                                                                                              >  2:617632/1‑95 (MQ=255)
                                                                             cTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTgg                                                                                                                              <  1:617632/95‑1 (MQ=255)
                                                                              tttCGGTTGTTGCTGGCCGCCGTGCTGGTGGGCCTGAGCATGCTCGACCTGCTCGCGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGTTCGCGGCCTTGCAGATTGGCCTGTGTGCCATCCGCGTCATCGGTCTGGGCCTGAGc                                                                      <  2:542345/150‑1 (MQ=255)
                                                                                tCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGACCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCg                                                                    >  2:646554/1‑150 (MQ=255)
                                                                                       ttgCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTc                                                             <  1:1718460/150‑1 (MQ=255)
                                                                                       ttgCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTc                                                             >  1:2179646/1‑150 (MQ=255)
                                                                                       ttgCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTc                                                             <  1:586287/150‑1 (MQ=255)
                                                                                       ttgCTGGCCTCCTTGCTGGTGGGCCTGAGCAGGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTc                                                             >  1:2246233/1‑150 (MQ=255)
                                                                                               ctcctTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCa                                                     <  2:2369201/150‑1 (MQ=255)
                                                                                                 cctTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCAcc                                                   <  1:885238/150‑1 (MQ=255)
                                                                                                  ctTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAgg                                                                                                                                                    >  1:444617/1‑52 (MQ=255)
                                                                                                  ctTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAgg                                                                                                                                                    <  2:444617/52‑1 (MQ=255)
                                                                                                   ttGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGt                                                 >  1:1790697/1‑150 (MQ=255)
                                                                                                      cTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCtt                                                                                                                    <  2:2097481/80‑1 (MQ=255)
                                                                                                      cTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCtt                                                                                                                    >  1:2097481/1‑80 (MQ=255)
                                                                                                         gtggGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGAc                                           >  1:340582/1‑150 (MQ=255)
                                                                                                          tggGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGc                                                                                                                  <  2:1654830/78‑1 (MQ=255)
                                                                                                          tggGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGc                                                                                                                  >  1:1654830/1‑78 (MQ=255)
                                                                                                             gCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGCTCGGCTTGTGGGCCACCCGCGCCATCGGCCTTGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGcc                                       <  1:1374828/150‑1 (MQ=255)
                                                                                                              ccTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCt                                      <  1:376987/150‑1 (MQ=255)
                                                                                                              ccTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCATACCGACGGCCt                                      >  1:353629/1‑150 (MQ=255)
                                                                                                               cTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCAACGTACCGACGGCCTg                                     >  2:673861/1‑150 (MQ=255)
                                                                                                                  aGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCGTGTGGTTCGTGGTCGCGGCCTTGGAGATTGGCCTGTTGGCCAACCGCGCCATCGGCCTGGGCCTGCGCCGTTCCTTCGCCCGCCACCGTACCGACGGCCTGAAc                                  <  2:1955176/150‑1 (MQ=255)
                                                                                                                    cATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCa                                <  2:289252/150‑1 (MQ=255)
                                                                                                                     aTGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCac                               <  2:1651260/150‑2 (MQ=255)
                                                                                                                     aTGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAg                               <  2:826490/150‑1 (MQ=255)
                                                                                                                     aTGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAg                               <  1:1524994/150‑1 (MQ=255)
                                                                                                                        cTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGt                                                                   >  2:1975741/1‑111 (MQ=255)
                                                                                                                        cTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGt                                                                   <  1:1975741/111‑1 (MQ=255)
                                                                                                                               tGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCAc                     >  1:1975214/1‑150 (MQ=255)
                                                                                                                                gCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACt                    >  1:2347205/1‑150 (MQ=255)
                                                                                                                                   cGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGc                 >  2:7507/1‑150 (MQ=255)
                                                                                                                                   cGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGc                 >  2:49408/1‑150 (MQ=255)
                                                                                                                                    gAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGc                                                                                                  <  1:1004617/68‑1 (MQ=255)
                                                                                                                                    gAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGc                                                                                                  >  2:1004617/1‑68 (MQ=255)
                                                                                                                                      gcgcTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGACAAc              <  1:673861/150‑1 (MQ=255)
                                                                                                                                           tgctgACCCTGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCaa               <  1:2150421/143‑1 (MQ=255)
                                                                                                                                           tgctgACCCTGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCaa               >  2:2150421/2‑144 (MQ=255)
                                                                                                                                                 ccAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGc                 <  1:1705643/136‑1 (MQ=255)
                                                                                                                                                 ccAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGc                 >  2:1705643/1‑136 (MQ=255)
                                                                                                                                                 ccAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGcc    <  2:1061778/149‑1 (MQ=255)
                                                                                                                                                 ccAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGcc    >  1:1061778/1‑149 (MQ=255)
                                                                                                                                                  cAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGa                                                                        <  2:426953/80‑1 (MQ=255)
                                                                                                                                                  cAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGa                                                                        >  1:426953/1‑80 (MQ=255)
                                                                                                                                                  cAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCAGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTg  >  2:1071950/1‑150 (MQ=255)
                                                                                                                                                    gggTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGt                                                                   >  1:2058553/1‑83 (MQ=255)
                                                                                                                                                    gggTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGt                                                                   <  2:2058553/83‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCCTTGCTGGTGGGCCTGAGCATGCTCGACCTGCCCGAGCGCTGGCTGACCAGGGTCAGCAGCCTGTGGTTCGTGGTCGCGGCCTTGCAGATTGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTG  >  NC_002947/3799841‑3800136

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: