Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,665,341 T→G 5.1% L69R (CTG→CGG) ‡ 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,3410TG 5.1% 391.1 / 8.8 138L69R (CTG→CGG) ‡PP_3229aliphatic sulfonate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (73/57);  new base G (5/2);  total (79/59)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

Only 100 of 142 total aligned reads displayed.

GCCAAGCCTTGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGC  >  NC_002947/3665192‑3665490
                                                                                                                                                     |                                                                                                                                                     
 ccAAGCCTTGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTg                                                                                                                                                      >  1:204260/1‑150 (MQ=255)
 ccAAGCCTTGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTg                                                                                                                                                      >  1:866128/1‑150 (MQ=255)
   aaGCCTTGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGcc                                                                                                                                                    >  2:1688480/1‑150 (MQ=255)
   aaGCCTTGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGcc                                                                                                                                                    >  2:1112238/1‑150 (MQ=255)
   aaGCCTTGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGAACTGGCACGAGTTCCCCAGCGGCCTGcc                                                                                                                                                    >  1:679637/1‑150 (MQ=255)
       cTTGGCAGCGCCGAGCAGACCCTGCGCATCGTCTATATGAAACCCTCGACCCTGCTGACCCTGCTGAAGGCCAGCGGCCCTCTCGAACAGCGCCTGCATGCCGAGGGCATTCGCGTGACCTGGAACGAGTTCCCCAGCGGCAGGCCgctg                                                                                                                                                <  1:1832556/150‑1 (MQ=255)
       cTTGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCACCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCgctg                                                                                                                                                >  1:2393253/1‑150 (MQ=255)
        ttGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCgctgc                                                                                                                                               <  1:888900/150‑1 (MQ=255)
        ttGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCgctgc                                                                                                                                               >  1:2401429/1‑150 (MQ=255)
             aGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAg                                                                                                                                          <  2:2135308/150‑1 (MQ=255)
              gcgcCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGCACCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCAACTCGATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAgg                                                                                                                                         <  1:1995004/150‑1 (MQ=255)
              gcgcCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAgg                                                                                                                                         >  2:2217556/1‑150 (MQ=255)
              gcgcCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGATGCTGGAgg                                                                                                                                         <  1:809718/150‑1 (MQ=255)
                gcCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGcc                                                                                                                                       >  2:1561425/1‑150 (MQ=255)
                  cGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCt                                                                                                                                     >  1:1215144/1‑150 (MQ=255)
                          cccTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTg                                                                                                                             <  1:1924155/150‑1 (MQ=255)
                           ccTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGctgctg                                                                                                                                             >  1:902458/1‑133 (MQ=255)
                           ccTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGctgctg                                                                                                                                             <  2:902458/133‑1 (MQ=255)
                           ccTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTgg                                                                                                                            >  2:185836/1‑150 (MQ=255)
                           ccTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTgg                                                                                                                            >  1:1877169/1‑150 (MQ=255)
                                 cATCGGCTACCAGAAATCCTCGATCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGCGGCCCTGAACCTTGGCGATGt                                                                                                                      >  2:1183026/1‑150 (MQ=255)
                                 cATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGt                                                                                                                      <  1:1203107/150‑1 (MQ=255)
                                  aTCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTc                                                                                                                     >  1:580477/1‑150 (MQ=255)
                                     ggCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGAc                                                                                                                  <  2:2116106/150‑1 (MQ=255)
                                     ggCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGAc                                                                                                                  <  1:2431769/150‑1 (MQ=255)
                                     ggCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGAc                                                                                                                  <  1:1605763/150‑1 (MQ=255)
                                      gCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGAcc                                                                                                                 >  2:1163312/1‑150 (MQ=255)
                                        tACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACctc                                                                                                               >  2:2173836/1‑150 (MQ=255)
                                        tACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACct                                                                                                                >  2:1013526/1‑149 (MQ=255)
                                        tACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACct                                                                                                                <  1:1013526/149‑1 (MQ=255)
                                          ccAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTg                                                                                                                             <  2:1728/134‑1 (MQ=255)
                                          ccAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCGGAACCTTg                                                                                                                             >  1:1728/1‑134 (MQ=255)
                                           cAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTcg                                                                                                            >  1:2409909/1‑150 (MQ=255)
                                           cAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCCAGGCGCTGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTcg                                                                                                            >  1:1332704/1‑150 (MQ=255)
                                            aGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTg                                                                                                                             >  1:40614/1‑132 (MQ=255)
                                            aGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTg                                                                                                                             <  2:40614/132‑1 (MQ=255)
                                                   cTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAAACTGGGCAATGTCGACCTCTCGGCGGACGt                                                                                                    >  1:2341561/1‑150 (MQ=255)
                                                    tCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGTCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGtt                                                                                                   >  1:1127753/1‑150 (MQ=255)
                                                    tCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGtt                                                                                                   >  2:1740372/1‑150 (MQ=255)
                                                      gaccctgctgaccctgctGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGc                                                                                                                                        <  1:2407092/111‑1 (MQ=255)
                                                          ctgctgaccctgctgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGac                                                                                             >  2:1159755/1‑150 (MQ=255)
                                                            gctgaccctgctgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGacac                                                                                           <  1:1056747/150‑1 (MQ=255)
                                                             ctgaccctgctgaAGGCCAGGGACACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACAcc                                                                                          >  2:308289/1‑150 (MQ=255)
                                                              tgaccctgctgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCg                                                                                         >  2:601654/1‑150 (MQ=255)
                                                                 ccctgctgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGc                                                                                      <  1:1688480/150‑1 (MQ=255)
                                                                  cctgctgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGa                                                                                                                   >  2:1924620/1‑120 (MQ=255)
                                                                  cctgctgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGa                                                                                                                   <  1:1924620/120‑1 (MQ=255)
                                                                  cctgctgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGcc                                                                                     >  1:1310284/1‑150 (MQ=255)
                                                                     gctgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGgt                                                                                  <  1:702650/150‑1 (MQ=255)
                                                                       tgaAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGgtgt                                                                                <  1:277962/150‑1 (MQ=255)
                                                                       tgaAGGCCAGGGGCACTCTCGAACAGCACCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGgtgt                                                                                <  2:2393253/150‑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)
                                                                              cAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCa                                                                         >  2:1294297/1‑150 (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)
                                                                                      ctctCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCCGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAg                                                                 >  2:1638804/1‑150 (MQ=255)
                                                                                       tctcGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAgg                                                                >  2:1309896/1‑150 (MQ=255)
                                                                                         tcGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAggcggc                                                                   >  1:1018644/1‑145 (MQ=255)
                                                                                         tcGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAggcggc                                                                   <  2:1018644/145‑1 (MQ=255)
                                                                                         tcGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGgcg                                                              >  2:1884843/1‑150 (MQ=255)
                                                                                           gAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAgg                                                                                                                                         <  2:2402666/73‑1 (MQ=255)
                                                                                            aaCAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGAc                                                                                                                  >  1:2370415/1‑95 (MQ=255)
                                                                                            aaCAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGAc                                                                                                                  <  2:2370415/95‑1 (MQ=255)
                                                                                             aCAGCGCCTGCAGGCCTAGGGGATACGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTCCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGACGGTGTTCACTCAGGCGGCAGGCGCCaa                                                          <  2:1215144/150‑1 (MQ=255)
                                                                                               aGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGt                                                        >  1:2155348/1‑150 (MQ=255)
                                                                                               aGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGAAGCTGCTGGAGGCACTGAACCGGGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGt                                                        >  1:290028/1‑150 (MQ=255)
                                                                                                    cTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGa                                                     <  2:265237/148‑1 (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                                    <  1:2173836/150‑1 (MQ=255)
                                                                                                                   aTTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGaa                                    <  1:2338723/150‑1 (MQ=255)
                                                                                                                    ttCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACATTTGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAg                                   >  1:427760/1‑150 (MQ=255)
                                                                                                                    ttCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCCGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACATATTTCGCCCGGGAAg                                   >  2:1237323/1‑150 (MQ=255)
                                                                                                                       gcgTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCAc                                <  2:2155348/150‑1 (MQ=255)
                                                                                                                              cTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGcc                         >  1:2133150/1‑150 (MQ=255)
                                                                                                                                 gCACGAGTCCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGCCGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGAc                      <  2:1944675/150‑1 (MQ=255)
                                                                                                                                   aCGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCg                    >  1:1149145/1‑150 (MQ=255)
                                                                                                                                   aCGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCCCCCGGGAAGCACCTTCGCCGACCg                    >  2:489940/1‑150 (MQ=255)
                                                                                                                                      aGTTCCCCAGAGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGccc                 >  2:615809/1‑150 (MQ=255)
                                                                                                                                       ctccagcAGCGGCAGGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCa                >  1:15287/7‑150 (MQ=255)
                                                                                                                                        tccagcAGCGGCAGGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGCCCTCCTGGAGGTCGTTGCCGACGCTATGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAg               >  1:2124528/6‑150 (MQ=255)
                                                                                                                                        tccagcAGCGGCAGGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAg               >  1:1505239/6‑150 (MQ=255)
                                                                                                                                          cagcAGCGGCAGGCCGCTGCTGGAGGCCCTGAACCTAGGAAATGTCGACCTCTCTGCGGGCGGTGCCGGCACCGGGCCGGGGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGc             >  2:2104129/4‑150 (MQ=255)
                                                                                                                                              aGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGACCTGCGCAACGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCCCCTTCGCCGCCCGCCCAGTCGATc         >  1:746352/1‑150 (MQ=255)
                                                                                                                                                cGGCCAGCCGCTGCTGGAGACCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTAATTCGCGCGGGAAGAACCTACGGCGACCGCCCAGGCGATCCt       >  2:1132935/1‑150 (MQ=255)
                                                                                                                                                   ccTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGt    <  2:1906116/150‑1 (MQ=255)
                                                                                                                                                    cTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGa                                                     >  1:1174924/1‑100 (MQ=255)
                                                                                                                                                    cTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTg   >  1:2114504/1‑150 (MQ=255)
                                                                                                                                                     tGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCaa                                                          >  2:1096145/1‑94 (MQ=255)
                                                                                                                                                     tGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGc  <  2:2146509/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCAAGCCTTGCCAGCGCCGAGCAGACCCTGCGCATCGGCTACCAGAAATCCTCGACCCTGCTGACCCTGCTGAAGGCCAGGGGCACTCTCGAACAGCGCCTGCAGGCCGAGGGCATTCGCGTGACCTGGCACGAGTTCCCCAGCGGCCTGCCGCTGCTGGAGGCCCTGAACCTTGGCAATGTCGACCTCTCGGCGGACGTTGCCGACACCGTGCCGGTGTTCACTCAGGCGGCAGGCGCCAAGTTGACCTATTTCGCCCGGGAAGCACCTTCGCCGACCGCCCAGGCGATCCTGGTGC  >  NC_002947/3665192‑3665490

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: ‑