Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,680,371 G→C 5.2% A183P (GCG→CCG) ‡ hom → homoserine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,680,3710GC 5.2% 1110.4 / 35.9 325A183P (GCG→CCG) ‡homhomoserine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (155/153);  new base C (8/9);  total (163/162)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 339 total aligned reads displayed.

CATCCCGGTGATCAAGGCCATCCGCGAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCA  >  NC_002947/1680223‑1680520
                                                                                                                                                    |                                                                                                                                                     
cATCCCGGTGATCAAGGCCATCCGCGAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGgc                                                                                                                                                      >  2:1667390/1‑150 (MQ=255)
  tCCCGGTGATCAAGGCCATCCGCGAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGgcgc                                                                                                                                                    >  1:3735687/1‑150 (MQ=255)
  tCCCGGTGATCAAGGCCATCCGCGAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGgcgc                                                                                                                                                    >  2:5612498/1‑150 (MQ=255)
  tCCCGGTGATCAAGGCCATCCGCGAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGgcgc                                                                                                                                                    >  1:4811297/1‑150 (MQ=255)
   cccGGTGATCAAGGCCATCCGCGAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCt                                                                                                                                                   <  2:3030903/150‑1 (MQ=255)
    ccGGTGATCAAGGCCATCCGCGAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTg                                                                                                                                                  <  2:245453/150‑1 (MQ=255)
    ccGGTGATCAAGGCCATCCGCGAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTg                                                                                                                                                  >  2:1898407/1‑150 (MQ=255)
     cGGTGATCAAGGCCATCCGCGAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTgg                                                                                                                                                 <  2:3735687/150‑1 (MQ=255)
      ggTGATCAAGGCCATCCGCGAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTggg                                                                                                                                                >  1:207453/1‑150 (MQ=255)
      ggTGATCAAGGCCATCCGCGAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTggg                                                                                                                                                >  1:3537144/1‑150 (MQ=255)
      ggTGATCAAGGCCATCCGCGAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTggg                                                                                                                                                <  1:873716/150‑1 (MQ=255)
        tGATCAAGGCCATCCGCGAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCt                                                                                                                                              >  2:5194718/1‑150 (MQ=255)
        tGATCAAGGCCATCCGCGAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCt                                                                                                                                              >  2:3186377/1‑150 (MQ=255)
          aTCAAGGCCATCCGCGAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGCGCCTGGGCttc                                                                                                                                            <  2:2578119/150‑3 (MQ=255)
          aTCAAGGCCATCCGCGAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGCGCCTGGGCttc                                                                                                                                            <  1:5648727/150‑3 (MQ=255)
           tCAAGGCCATCCGCGCAGGCCTGTCGGCCAACCGCATCATCTGGCTGGCCTGCATGATCCAGGGCACCGTCATCTTGATCCTCACCGACATGCGTGAGTAGGGCGGTGCCTTCCCGGTCGTGCTGGCTGAAGCCCAGCGGCTGGGCttca                                                                                                                                           <  1:1171005/150‑4 (MQ=255)
                gCCATCCGCGAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGaa                                                                                                                                      <  1:5223234/150‑1 (MQ=255)
                    tCCGCGAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGccg                                                                                                                                  >  1:1956493/1‑150 (MQ=255)
                      cgcgAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGccgac                                                                                                                                <  2:868577/150‑1 (MQ=255)
                       gcgAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGccgacc                                                                                                                               >  1:6695428/1‑150 (MQ=255)
                       gcgAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGccgacc                                                                                                                               <  2:5858556/150‑1 (MQ=255)
                          aaGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGccgacccga                                                                                                                            <  1:1823692/150‑1 (MQ=255)
                               cTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTc                                                                                                                       <  1:4285251/150‑1 (MQ=255)
                                tGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCg                                                                                                                      >  1:3549089/1‑150 (MQ=255)
                                 gTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGa                                                                                                                     >  2:4038691/1‑150 (MQ=255)
                                  tCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGAc                                                                                                                    >  2:728927/1‑150 (MQ=255)
                                    ggCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGt                                                                                                                  <  1:4161300/150‑1 (MQ=255)
                                    ggCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGt                                                                                                                  <  1:4596969/150‑1 (MQ=255)
                                     gCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTg                                                                                                                 >  1:4226755/1‑150 (MQ=255)
                                     gCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTg                                                                                                                 >  1:5083487/1‑150 (MQ=255)
                                     gCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTg                                                                                                                 >  2:84447/1‑150 (MQ=255)
                                       cAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGa                                                                                                               >  2:3527938/1‑150 (MQ=255)
                                        aaCCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGaa                                                                                                              <  2:3265344/150‑1 (MQ=255)
                                        aaCCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGaa                                                                                                              <  1:1632101/150‑1 (MQ=255)
                                             cATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCAt                                                                                                         <  1:4043554/150‑1 (MQ=255)
                                             cATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTACTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCAt                                                                                                         >  1:6126896/1‑150 (MQ=255)
                                              aTCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATc                                                                                                        <  1:1104929/150‑1 (MQ=255)
                                                    tggctggcCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCt                                                                                                  >  2:669603/1‑150 (MQ=255)
                                                    tggctggcCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACACCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCt                                                                                                  >  1:747084/1‑150 (MQ=255)
                                                      gctggcCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTgc                                                                                                <  1:5981433/150‑1 (MQ=255)
                                                      gctggcCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTgc                                                                                                <  1:6310381/150‑1 (MQ=255)
                                                            cGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACaa                                                                                          >  2:2714142/1‑150 (MQ=255)
                                                            cGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACaa                                                                                          <  1:3448712/150‑1 (MQ=255)
                                                             ggCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAg                                                                                         <  1:3606483/150‑1 (MQ=255)
                                                                 tcatcaACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGa                                                                                     <  2:3700596/150‑1 (MQ=255)
                                                                    tcaACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGa                                                                                  >  2:5699421/1‑150 (MQ=255)
                                                                       aCGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTc                                                                               <  1:807759/150‑1 (MQ=255)
                                                                        cGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCt                                                                              <  1:4485587/150‑1 (MQ=255)
                                                                             ccGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGt                                                                         <  2:4066218/150‑1 (MQ=255)
                                                                              cGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTc                                                                        >  2:5745493/1‑150 (MQ=255)
                                                                              cGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTc                                                                        >  1:3918135/1‑150 (MQ=255)
                                                                               ggCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTcc                                                                       >  2:1562639/1‑150 (MQ=255)
                                                                               ggCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTcc                                                                       >  2:2325715/1‑150 (MQ=255)
                                                                                gCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCa                                                                      >  2:4693045/1‑150 (MQ=255)
                                                                                  aaCTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATc                                                                    >  2:423019/1‑150 (MQ=255)
                                                                                  aaCTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATc                                                                    <  1:2996685/150‑1 (MQ=255)
                                                                                         tCCTCACCGACAGGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCg                                                             <  1:5001612/150‑1 (MQ=255)
                                                                                         tCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCg                                                             <  1:5431919/150‑1 (MQ=255)
                                                                                           cTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGt                                                           >  2:5244202/1‑150 (MQ=255)
                                                                                           cTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGt                                                           >  2:4275393/1‑150 (MQ=255)
                                                                                            tCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTa                                                          >  1:4663772/1‑150 (MQ=255)
                                                                                            tCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTa                                                          <  2:4500169/150‑1 (MQ=255)
                                                                                            tCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTa                                                          >  1:981580/1‑150 (MQ=255)
                                                                                              aCCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATc                                                        <  1:3988359/150‑1 (MQ=255)
                                                                                                 gAAATGCGTGAGCAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCg                                                     <  2:4576626/150‑1 (MQ=255)
                                                                                                  aaaTGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGc                                                    <  1:2921622/150‑1 (MQ=255)
                                                                                                      gCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCa                                                <  2:4811297/150‑1 (MQ=255)
                                                                                                          gagaAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTc                                            <  2:2893705/150‑1 (MQ=255)
                                                                                                          gagaAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTc                                            <  2:5417230/150‑1 (MQ=255)
                                                                                                          gagaAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTc                                            <  1:5804831/150‑1 (MQ=255)
                                                                                                            gaAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGa                                          >  2:6467835/1‑150 (MQ=255)
                                                                                                             aaGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGAc                                         >  2:1267534/1‑150 (MQ=255)
                                                                                                                ggCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAg                                      <  2:5143154/150‑1 (MQ=255)
                                                                                                                   cGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGcc                                   <  1:2785034/150‑1 (MQ=255)
                                                                                                                        cTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTacac                              <  2:1135822/150‑1 (MQ=255)
                                                                                                                         ttCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACAcc                             >  2:663224/1‑150 (MQ=255)
                                                                                                                         ttCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACAcc                             >  2:2774361/1‑150 (MQ=255)
                                                                                                                            ccGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGaa                          <  1:2703876/150‑1 (MQ=255)
                                                                                                                            ccGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGaa                          <  1:1667390/150‑1 (MQ=255)
                                                                                                                             cGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAg                         <  1:4419499/150‑1 (MQ=255)
                                                                                                                              ggACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAgg                        >  1:6565460/1‑150 (MQ=255)
                                                                                                                               gACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGt                       <  2:417003/150‑1 (MQ=255)
                                                                                                                               gACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGt                       <  2:535043/150‑1 (MQ=255)
                                                                                                                               gACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGt                       <  1:5612498/150‑1 (MQ=255)
                                                                                                                                aCGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTa                      >  2:6248606/1‑150 (MQ=255)
                                                                                                                                    gctggctgAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCAc                  <  1:1282671/150‑1 (MQ=255)
                                                                                                                                      tggctgAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCAccc                >  2:1775613/1‑150 (MQ=255)
                                                                                                                                          tgAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTCTCACCCAGCt            >  2:3905204/1‑150 (MQ=255)
                                                                                                                                          tgAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCt            <  1:3102657/150‑1 (MQ=255)
                                                                                                                                          tgAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCt            >  2:6429895/1‑150 (MQ=255)
                                                                                                                                           gtaGCCCAGCGCCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTg           >  2:5061760/3‑150 (MQ=255)
                                                                                                                                           gtaGCCCAGCGCCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTg           >  2:6284909/3‑150 (MQ=255)
                                                                                                                                           gAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTg           <  1:4038691/150‑1 (MQ=255)
                                                                                                                                            aaGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGa          <  2:1981995/150‑1 (MQ=255)
                                                                                                                                              gCCTAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACa        >  2:597139/1‑150 (MQ=255)
                                                                                                                                                ccAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAAc      >  1:2964802/1‑150 (MQ=255)
                                                                                                                                                ccAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAAc      >  2:206228/1‑150 (MQ=255)
                                                                                                                                                 cAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACg     >  1:6796588/1‑150 (MQ=255)
                                                                                                                                                 cAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACg     >  1:2344965/1‑150 (MQ=255)
                                                                                                                                                   ggCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGc   <  2:6065513/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CATCCCGGTGATCAAGGCCATCCGCGAAGGCCTGTCGGCCAACCGCATCAACTGGCTGGCCGGCATCATCAACGGCACCGGCAACTTCATCCTCACCGAAATGCGTGAGAAGGGCCGTGCCTTCCCGGACGTGCTGGCTGAAGCCCAGGCGCTGGGCTACGCCGAAGCCGACCCGACCTTCGACGTGGAAGGCATCGACGCTGCGCACAAGCTGACGATTCTGGCGTCCATCGCGTTCGGTATCCCGCTGCAGTTCGACAAGGCCTACACCGAAGGTATCACCCAGCTGACAACGGCA  >  NC_002947/1680223‑1680520

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