Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,665,340 C→A 5.6% L69M (CTG→ATG) ‡ PP_3229 → aliphatic sulfonate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,665,3400CA 5.6% 433.5 / 7.2 134L69M (CTG→ATG) ‡PP_3229aliphatic sulfonate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (71/55);  new base A (6/2);  total (77/57)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

TGCCAAGCCTTGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTG  >  NC_002947/3665191‑3665489
                                                                                                                                                     |                                                                                                                                                     
tgcCAAGCCTTGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGTGTAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGcc                                                                                                                                                       >  1:167025/1‑150 (MQ=255)
 gcCAAGCCTTGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCt                                                                                                                                                      >  2:389367/1‑150 (MQ=255)
  ccAAGCCTTGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTg                                                                                                                                                     >  1:204260/1‑150 (MQ=255)
  ccAAGCCTTGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTg                                                                                                                                                     >  1:866128/1‑150 (MQ=255)
   cAAGCCTTGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTTACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGc                                                                                                                                                    <  2:1903779/150‑1 (MQ=255)
    aaGCCTTGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGcc                                                                                                                                                   >  2:1688480/1‑150 (MQ=255)
    aaGCCTTGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGcc                                                                                                                                                   >  2:1112238/1‑150 (MQ=255)
        cTTGGCAGCGCCGAGCAGACCCTGCGCATCGTCTATATGAAACCCTCGACCCTGCTGACCCTGCTGAAGGCCAGCGGCCCTCTCGAACAGCGCCTGCATGCCGAGGGCATTCGCGTGACCTGGAACGAGTTCCCCAGCGGCAGGCCgctg                                                                                                                                               <  1:1832556/150‑1 (MQ=255)
         ttGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCgctgc                                                                                                                                              >  1:2401429/1‑150 (MQ=255)
               gcgcCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAgg                                                                                                                                        >  1:2370503/1‑150 (MQ=255)
                 gcCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGcc                                                                                                                                      >  2:1561425/1‑150 (MQ=255)
                 cccGAGCAGACCCTGCGCATCGGCTCCCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGcc                                                                                                                                      <  1:2177838/149‑1 (MQ=255)
                   cGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCt                                                                                                                                    >  1:1215144/1‑150 (MQ=255)
                         gACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCt                                                                                                                              >  2:1039404/1‑150 (MQ=255)
                           cccTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTg                                                                                                                            >  1:2116106/1‑150 (MQ=255)
                           cccTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCCGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTg                                                                                                                            <  1:835155/150‑1 (MQ=255)
                            ccTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGctgctg                                                                                                                                            >  1:902458/1‑133 (MQ=255)
                            ccTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGctgctg                                                                                                                                            <  2:902458/133‑1 (MQ=255)
                            ccTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTgg                                                                                                                           >  1:1877169/1‑150 (MQ=255)
                            ccTGCGCATCGGCTACCAGAAATCCCCGACCCTGCTGACCCTGCGTAAGGCCAGGGGCACTCTCGAACAGCGATTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTgg                                                                                                                           <  1:1342288/150‑1 (MQ=255)
                               gcgcATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCaa                                                                                                                        <  1:1645140/150‑1 (MQ=255)
                                  cATCGGCTACCAGAAATCCTCGATCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGCGGCCCTGAACCTTGGCGATGt                                                                                                                     >  2:1183026/1‑150 (MQ=255)
                                   aTCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCGGAGGGCATTCGCGTGGCCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTAGAAGCTCGGAACCTTGGCAATGTc                                                                                                                    >  1:1944675/1‑150 (MQ=255)
                                   aTCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTc                                                                                                                    >  1:580477/1‑150 (MQ=255)
                                      ggCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGAc                                                                                                                 <  1:2431769/150‑1 (MQ=255)
                                      ggCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGAc                                                                                                                 <  2:2116106/150‑1 (MQ=255)
                                       gCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGAcc                                                                                                                >  2:1163312/1‑150 (MQ=255)
                                         tACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACctc                                                                                                              >  2:2173836/1‑150 (MQ=255)
                                         tACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACct                                                                                                               <  1:1013526/149‑1 (MQ=255)
                                         tACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACct                                                                                                               >  2:1013526/1‑149 (MQ=255)
                                           ccAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTg                                                                                                                            <  2:1728/134‑1 (MQ=255)
                                           ccAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCGGAACCTTg                                                                                                                            >  1:1728/1‑134 (MQ=255)
                                            cAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCCAGGCGCTGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTcg                                                                                                           >  1:1332704/1‑150 (MQ=255)
                                             aGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTg                                                                                                                            >  1:40614/1‑132 (MQ=255)
                                             aGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTg                                                                                                                            <  2:40614/132‑1 (MQ=255)
                                                  tCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGAc                                                                                                     <  1:389367/150‑1 (MQ=255)
                                                     tCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGtt                                                                                                  >  2:1740372/1‑150 (MQ=255)
                                                     tCGACCCTGCTGACCCTGCTGAAGACCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGtt                                                                                                  >  2:1620614/1‑150 (MQ=255)
                                                       gaccctgctgaccctgctGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGc                                                                                                                                       <  1:2407092/111‑1 (MQ=255)
                                                       gaccctgctgaccctgctGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGc                                                                                                                                       >  2:2407092/1‑111 (MQ=255)
                                                        accctgctgaccctgctgAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGcc                                                                                               <  1:506673/150‑1 (MQ=255)
                                                        accctgctgaccctgctgAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGcc                                                                                               <  1:1561425/150‑1 (MQ=255)
                                                           ctgctgaccctgctgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGac                                                                                            >  2:1159755/1‑150 (MQ=255)
                                                            tgctgaccctgctgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGaca                                                                                           <  1:1134323/150‑1 (MQ=255)
                                                             gctgaccctgctgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGacac                                                                                          <  1:1056747/150‑1 (MQ=255)
                                                               tgaccctgctgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCg                                                                                        >  2:601654/1‑150 (MQ=255)
                                                                gaccctgctgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGt                                                                                       >  2:2154839/1‑150 (MQ=255)
                                                                   cctgctgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGa                                                                                                                  >  2:1924620/1‑120 (MQ=255)
                                                                   cctgctgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGa                                                                                                                  <  1:1924620/120‑1 (MQ=255)
                                                                   cctgctgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGcc                                                                                    >  1:1906116/1‑150 (MQ=255)
                                                                   cctgctgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGcc                                                                                    >  1:2146509/1‑150 (MQ=255)
                                                                    ctgctgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCg                                                                                   <  2:2370503/150‑1 (MQ=255)
                                                                      gctgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGgt                                                                                 <  1:702650/150‑1 (MQ=255)
                                                                        tgaAGGCCAGGGGCACTCTCGAACAGCACCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGgtgt                                                                               <  2:2393253/150‑1 (MQ=255)
                                                                          aaGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGAc                                                                                                                 <  1:1182884/114‑1 (MQ=255)
                                                                          aaGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGAc                                                                                                                 >  2:1182884/1‑114 (MQ=255)
                                                                           aGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCa                                                                            <  2:1310284/150‑1 (MQ=255)
                                                                           aGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCGGGAGGCCCTGAACCGTGGCAATGGCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCa                                                                            <  2:866128/150‑1 (MQ=255)
                                                                             gCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACt                                                                          <  2:1612809/150‑1 (MQ=255)
                                                                                 ggggCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAgg                                                                      <  2:2341561/150‑1 (MQ=255)
                                                                                  gggCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAggc                                                                     <  1:1163312/150‑1 (MQ=255)
                                                                                      aCTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCa                                                                 <  1:1620614/150‑1 (MQ=255)
                                                                                      aCTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCa                                                                 >  2:433880/1‑150 (MQ=255)
                                                                                       ctctCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAg                                                                <  2:167025/150‑1 (MQ=255)
                                                                                        tctcGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAgg                                                               >  2:1309896/1‑150 (MQ=255)
                                                                                         ctcGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGgc                                                              <  2:2327456/150‑1 (MQ=255)
                                                                                          tcGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGCGTTAAACTGTAGCATGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGgcg                                                             >  2:1797488/1‑150 (MQ=255)
                                                                                          tcGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAggcggc                                                                  >  1:1018644/1‑145 (MQ=255)
                                                                                          tcGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGgcg                                                             >  2:1884843/1‑150 (MQ=255)
                                                                                            gAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAgg                                                                                                                                        >  1:2402666/1‑73 (MQ=255)
                                                                                             aaCAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGAc                                                                                                                 >  1:2370415/1‑95 (MQ=255)
                                                                                              aCAGCGCCTGCAGGCCTAGGGGATACGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTCCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGACGGTGTTCACTCAGGCGGCAGGCGCCaa                                                         <  2:1215144/150‑1 (MQ=255)
                                                                                                aGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGAAGCTGCTGGAGGCACTGAACCGGGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGt                                                       >  1:290028/1‑150 (MQ=255)
                                                                                                 gcgcCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGtt                                                      >  1:706344/1‑150 (MQ=255)
                                                                                                     cTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGa                                                    <  2:265237/148‑1 (MQ=255)
                                                                                                     cTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGa                                                    >  1:265237/1‑148 (MQ=255)
                                                                                                      tGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCAGGCCgctggg                                                                                                                                             <  1:1132523/58‑3 (MQ=38)
                                                                                                         aGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTAtt                                              >  1:135191/1‑150 (MQ=255)
                                                                                                          ggCCGAGGGCATTCGCGTGACCTGGCACGAGGTCCCCAGCGGCCTGCCGCTGCTGGAGGCGCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCGCTCAGGCGGCAGGCGCCAAGTTGACCTAtgt                                             >  2:832657/1‑148 (MQ=255)
                                                                                                           gCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTc                                            <  2:2401429/150‑1 (MQ=255)
                                                                                                                 ggCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGCCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCgg                                      >  2:914723/1‑150 (MQ=255)
                                                                                                                   cATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGa                                    <  1:1884843/150‑1 (MQ=255)
                                                                                                                    aTTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGaa                                   >  2:2134124/1‑150 (MQ=255)
                                                                                                                    aTTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGaa                                   <  1:2173836/150‑1 (MQ=255)
                                                                                                                    aTTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGaa                                   <  1:2338723/150‑1 (MQ=255)
                                                                                                                     ttCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCCGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACATATTTCGCCCGGGAAg                                  >  2:1237323/1‑150 (MQ=255)
                                                                                                                        gcgTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCAc                               <  2:2155348/150‑1 (MQ=255)
                                                                                                                            gACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGAAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTc                           >  1:889160/1‑150 (MQ=255)
                                                                                                                               cTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGcc                        >  1:2133150/1‑150 (MQ=255)
                                                                                                                                 ggCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGGTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGa                      <  2:679637/150‑1 (MQ=255)
                                                                                                                                  gCACGAGTCCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGCCGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGAc                     <  2:1944675/150‑1 (MQ=255)
                                                                                                                                    aCGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCCGTGAGACCGGAGGCGCCAAGTTGACTTATTTCGCCCGGGAAGCACCTTCGCCGACCg                   >  2:1243184/1‑150 (MQ=255)
                                                                                                                                    aCGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCCCCCGGGAAGCACCTTCGCCGACCg                   >  2:489940/1‑150 (MQ=255)
                                                                                                                                        ctccagcAGCGGCAGGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCa               >  1:15287/7‑150 (MQ=255)
                                                                                                                                              cAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGAt         <  1:1309896/150‑1 (MQ=255)
                                                                                                                                               aGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGACCTGCGCAACGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCCCCTTCGCCGCCCGCCCAGTCGATc        >  1:746352/1‑150 (MQ=255)
                                                                                                                                                 cGGCCAGCCGCTGCTGGAGACCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTAATTCGCGCGGGAAGAACCTACGGCGACCGCCCAGGCGATCCt      >  2:1132935/1‑150 (MQ=255)
                                                                                                                                                     cTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGa                                                    >  1:1174924/1‑100 (MQ=255)
                                                                                                                                                     cTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTg  >  2:2063850/1‑150 (MQ=255)
                                                                                                                                                     cTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTg  >  1:2114504/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCCAAGCCTTGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTG  >  NC_002947/3665191‑3665489

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: