Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,587,415 A→C 5.1% K259Q (AAG→CAG)  PP_2263 → sugar ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,587,4150AC 5.1% 1350.1 / 41.8 393K259Q (AAG→CAG) PP_2263sugar ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (179/193);  new base C (8/12);  total (187/206)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

Only 100 of 396 total aligned reads displayed.

ACCTCGGTGAATGCCAAGCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGCTTGAACTGCGTCGCCCCTGCGTAGAGCGACGG  >  NC_002947/2587268‑2587564
                                                                                                                                                   |                                                                                                                                                     
 ccTCGGTGAATGCCAAGCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAgg                                                                                                                                                    <  2:3688175/150‑1 (MQ=255)
   tCGGTGAATGCCAAGCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGc                                                                                                                                                  >  2:4279117/1‑150 (MQ=255)
      gTGAATGCCAAGCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTc                                                                                                                                               <  1:3187824/150‑1 (MQ=255)
      gTGAATGCCAAGCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTc                                                                                                                                               <  2:1877101/150‑1 (MQ=255)
           tGCCAAGCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCt                                                                                                                                          >  2:2682280/1‑150 (MQ=255)
            gCCAAGCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTg                                                                                                                                         <  2:6871284/150‑1 (MQ=255)
            gCCAAGCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTg                                                                                                                                         <  2:7847195/150‑1 (MQ=255)
             ccAAGCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTTGCGGCGGCGGTGGTGCTTACCATCTTGCCGGGCATTCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCATGGGCTTCGCCCTgg                                                                                                                                        <  2:4807405/150‑1 (MQ=255)
             ccAAGCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTgg                                                                                                                                        >  1:1245641/1‑150 (MQ=255)
              cAAGCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTggg                                                                                                                                       >  2:2193838/1‑150 (MQ=255)
               aaGCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCCAGGGCGAAGCCCTTGGc                                                                                                                                      <  1:6904550/150‑1 (MQ=255)
               aaGCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCCAGGGCGAAGCCCTTGGc                                                                                                                                      <  2:4063360/150‑1 (MQ=255)
               aaGCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGc                                                                                                                                      >  2:6387951/1‑150 (MQ=255)
                aGCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCCAGGGCGAAGCCCTTGGcc                                                                                                                                     <  2:7303397/150‑1 (MQ=255)
                aGCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGcc                                                                                                                                     >  1:8131105/1‑150 (MQ=255)
                aGCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGcc                                                                                                                                     >  2:1387982/1‑150 (MQ=255)
                 gCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGTCTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGTCTTCGCCCTGGGCCg                                                                                                                                    <  1:6685931/150‑1 (MQ=255)
                 gCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCCAGGGCGAAGCCCTTGGcca                                                                                                                                    <  2:3854970/150‑2 (MQ=255)
                 gCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCg                                                                                                                                    >  1:6786926/1‑150 (MQ=255)
                  ccGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCgg                                                                                                                                   <  2:2161887/150‑1 (MQ=255)
                  ccGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGCGGCCCAGGGCGAAGCCCTTGGccac                                                                                                                                   <  1:7024830/150‑3 (MQ=255)
                          cgcgGTGATGACCCGCACGGTGTCGGCATCGGGCATCGCCTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGGGTgagg                                                                                                                           >  1:5258697/1‑150 (MQ=255)
                          cgcgGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTgagg                                                                                                                           >  1:958057/1‑150 (MQ=255)
                           gcgGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTgagga                                                                                                                          >  1:2838068/1‑150 (MQ=255)
                            cgGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCGATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTgaggag                                                                                                                         >  2:2902869/1‑150 (MQ=255)
                            cgGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTgaggag                                                                                                                         <  1:3908522/150‑1 (MQ=255)
                              gTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAggg                                                                                                                       <  2:6341876/150‑1 (MQ=255)
                                  tgccccGCACGGTGGCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCGGGCCGGGCATGCTGGTGCTCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAgc                                                                                                                   <  1:6436091/147‑1 (MQ=255)
                                  tgaCCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAgc                                                                                                                   <  1:2949641/150‑1 (MQ=255)
                                  tgaCCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAgc                                                                                                                   <  2:6749650/150‑1 (MQ=255)
                                   gaCCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAgcg                                                                                                                  >  2:5324116/1‑150 (MQ=255)
                                      ccGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCAt                                                                                                               >  2:7380894/1‑150 (MQ=255)
                                        gCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATgg                                                                                                             >  2:8395953/1‑150 (MQ=255)
                                            ggTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTg                                                                                                         <  2:6200582/150‑1 (MQ=255)
                                            ggTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTg                                                                                                         <  1:2331520/150‑1 (MQ=255)
                                                tCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATg                                                                                                     <  2:1149749/150‑1 (MQ=255)
                                                tCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATg                                                                                                     <  1:3455720/150‑1 (MQ=255)
                                                tCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATg                                                                                                     >  2:8673204/1‑150 (MQ=255)
                                                    catcgggcatcgACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCt                                                                                                 <  1:1396252/150‑1 (MQ=255)
                                                       cgggcatcgACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGa                                                                                              <  2:5650449/150‑1 (MQ=255)
                                                        gggcatcgACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGAc                                                                                             >  2:2454554/1‑150 (MQ=255)
                                                        gggcatcgACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGAc                                                                                             >  2:8162983/1‑150 (MQ=255)
                                                             tcgACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGc                                                                                        >  1:5677687/1‑150 (MQ=255)
                                                             tcgACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGc                                                                                        <  1:1039741/150‑1 (MQ=255)
                                                                  tGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACa                                                                                   >  1:7810272/1‑150 (MQ=255)
                                                                    gggTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAgc                                                                                 <  1:5955493/150‑1 (MQ=255)
                                                                    gggTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAgc                                                                                 >  1:7477060/1‑150 (MQ=255)
                                                                       tgtgCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCt                                                                              <  2:5174335/150‑1 (MQ=255)
                                                                       tgtgCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCt                                                                              <  2:1283540/150‑1 (MQ=255)
                                                                            tGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGttct                                                                         >  2:8780923/1‑150 (MQ=255)
                                                                              gcggcggcgGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGttcttc                                                                       >  1:2487653/1‑150 (MQ=255)
                                                                              gcggcggcgGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGttcttc                                                                       <  2:7247219/150‑1 (MQ=255)
                                                                              gcggcggcgGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGttcttc                                                                       <  2:257723/150‑1 (MQ=255)
                                                                              gcggcggcgGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGttcttc                                                                       <  2:7505751/150‑1 (MQ=255)
                                                                                      gggggTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCtgt                                                               <  2:4646826/150‑1 (MQ=255)
                                                                                      gggggTGCTGACCATCGGGCCGGGCATGCGGGGGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCtgt                                                               <  2:7112554/150‑1 (MQ=255)
                                                                                        gggTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTgg                                                             <  2:4583258/150‑1 (MQ=255)
                                                                                        gggTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTgg                                                             <  2:6461136/150‑1 (MQ=255)
                                                                                        gggTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTgg                                                             <  2:279671/150‑1 (MQ=255)
                                                                                         ggTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTggg                                                            >  1:5492910/1‑150 (MQ=255)
                                                                                           tGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGAGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGcc                                                          >  1:1860693/1‑150 (MQ=255)
                                                                                            gCTGACCATCATGCCGGGCATGCTGGTGATCTGGGTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCt                                                         >  2:582173/1‑150 (MQ=255)
                                                                                             cTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTg                                                        >  1:6251763/1‑150 (MQ=255)
                                                                                              tGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGt                                                       <  2:3134657/150‑1 (MQ=255)
                                                                                              tGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGt                                                       >  2:5954361/1‑150 (MQ=255)
                                                                                               gACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGtt                                                      >  1:2739133/1‑150 (MQ=255)
                                                                                                aCCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTtg                                                     >  1:5771992/1‑150 (MQ=255)
                                                                                                  cATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTtgct                                                   >  1:1977168/1‑150 (MQ=255)
                                                                                                  cATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTtgct                                                   >  1:2428163/1‑150 (MQ=255)
                                                                                                    tCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCACTCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTtgctgc                                                 <  1:4961196/150‑1 (MQ=255)
                                                                                                    tCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTtgctgc                                                 >  1:3027951/1‑150 (MQ=255)
                                                                                                    tCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTtgctgc                                                 <  1:7121936/150‑1 (MQ=255)
                                                                                                    tCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTtgctgc                                                 <  1:4035324/150‑1 (MQ=255)
                                                                                                      ttGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGctgctg                                               <  2:4415311/150‑1 (MQ=255)
                                                                                                       tGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTgg                                              <  1:8447901/150‑1 (MQ=255)
                                                                                                        gCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTggg                                             >  2:1221521/1‑150 (MQ=255)
                                                                                                           gggCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCAt                                          <  1:8673204/150‑1 (MQ=255)
                                                                                                           gggCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCAt                                          >  2:5043356/1‑150 (MQ=255)
                                                                                                           gggCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCAt                                          >  1:4912784/1‑150 (MQ=255)
                                                                                                             gCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGa                                        <  1:4400723/150‑1 (MQ=255)
                                                                                                                  cTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGGTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTg                                   >  2:8516279/1‑150 (MQ=255)
                                                                                                                   tGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGc                                  >  2:5392724/1‑150 (MQ=255)
                                                                                                                   tGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGc                                  >  1:6604200/1‑150 (MQ=255)
                                                                                                                   tGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGc                                  <  2:1534772/150‑1 (MQ=255)
                                                                                                                    ggTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGCt                                 >  2:4755474/1‑150 (MQ=255)
                                                                                                                     gTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGCtt                                >  1:3657830/1‑150 (MQ=255)
                                                                                                                     gTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGCtt                                <  2:6606089/150‑1 (MQ=255)
                                                                                                                         tCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGCTTGAAc                            <  2:6654854/150‑1 (MQ=255)
                                                                                                                         tCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGCTTGAAc                            >  1:2570983/1‑150 (MQ=255)
                                                                                                                          cTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGCTTGAACt                           <  1:820677/150‑1 (MQ=255)
                                                                                                                             gTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGCTTGAACTGCg                        >  1:1688383/1‑150 (MQ=255)
                                                                                                                                       aaTCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGCTTGAACTGCGTCGCCCCTGc              <  1:2170524/150‑1 (MQ=255)
                                                                                                                                       aaTCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGCTTGAACTGCGTCGCCCCTGc              <  1:4279117/150‑1 (MQ=255)
                                                                                                                                       aaTCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGCTTGAACTGCGTCGCCCCTGc              >  2:3325373/1‑150 (MQ=255)
                                                                                                                                         acacccggCCCAGGGCGAAGCCCTTGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGCTTGAACTGCGTCGCCCCTGCGt            >  1:4695382/7‑150 (MQ=255)
                                                                                                                                          cATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGCTTGAACTGCGTCGCCCCTGCGTa           <  1:1014478/150‑1 (MQ=255)
                                                                                                                                           aTGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGGTGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGCTTGAACTGCGTCGCCCCTGCGTag          >  1:5312226/1‑150 (MQ=255)
                                                                                                                                                  cAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGCTTGAACTGCGTCGCCCCTGCGTAGAGCGACg   >  2:6573302/1‑150 (MQ=255)
                                                                                                                                                  cAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGCTTGAACTGCGTCGCCCCTGCGTAGAGCGACg   >  1:3792086/1‑150 (MQ=255)
                                                                                                                                                   aaGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGCTTGAACTGCGTCGCCCCTGCGTAGAGCGACgg  >  1:7588255/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACCTCGGTGAATGCCAAGCCGATTGCCGCGGTGATGACCCGCACGGTGTCGGCATCGGGCATCGACTGGGGTGTGCTGGCGGCGGCGGGGGTGCTGACCATCTTGCCGGGCATGCTGGTGATCTGGTTCGTGCGCAATCATGTGGCCAAGGGCTTCGCCCTGGGCCGGGTGTGAGGAGGGCAGCGCATGGAGTGGATGGCCTGGACCCTGCCGACAGCGCTGTTCTTCGCCGCTGTGGGCCTGTTGCTGCTGGGCATGACCCTGCTTGAACTGCGTCGCCCCTGCGTAGAGCGACGG  >  NC_002947/2587268‑2587564

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