Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,171,910 C→A 6.1% A150S (GCG→TCG)  glmS ← L‑glutamine/D‑fructose‑6‑phosphate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,171,9100CA 6.1% 253.4 / 3.4 80A150S (GCG→TCG) glmSL‑glutamine/D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (30/45);  new base A (2/3);  total (32/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

AGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCGATGATGC  >  NC_002947/6171765‑6172056
                                                                                                                                                 |                                                                                                                                                   
aGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCgccgcc                                                                                                                                                <  2:286833/150‑1 (MQ=255)
   gcCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCACAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTa                                                                                                                                             >  1:388198/1‑150 (MQ=255)
    ccAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTaa                                                                                                                                            <  1:1247857/150‑1 (MQ=255)
     cAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATCTGGGGTTCGGTACTACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAg                                                                                                                                           <  2:37456/150‑1 (MQ=255)
     cAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAg                                                                                                                                           <  2:1204723/150‑1 (MQ=255)
     cAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAg                                                                                                                                           <  2:848720/150‑1 (MQ=255)
     cAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAg                                                                                                                                           >  1:187578/1‑150 (MQ=255)
      aGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAgg                                                                                                                                          <  2:12149/150‑1 (MQ=255)
      aGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAgg                                                                                                                                          <  1:1273452/150‑1 (MQ=255)
           tacgcgcGCCAAGAATTTCTCTCCAAGGCCCAGGCCGATACCCAGCGGGCTCCCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTc                                                                                                                                     <  2:945009/144‑1 (MQ=255)
                  gCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAgg                                                                                                                              >  2:1250486/1‑150 (MQ=255)
                    cAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTc                                                                                                                            >  2:516891/1‑150 (MQ=255)
                    cAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTc                                                                                                                            >  1:881921/1‑150 (MQ=255)
                      ggAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTg                                                                                                                          <  1:510776/150‑1 (MQ=255)
                         aGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGaa                                                                                                                       <  1:612661/150‑1 (MQ=255)
                          gTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAAt                                                                                                                      <  2:467701/150‑1 (MQ=255)
                              tcGCCAAGGCCCAGGCCGATAACCAGCTGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGctc                                                                                                                  >  2:1196912/1‑150 (MQ=255)
                              tcGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGctc                                                                                                                  <  2:197339/150‑1 (MQ=255)
                              tcGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGctc                                                                                                                  <  2:197351/150‑1 (MQ=255)
                                  cAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCa                                                                                                              >  2:1068789/1‑150 (MQ=255)
                                  cAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCa                                                                                                              >  2:1068676/1‑150 (MQ=255)
                                   aaGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAg                                                                                                             <  1:1068676/150‑1 (MQ=255)
                                   aaGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAg                                                                                                             <  1:1068789/150‑1 (MQ=255)
                                    aGGCCCAGGCCGATAACCAGCGGGGTCCCTCGGCGTGCTGCCACCAGGTGGTCCGGCTGCTTGGCACTGATCAGCGACAGGCCGTAAGCGACGTGCTCGGCCATCAGTTCCTGTTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAgg                                                                                                            <  1:247900/150‑1 (MQ=255)
                                     ggCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAggg                                                                                                           <  2:627014/150‑1 (MQ=255)
                                        ccAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCATCAgtgagt                                                                                                        >  1:212733/1‑145 (MQ=255)
                                           ggCCGATAACCAGCGGGCTACCGCTGCGTGCTTCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTgtggt                                                                                                     <  2:225501/150‑1 (MQ=255)
                                           ggCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTACGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTgtggt                                                                                                     <  1:99747/150‑1 (MQ=255)
                                              cGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGACGCCCTTAAg                                                                                                                                           >  2:694334/1‑110 (MQ=255)
                                              cGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGACGCCCTTAAg                                                                                                                                           <  1:694334/110‑1 (MQ=255)
                                               gATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGAt                                                                                                 >  1:1290985/1‑150 (MQ=255)
                                                aTAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATc                                                                                                <  2:1290985/150‑1 (MQ=255)
                                                    ccAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCAAGGCCGTAAGCGCCGGGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGt                                                                                            <  1:50792/150‑1 (MQ=255)
                                                     cAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTg                                                                                           >  2:767879/1‑150 (MQ=255)
                                                     cAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGACGCCCTTAAg                                                                                                                                           >  1:244016/1‑103 (MQ=255)
                                                     cAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGACGCCCTTAAg                                                                                                                                           <  2:244016/103‑1 (MQ=255)
                                                       gCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGa                                                                                         <  2:388198/150‑1 (MQ=255)
                                                          ggCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGa                                                                                      <  1:1227587/150‑1 (MQ=255)
                                                          ggCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGa                                                                                      <  1:1250486/150‑1 (MQ=255)
                                                             tACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAg                                                                                              >  2:1025635/1‑139 (MQ=255)
                                                             tACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAg                                                                                              <  1:1025635/139‑1 (MQ=255)
                                                                cACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACtt                                                                                >  2:29343/1‑150 (MQ=255)
                                                                  cTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAgg                                                                                             <  1:630868/135‑1 (MQ=255)
                                                                  cTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAgg                                                                                             >  2:630868/1‑135 (MQ=255)
                                                                        gctgcCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCg                                                                        <  1:301163/150‑1 (MQ=255)
                                                                                         cGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGa                                                                >  2:1025802/1‑141 (MQ=255)
                                                                                         cGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGa                                                                <  1:1025802/141‑1 (MQ=255)
                                                                                             tgctTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAg                                                   <  1:516891/150‑1 (MQ=255)
                                                                                             ggcctGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAg                                                   >  2:756044/5‑150 (MQ=255)
                                                                                               ctTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGcc                                                 >  1:1364528/1‑150 (MQ=255)
                                                                                                ttGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGAc                                                                                        >  2:422425/1‑110 (MQ=255)
                                                                                                ttGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGAc                                                                                        <  1:422425/110‑1 (MQ=255)
                                                                                                ttGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGAGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCg                                                                              <  2:710412/120‑1 (MQ=255)
                                                                                                ttGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGAGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCg                                                                              >  1:710412/1‑120 (MQ=255)
                                                                                                   gCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAgg                                             <  2:187578/150‑1 (MQ=255)
                                                                                                    cACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGc                                            >  2:84455/1‑150 (MQ=255)
                                                                                                    cACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGc                                            >  2:33223/1‑150 (MQ=255)
                                                                                                     aCTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGcc                                           >  2:285606/1‑150 (MQ=255)
                                                                                                       tGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCt                                         >  1:1203928/1‑150 (MQ=255)
                                                                                                        gATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAg                                                                                              >  2:1191197/1‑96 (MQ=255)
                                                                                                        gATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAg                                                                                              <  1:1191197/96‑1 (MQ=255)
                                                                                                         aTCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTc                                       <  1:84455/150‑1 (MQ=255)
                                                                                                         aTCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTc                                       <  2:1364528/150‑1 (MQ=255)
                                                                                                           cAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAg                                     <  1:33223/150‑1 (MQ=255)
                                                                                                           cAGCGCCAGGCCGTAAGCGCCGTGCACGCGCCTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATATCCCAGGCCCTTCAg                                     <  2:881921/150‑1 (MQ=255)
                                                                                                            aGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGc                                    <  1:1196912/150‑1 (MQ=255)
                                                                                                                ccAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGGCGGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCtt                                >  1:12713/1‑150 (MQ=255)
                                                                                                                    gCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTcgcg                            >  2:203512/1‑150 (MQ=255)
                                                                                                                     ccGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAgg                                                                                             >  1:143124/1‑84 (MQ=255)
                                                                                                                     ccGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAgg                                                                                             <  2:143124/84‑1 (MQ=255)
                                                                                                                     ccGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGaa                                                          >  1:418229/1‑119 (MQ=255)
                                                                                                                     ccGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGaa                                                          <  2:418229/119‑1 (MQ=255)
                                                                                                                                gTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCAAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATg                >  1:575043/1‑150 (MQ=255)
                                                                                                                                 tGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCGGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATgg               <  1:659676/150‑1 (MQ=255)
                                                                                                                                     ggCGCTTCACCGCCGCCGTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTtct           <  2:227046/150‑1 (MQ=255)
                                                                                                                                        gcTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCga        <  1:497777/150‑1 (MQ=255)
                                                                                                                                           tgaCCGACGCCCTTAAGGCGGCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGct                                   <  2:443326/118‑1 (MQ=255)
                                                                                                                                           tCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCgatga     <  2:1203928/150‑1 (MQ=255)
                                                                                                                                            cACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCgatgat    <  1:756044/150‑1 (MQ=255)
                                                                                                                                             aCCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGc                                                  >  1:843173/1‑103 (MQ=255)
                                                                                                                                             aCCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGc                                                  <  2:843173/103‑1 (MQ=255)
                                                                                                                                              ccgccgCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCGATGATGc  >  1:1090495/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTCATGGTTCTCGATGATGC  >  NC_002947/6171765‑6172056

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