Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,582,862 C→T 6.0% intergenic (+78/+31) PP_4057 → / ← glgB membrane protein/1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,582,8620CT 6.0% 937.5 / 29.6 297intergenic (+78/+31)PP_4057/glgBmembrane protein/1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base C (148/131);  new base T (15/3);  total (163/134)
Fisher's exact test for biased strand distribution p-value = 1.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 298 total aligned reads displayed.

CAGCAGTGATGTGGATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGAC  >  NC_002947/4582715‑4583010
                                                                                                                                                   |                                                                                                                                                     
cagccattatcgctggccctgaacctgccgccgttgggggtgttgatcatgaagcccgcctgatcggttggctgtactggccctttcgcgggtaaACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGtctc                                                                                                                                                    >  2:5407653‑M1/96‑150 (MQ=255)
   ccattatcgctggccctgaacctgccgccgttgggggtgttgatcatgaagcccgcctgatcggttggctgtactggccccttcgcgggtaaACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCcg                                                                                                                                                 >  2:6362962‑M1/93‑150 (MQ=255)
   cagTGATGTGGATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCcg                                                                                                                                                 >  2:4224350/1‑150 (MQ=255)
   cagTGATGTGGATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCcg                                                                                                                                                 >  2:2041677/1‑150 (MQ=255)
    agTGATGTGGATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCcgc                                                                                                                                                <  1:94350/150‑1 (MQ=255)
    agTGATGTGGATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACACGCTCCTAAACGTGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACCCGTCTccccc                                                                                                                                                >  1:3522039/1‑148 (MQ=255)
        aTGTGGATGACAATGACATGAAATGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGaaa                                                                                                                                            >  2:3434429/1‑150 (MQ=255)
         tgtgGATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAg                                                                                                                                           <  1:455997/150‑1 (MQ=255)
           tgGATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAggg                                                                                                                                         <  1:1873402/150‑1 (MQ=255)
            ggATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGc                                                                                                                                        >  2:4573037/1‑150 (MQ=255)
             gATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGcc                                                                                                                                       >  1:4427204/1‑150 (MQ=255)
              gggccctgaacctgccgccgttgggggtgttgatcatgaagcccgcctgatcggttggctgtacgggcccttgcgcgggagACCCGTTCCTACCCGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCa                                                                                                                                      >  2:5862947‑M1/82‑150 (MQ=255)
              aTGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTTTACCGCGAAAGGGCCa                                                                                                                                      >  1:2814843/1‑150 (MQ=255)
               tGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAg                                                                                                                                     >  1:4670776/1‑150 (MQ=255)
               tGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAg                                                                                                                                     >  2:6727621/1‑150 (MQ=255)
                 cctgaacctgccgccgttgggggtgttgatcatgaagcccgcctgatcggtaggcagggccggccccttcgcgggtaaACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTa                                                                                                                                   >  1:550150‑M1/79‑150 (MQ=255)
                 aCAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGAAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTa                                                                                                                                   <  2:3538257/150‑1 (MQ=255)
                   aaTGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACa                                                                                                                                 >  2:298444/1‑150 (MQ=255)
                   aaTGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACa                                                                                                                                 >  1:5008660/1‑150 (MQ=255)
                    aTGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAg                                                                                                                                >  2:1713092/1‑150 (MQ=255)
                         cTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAAc                                                                                                                           <  1:3143778/150‑1 (MQ=255)
                          tGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAAcc                                                                                                                          >  2:1437379/1‑150 (MQ=255)
                           gAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCg                                                                                                                         >  2:22894/1‑150 (MQ=255)
                               gggTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCa                                                                                                                     <  2:1143775/150‑1 (MQ=255)
                                 gTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAgg                                                                                                                   >  1:86229/1‑150 (MQ=255)
                                  ttGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGc                                                                                                                  >  2:5199038/1‑150 (MQ=255)
                                  ttGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGc                                                                                                                  <  2:590048/150‑1 (MQ=255)
                                     aTTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCggg                                                                                                               <  2:906313/150‑1 (MQ=255)
                                          cAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCa                                                                                                          <  1:3594393/150‑1 (MQ=255)
                                           aGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCAt                                                                                                         >  1:2611824/1‑150 (MQ=255)
                                             ccTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGa                                                                                                       <  1:4039388/150‑1 (MQ=255)
                                               ttGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATc                                                                                                     <  1:5153156/150‑1 (MQ=255)
                                                  gggTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAac                                                                                                  <  2:2477717/150‑1 (MQ=255)
                                                        ggATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACAccccc                                                                                            >  2:2422451/1‑150 (MQ=255)
                                                             cGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAAcgg                                                                                       <  2:665428/150‑1 (MQ=255)
                                                              ggtggtAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAAcggc                                                                                      >  1:2116921/1‑150 (MQ=255)
                                                              ggtggtAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAAcggc                                                                                      <  1:1617789/150‑1 (MQ=255)
                                                                tggtAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCGCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAAcggcgg                                                                                    <  1:619026/150‑1 (MQ=255)
                                                                tggtAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCAACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAAcggcgg                                                                                    >  1:1153391/1‑150 (MQ=255)
                                                                   tAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAg                                                                                 <  2:2893010/150‑1 (MQ=255)
                                                                     ggCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGt                                                                               >  2:1565029/1‑150 (MQ=255)
                                                                       cAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTc                                                                             <  2:6415396/150‑1 (MQ=255)
                                                                        aGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCAt                                                                                                         <  2:4942208/121‑1 (MQ=255)
                                                                        aGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCa                                                                            <  1:5958696/150‑1 (MQ=255)
                                                                         gggCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAg                                                                           <  2:2134032/150‑1 (MQ=255)
                                                                         gggCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAg                                                                           <  1:5726885/150‑1 (MQ=255)
                                                                          ggccggccCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAgg                                                                          <  1:2980892/150‑1 (MQ=255)
                                                                          ggccggccCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAgg                                                                          <  1:2362044/150‑1 (MQ=255)
                                                                               gccCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAg                                                                     >  1:5406806/1‑150 (MQ=255)
                                                                               gccCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAg                                                                     >  1:461777/1‑150 (MQ=255)
                                                                                 ccctttcgcgggagacccgttcctacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtTTACCCGCGAAGGGGCCGGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGc                                                                    >  2:4900737‑M2/67‑150 (MQ=255)
                                                                                  ccTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGa                                                                  <  1:1712521/150‑1 (MQ=255)
                                                                                    ttCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATa                                                                >  2:734370/1‑150 (MQ=255)
                                                                                    ggagcgGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATa                                                                >  1:5062764/4‑150 (MQ=255)
                                                                                     tCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATaa                                                               <  2:126350/150‑1 (MQ=255)
                                                                                      cgcgGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAAt                                                              <  2:548483/150‑1 (MQ=255)
                                                                                      cgcgGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAAt                                                              <  1:5564238/150‑1 (MQ=255)
                                                                                          gggagacccgttcctacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtTTACCCGCGAAGGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATCATGGc                                                           >  1:5855675‑M2/58‑150 (MQ=255)
                                                                                          gggagacccgttcctacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtTTACCCGCGAAGGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGc                                                           >  1:1946504‑M2/58‑150 (MQ=255)
                                                                                           gTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTg                                                         >  2:1657131/1‑150 (MQ=255)
                                                                                            tAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAAt                                                              >  2:5581088/1‑144 (MQ=255)
                                                                                            tAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGc                                                        <  2:2682632/150‑1 (MQ=255)
                                                                                             aaaCCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCCGCGATAATGGCTGcc                                                       >  2:4349631/1‑150 (MQ=255)
                                                                                                CCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGt                                                    >  2:6481832‑M1/1‑150 (MQ=255)
                                                                                                CCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGCTCAGGGCCAGCGATAATGGCTGCCCGt                                                    <  2:4830747‑M1/150‑1 (MQ=255)
                                                                                                 ccgttcctacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGtg                                                   >  1:3301113‑M2/51‑150 (MQ=255)
                                                                                                     tCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGtggctg                                               >  1:3886551/1‑150 (MQ=255)
                                                                                                      ccTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGtggctgg                                              <  1:4523976/150‑1 (MQ=255)
                                                                                                       ctacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtTTACCCGCGAAGGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGtggctgg                                              >  1:4817531‑M2/45‑150 (MQ=255)
                                                                                                        tacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGcc                                            >  2:5831863‑M2/44‑150 (MQ=255)
                                                                                                        tacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGcc                                            >  1:1264281‑M2/44‑150 (MQ=255)
                                                                                                         acacGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCa                                           >  1:1395199/1‑150 (MQ=255)
                                                                                                          cacGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCAc                                          >  2:6147322/1‑150 (MQ=255)
                                                                                                           gcaggcaccgcacccggcggtgacccgtgtaggagcgggtTTACCCGCGAAGGGGCCGGCCCTGCCTACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCAc                                          >  1:707632‑M2/41‑150 (MQ=255)
                                                                                                             gggTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCt                                                          >  1:3057173/1‑131 (MQ=255)
                                                                                                             gggTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCt                                                          <  2:3057220/131‑1 (MQ=255)
                                                                                                             gggTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTc                                       >  1:533896/1‑150 (MQ=255)
                                                                                                             gggTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTc                                       <  1:4142243/150‑1 (MQ=255)
                                                                                                             gggTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGACAATGGCTGCCCGTGGCTGGCCACTTc                                       <  1:802150/150‑1 (MQ=255)
                                                                                                             aggcaccgcacccggcggtgacccgcgtaggaacgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTc                                       >  2:5303887‑M2/39‑150 (MQ=255)
                                                                                                                   ccgccgGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACt                                 <  2:5028875/150‑1 (MQ=255)
                                                                                                                     gccgGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTgg                               >  1:3790257/1‑150 (MQ=255)
                                                                                                                      ccgGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAACGGCTGCCCGTGGCTGGCCACTTCGTCACTGGc                              >  1:2417108/1‑150 (MQ=255)
                                                                                                                         ggtgcggtgcCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCAc                           <  2:3700930/150‑1 (MQ=255)
                                                                                                                                 gcCTGCGTAGGAACGGGTTTACCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCga                   <  2:2814843/150‑1 (MQ=255)
                                                                                                                                    tGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAAAGGCGGCAGGTTCAGGGCCAGCGATAATGGCTCCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGAt                >  2:1327693/1‑150 (MQ=255)
                                                                                                                                     gCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGAtt               <  1:4224350/150‑1 (MQ=255)
                                                                                                                                       gTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGc             <  2:1786687/150‑1 (MQ=255)
                                                                                                                                       gTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGc             <  2:1200375/150‑1 (MQ=255)
                                                                                                                                        tAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGcc            <  2:2343937/150‑1 (MQ=255)
                                                                                                                                          ggagcgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGa          <  1:5831863‑M2/141‑1 (MQ=255)
                                                                                                                                          ggAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGa          >  1:3957486/1‑150 (MQ=255)
                                                                                                                                            aaCGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACa        <  1:5328533/150‑1 (MQ=255)
                                                                                                                                             aCGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACAt       <  1:1493090/150‑1 (MQ=255)
                                                                                                                                               gggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTc     >  2:4022548‑M2/5‑150 (MQ=255)
                                                                                                                                               gggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTc     <  1:22894‑M2/146‑1 (MQ=255)
                                                                                                                                                 gtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGa   >  2:1504565‑M2/3‑150 (MQ=255)
                                                                                                                                                 gtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGa   >  1:1060257‑M2/3‑150 (MQ=255)
                                                                                                                                                  tCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGAc  >  2:4136967‑M2/2‑150 (MQ=255)
                                                                                                                                                  tCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGAc  <  1:6507772‑M2/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CAGCAGTGATGTGGATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGAC  >  NC_002947/4582715‑4583010

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: