Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,513,781 A→T 9.8% M109K (ATG→AAG)  PP_2203 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,513,7810AT 9.8% 867.8 / 44.5 324M109K (ATG→AAG) PP_2203hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (187/93);  new base T (12/24);  total (200/124)
Fisher's exact test for biased strand distribution p-value = 1.70e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 326 total aligned reads displayed.

GGCGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTTGACG  >  NC_002947/2513634‑2513930
                                                                                                                                                   |                                                                                                                                                      
 gcggcACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGc                                                                                                                                                     <  2:943991/150‑1 (MQ=255)
   ggcACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGccc                                                                                                                                                   >  2:4918668/1‑150 (MQ=255)
   ggcACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGccc                                                                                                                                                   <  1:3574878/150‑1 (MQ=255)
     cACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCgg                                                                                                                                                 >  1:4388011/1‑150 (MQ=255)
      aCCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGc                                                                                                                                               <  2:2358940/150‑1 (MQ=255)
      aCCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGc                                                                                                                                               <  1:719367/150‑1 (MQ=255)
      aCCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGc                                                                                                                                               >  1:149320/1‑150 (MQ=255)
         aaaGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCATGCCCATGGGCTCGCCCGGCatg                                                                                                                                            <  2:4391304/150‑3 (MQ=255)
           aGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGc                                                                                                                                          <  2:6569613/150‑1 (MQ=255)
              cgcTGAGAAACCGCGCCCCTAGGCTGAGCATCCGTAGTGGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTTGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCTTGGGCATGCCCGGCACCGCCAg                                                                                                                                       <  2:6047904/150‑1 (MQ=255)
                 tGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGcc                                                                                                                                    >  2:567080/1‑150 (MQ=255)
                 tGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGcc                                                                                                                                    >  2:2732431/1‑150 (MQ=255)
                    gaAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCtt                                                                                                                                 >  2:5235673/1‑150 (MQ=255)
                     aaaCAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGTTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACAGCCAGGCCCTTg                                                                                                                                >  1:83324/1‑150 (MQ=255)
                     aaaCAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTg                                                                                                                                >  2:2792287/1‑150 (MQ=255)
                     aaaCAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTg                                                                                                                                >  1:1578388/1‑150 (MQ=255)
                       accgcgcCCCTCGGCTGAGCATCATTTGTTGGCGACCAGGGTGTCGTGGCCATCCTGGGTCACCCCGTTGACCTGTTTGGCATCCTTGTGTTCACCCATCTCCATGCCGGGGGAGCCCTTGGGCTTGCCCGGCACCGCCAGGCCCTTGAg                                                                                                                              <  2:1039190/147‑1 (MQ=255)
                         aGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGt                                                                                                                            >  1:2240169/1‑150 (MQ=255)
                         aGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGt                                                                                                                            >  2:3815862/1‑150 (MQ=255)
                          gcgcCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTc                                                                                                                           >  1:5359270/1‑150 (MQ=255)
                           cgcCCCTCGGCTGAGCATCAGTAGTGGGCGACCACGGGGTCCTGGCCAGCCTGGGTCACCCCTATGCCCTGTTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGCGCCCATGGGCAGGCCCGGCACCGCCAGGCCCTTGAGGTCg                                                                                                                          <  2:7152520/150‑1 (MQ=255)
                              cccTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTg                                                                                                                       >  2:7579/1‑150 (MQ=255)
                              cccTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTGGAGGTCGTTg                                                                                                                       >  2:373738/1‑150 (MQ=255)
                                      tGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGc                                                                                                               >  1:3856468/1‑150 (MQ=255)
                                      tGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGc                                                                                                               >  1:2879751/1‑150 (MQ=255)
                                        cgCATAAGTTGTGGGCGACCCGGGGGTCCTGGCCAGCCGGGGCCACCCCGATGACCTGTTAGGCATCGTTGTGGCCACCCATCCCCAGGCCGGGCGAGCCCTGGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCca                                                                                                             <  1:2371355/149‑1 (MQ=255)
                                          cATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCcagc                                                                                                           >  1:4965451/1‑150 (MQ=255)
                                           aaccgtagtGGGCCCCCCCGGGGCCCTGGCCCGCCGGTGACCCCCCGATGCCCGGGTAGGCACCCTTGGGGACACCCATCTCCCTGCGGGGCGCCCACTTGGGCTTGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCcagca                                                                                                          <  1:508313/146‑1 (MQ=255)
                                            tCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCATGCCCATGGGCTCGCCCg                                                                                                                                                  <  2:4302032/110‑1 (MQ=255)
                                             cAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCggg                                                                                                        >  1:5810575/1‑150 (MQ=255)
                                                tagtTGGCGCCCACGGGGGCCTGGCCAGCCTGGGCCACCCCGAGGCCCGGGTCGGCAGCCTGGGGGTCACCCATCTCCCGGCGGGGCGCGCCCCGGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGgtgc                                                                                                     <  2:4141770/150‑1 (MQ=255)
                                                tagtTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGgcgc                                                                                                     >  1:7414881/1‑150 (MQ=255)
                                                  gtTGGCGACCAGGTTGTCCTGGCCTTCGTGGGTCACCCCGTTGACCTGTTAGGAAGCCTTGGGGACACCCATCTCCATGCCGGGCGAGCCCTGGGGCAGGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCAc                                                                                                   <  1:6049986/150‑1 (MQ=255)
                                                  gtTGGCGACCAGGGTGTCCTGGCCTTCGTGGGTCACCCCGATGACCTGGTTGGCATCCTTGTGGTCACCAATCGCCATGCCGGGGGAGCCCTTGGGCTTGCCGGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCAc                                                                                                   <  2:4041268/150‑1 (MQ=255)
                                                        gACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTc                                                                                             >  2:116998/1‑150 (MQ=255)
                                                            aCGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTcggcc                                                                                         <  2:5995331/150‑1 (MQ=255)
                                                            aCGGTGCCGTGGCCAGCCGGGGGCCCCCCCATGCCCTGGTAGGCAGCCTGGGGGTCACCCAACCCCATGCGGGGCGCCCCATTGGGCTTGCCGGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTcggcc                                                                                         <  1:2842957/150‑1 (MQ=255)
                                                              ggTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTcggccgg                                                                                       >  2:3445823/1‑150 (MQ=255)
                                                               gtgtCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTcggccggc                                                                                      <  2:4965451/150‑1 (MQ=255)
                                                                  tCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACg                                                                                   >  2:6844655/1‑150 (MQ=255)
                                                                     tGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTgg                                                                                >  2:1689846/1‑150 (MQ=255)
                                                                     tGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTgg                                                                                <  2:4857797/150‑1 (MQ=255)
                                                                     tGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTgg                                                                                >  2:213313/1‑150 (MQ=255)
                                                                     tGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTgg                                                                                >  2:7043391/1‑150 (MQ=255)
                                                                      ggCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGc                                                                               >  1:4304833/1‑150 (MQ=255)
                                                                         cATCCTGGGTCACCCCGATGACCTGGTCGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGGTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCtt                                                                            >  1:3266096/1‑150 (MQ=255)
                                                                          aTCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTc                                                                           >  1:567928/1‑150 (MQ=255)
                                                                           tCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTcg                                                                          >  1:4126959/1‑150 (MQ=255)
                                                                            ccTGGGTCACCCCGATGACCTGGTTGGCATCCTTGTGGTCACCCATGTCCTTGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTcga                                                                         <  2:1628694/150‑1 (MQ=255)
                                                                            ccTGGGCCACCCCGATGACGTGGTTGGCATCCTTGGGGACACCCATCTCCAGGCCGGGCGAGCCACTGGGCTTGCCCGCCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTcga                                                                         <  1:1141127/150‑1 (MQ=255)
                                                                                ggTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTcgacgac                                                                     >  1:2500995/1‑149 (MQ=255)
                                                                                     ccccGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTgc                                                                >  1:6406754/1‑150 (MQ=255)
                                                                                      cccGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGTATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTgcc                                                               >  2:5468231/1‑150 (MQ=255)
                                                                                      cccGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTgcc                                                               >  1:6973919/1‑150 (MQ=255)
                                                                                      cccGATGACCTGGTAGGCATCCGTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCCCGTGGCCTTCGACGAACTTgcc                                                               >  2:6314416/1‑150 (MQ=255)
                                                                                           ggccgtGGTCGGCATCCTTGTGGTCCCCCAACTCCATGCCGGGCGAGCCCTGGGGCTGGCCCGCCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGa                                                          <  1:4099133/145‑1 (MQ=255)
                                                                                              ccTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCGCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCa                                                       >  2:5847290/1‑150 (MQ=255)
                                                                                                tGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCAcc                                                     >  2:3634563/1‑150 (MQ=255)
                                                                                                   ttggCTTCTTTTTGGACACCAATCTCCTGGCCGGGCGACCCCTTGGGCTGCCCGGCCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCg                                                  <  1:5200429/148‑1 (MQ=255)
                                                                                                     ggaaTCCTTGTGGTCACCCATCTCCTTGCCGGGCGAGCCCAGGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGgt                                                <  1:2351213/147‑1 (MQ=255)
                                                                                                            ttGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCac                                         >  2:2076438/1‑150 (MQ=255)
                                                                                                            ttGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCac                                         >  1:6236297/1‑150 (MQ=255)
                                                                                                                ggTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCacgacg                                     >  2:6468470/1‑150 (MQ=255)
                                                                                                                 gTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGc                                    >  2:1308575/1‑150 (MQ=255)
                                                                                                                      ccATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGc                               >  1:6898454/1‑150 (MQ=255)
                                                                                                                      ccATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGc                               >  2:3055944/1‑150 (MQ=255)
                                                                                                                      ccATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGc                               >  1:2124610/1‑150 (MQ=255)
                                                                                                                      ccATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGc                               >  1:4677550/1‑150 (MQ=255)
                                                                                                                      ccATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTCCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGc                               >  2:6877024/1‑150 (MQ=255)
                                                                                                                        aTCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTGGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGCTCACCCCGGTGTGGCCCGACGCAAGGCGt                             >  2:4745185/1‑150 (MQ=255)
                                                                                                                        aTCTCCAGGCGGGGCGCGCCCTGGGGCAGGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGt                             <  2:703308/150‑1 (MQ=255)
                                                                                                                          ctcCATGCCGGGCGAGCCCATGGGCAGGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTgg                           <  2:3847769/150‑1 (MQ=255)
                                                                                                                          ctcCAGGCGGGGCTAGCCCAGGGCCAGCCCCGCCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTgg                           <  1:689945/150‑1 (MQ=255)
                                                                                                                             ggtGCCGGGCATGCCCATGGGCTCGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGgcgc                        >  1:2187744/3‑150 (MQ=255)
                                                                                                                             cctGCCGGGCGAGCCCAGGGGCAGGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGgcgc                        <  2:6406754/148‑1 (MQ=255)
                                                                                                                               tGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCa                      >  1:1011624/1‑150 (MQ=255)
                                                                                                                               tGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCa                      >  1:4739693/1‑150 (MQ=255)
                                                                                                                                 ccGGGCATGCCCATGGGCTCGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGGTTGGCCAGCAGGCGCACATGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGAAGCAAGGCGTGGCGCCACg                    >  1:7275249/1‑150 (MQ=255)
                                                                                                                                   gggcatgCCCATGGGCTCGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGcc                  >  1:6544756/7‑150 (MQ=255)
                                                                                                                                   gggCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGcc                  >  2:3980394/1‑150 (MQ=255)
                                                                                                                                    ggcatgCCCATGGGCTCGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCg                 >  2:247404/6‑150 (MQ=255)
                                                                                                                                    ggCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGcct                 <  2:3581492/150‑2 (MQ=255)
                                                                                                                                    ggCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCg                 >  1:1738122/1‑150 (MQ=255)
                                                                                                                                    ggCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCg                 >  2:4625606/1‑150 (MQ=255)
                                                                                                                                     gCGACCCCATGGCCAGGCCCGCCACCGCCAGGCCCTTGAGGTCGTTGCGATTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGa                <  1:499241/150‑1 (MQ=255)
                                                                                                                                      cGAGCCCAGGGGCCTGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAg               <  2:3289165/150‑1 (MQ=255)
                                                                                                                                      cGAGCCCAGGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGCGTGGCGCCACGCCTTGCGTCGTGCCACACCggg               <  2:4783201/150‑3 (MQ=38)
                                                                                                                                       gAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAgg              >  1:3826572/1‑150 (MQ=255)
                                                                                                                                        cgCCCAGGGGCAGGCCCGCCACCGCCCGGCCCTGGCGGTCGTGGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGc             <  1:4918904/149‑1 (MQ=255)
                                                                                                                                         gCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGcctt                <  1:3704393/146‑3 (MQ=255)
                                                                                                                                         gCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCg            >  2:877230/1‑150 (MQ=255)
                                                                                                                                           ccATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGtt          >  1:3058210/1‑150 (MQ=255)
                                                                                                                                           ccATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGtt          >  2:5641551/1‑150 (MQ=255)
                                                                                                                                           ccATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGtt          >  1:5463475/1‑150 (MQ=255)
                                                                                                                                           ccATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGCGTGGCGCCACGCCTTGCGTCGTGCCACACCg                 <  2:1525886/143‑1 (MQ=38)
                                                                                                                                            cATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTg         <  2:2945370/150‑1 (MQ=255)
                                                                                                                                              ggggCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGTGTGGCGCCACGCCGAGGCgccacg       <  1:4867330/149‑6 (MQ=255)
                                                                                                                                               gggCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCacgacg                                     <  2:1870651/119‑1 (MQ=255)
                                                                                                                                               gggCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCacgacg                                     >  1:1870651/1‑119 (MQ=255)
                                                                                                                                                ggCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGGGTGGCACGACGCGAGGCGTGGCGCCACGCCGAGGCGTTGCTTg     >  2:720932/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GGCGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTTGACG  >  NC_002947/2513634‑2513930

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