Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,818,468 T→C A257A (GCA→GCG celB ← PTS system cellobiose‑specific enzyme IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,818,4680TC84.3% 591.7 / 100.6 261A257A (GCA→GCGcelBPTS system cellobiose‑specific enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base T (24/17);  new base C (129/91);  total (153/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 266 total aligned reads displayed.

ACCGGCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGAT  >  NC_000913/1818378‑1818558
                                                                                          |                                                                                          
aCCGGCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAgcgc                                                                                          >  1:777715/1‑93 (MQ=255)
 ccGGCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGcgcg                                                                                         >  1:490485/1‑93 (MQ=255)
 ccGGCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGcgcg                                                                                         >  1:39848/1‑93 (MQ=255)
 ccGGCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGcgcg                                                                                         <  2:6746840/93‑1 (MQ=255)
   ggCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGt                                                                                       <  2:7128271/93‑1 (MQ=255)
    gCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTc                                                                                      <  1:6240241/93‑1 (MQ=255)
      tGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTcag                                                                                    >  2:636725/1‑93 (MQ=255)
      tGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTcag                                                                                    >  1:7052538/1‑93 (MQ=255)
       gCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTcagc                                                                                   <  2:4398349/93‑1 (MQ=255)
       gCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTcagc                                                                                   <  2:1270513/93‑1 (MQ=255)
       gCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTcagc                                                                                   <  2:16451/93‑1 (MQ=255)
             gccgcTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTcagcgccagc                                                                             >  1:323135/1‑93 (MQ=255)
              ccgcTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTcagcgccagcg                                                                            >  2:119625/1‑93 (MQ=255)
               cgcTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTcagcgccagcgc                                                                           >  1:5962555/1‑93 (MQ=255)
                gcTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAgcgcg                                                                          >  1:1195553/1‑93 (MQ=255)
                 cTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAgcgcgc                                                                         >  1:2409210/1‑93 (MQ=255)
                  ttCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGcc                                                                        >  1:7398172/1‑93 (MQ=255)
                  ttCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGcc                                                                        >  2:3454418/1‑93 (MQ=255)
                   tCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCa                                                                       <  2:6902245/93‑1 (MQ=255)
                    cGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCAt                                                                      <  1:5086838/93‑1 (MQ=255)
                      aCGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGa                                                                    <  2:824373/93‑1 (MQ=255)
                      aCGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGa                                                                    >  1:7454839/1‑93 (MQ=255)
                       cGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGaa                                                                   <  1:2805311/93‑1 (MQ=255)
                       cGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGaa                                                                   >  1:6853541/1‑93 (MQ=255)
                        ggAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAAt                                                                  >  2:2117455/1‑93 (MQ=255)
                        ggAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAAt                                                                  >  1:2170255/1‑93 (MQ=255)
                        ggAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAAt                                                                  >  1:500125/1‑93 (MQ=255)
                         gAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATa                                                                 >  2:1788650/1‑93 (MQ=255)
                           aCCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATAcc                                                               >  1:2323640/1‑93 (MQ=255)
                            ccATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCg                                                              >  1:7211810/1‑93 (MQ=255)
                             cATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCga                                                             >  1:7231625/1‑93 (MQ=255)
                             cATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCga                                                             >  1:4032355/1‑93 (MQ=255)
                             cATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCga                                                             >  1:3029740/1‑93 (MQ=255)
                              atatTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCgaa                                                            >  2:362218/1‑93 (MQ=255)
                              atatTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCgaa                                                            <  2:4893367/93‑1 (MQ=255)
                               tatTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCgaag                                                           >  2:2851708/1‑93 (MQ=255)
                               tatTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCgaaa                                                           >  1:933591/1‑92 (MQ=255)
                                 ttGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGCTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCgaagaa                                                         >  1:5204369/1‑93 (MQ=255)
                                     tGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCag                                                     >  1:3495/1‑93 (MQ=255)
                                      gATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCaga                                                    <  2:5534548/93‑1 (MQ=255)
                                      gATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCaga                                                    >  2:5633169/1‑93 (MQ=255)
                                      gATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCaga                                                    <  2:1019997/93‑1 (MQ=255)
                                          ggTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAg                                                >  1:6031286/1‑93 (MQ=255)
                                             cgcgATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTgg                                             <  2:7455742/93‑1 (MQ=255)
                                             cgcgATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTgg                                             >  1:6673163/1‑93 (MQ=255)
                                               cgATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGaa                                           <  2:5147148/93‑1 (MQ=255)
                                               cgATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGaa                                           >  1:1905422/1‑93 (MQ=255)
                                               cgATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGaa                                           <  2:5787553/93‑1 (MQ=255)
                                                gATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAAc                                          >  1:1010816/1‑93 (MQ=255)
                                                 atatTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACa                                         <  2:7271698/93‑1 (MQ=255)
                                                   atTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACaaa                                       <  2:5566605/93‑1 (MQ=255)
                                                    ttttCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAg                                      >  2:3248252/1‑93 (MQ=255)
                                                    ttttCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAg                                      <  2:2238935/93‑1 (MQ=255)
                                                    ttttCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAg                                      >  1:3560358/1‑93 (MQ=255)
                                                     tttCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGa                                     <  2:1094810/93‑1 (MQ=255)
                                                      ttCCAGTGCCCACGGCGTCATAATGCCATTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGAt                                    >  1:5661314/1‑93 (MQ=255)
                                                         cAGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATcac                                 <  1:6628984/93‑1 (MQ=255)
                                                          aGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATcaca                                >  1:4407835/1‑93 (MQ=255)
                                                          aGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATcaca                                <  1:1022596/93‑1 (MQ=255)
                                                          aGTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATcaca                                <  2:6312588/93‑1 (MQ=255)
                                                           gagCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACAt                               >  1:4790883/3‑93 (MQ=255)
                                                           gTGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACAt                               >  1:1630217/1‑93 (MQ=255)
                                                            tGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATa                              <  1:269442/93‑1 (MQ=255)
                                                            tGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATa                              <  1:6326374/93‑1 (MQ=255)
                                                            tGCCCACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATa                              <  1:1804886/93‑1 (MQ=255)
                                                               ccACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGgc                           >  1:6094327/1‑93 (MQ=255)
                                                                cACGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGgcc                          <  1:4736119/93‑1 (MQ=255)
                                                                 aCGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGgccc                         >  1:3730524/1‑93 (MQ=255)
                                                                  cGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGgccca                        >  2:2665280/1‑93 (MQ=255)
                                                                  cGGCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGgccca                        >  2:3178044/1‑93 (MQ=255)
                                                                    gCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGgcccagc                      <  1:6764660/93‑1 (MQ=255)
                                                                    gCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGgcccagc                      >  1:4096699/1‑93 (MQ=255)
                                                                    gCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGgcccagc                      <  1:112473/93‑1 (MQ=255)
                                                                    gCGTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGgcccagc                      >  2:6004098/1‑93 (MQ=255)
                                                                      gTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGgcccagccc                    >  1:6486967/1‑93 (MQ=255)
                                                                      gTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGgcccagccc                    >  1:5308061/1‑93 (MQ=255)
                                                                      gTCATAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGgcccagccc                    <  1:5899254/93‑1 (MQ=255)
                                                                        cATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCccac                  >  1:541551/1‑93 (MQ=255)
                                                                         aTAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCccacc                 >  2:5219679/1‑93 (MQ=255)
                                                                         aTAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCccacc                 <  2:1912656/93‑1 (MQ=255)
                                                                          tAATGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCccacca                <  2:4209309/93‑1 (MQ=255)
                                                                           aaTGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCcaccac               >  2:5013337/1‑93 (MQ=255)
                                                                            aTGCCGTTGTCCAGCGCGGTCAGCGCCGGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACg              >  2:6428335/1‑93 (MQ=255)
                                                                             tGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGc             <  1:1263383/93‑1 (MQ=255)
                                                                             tGCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGc             <  2:6633616/93‑1 (MQ=255)
                                                                              gCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCt            >  1:1313470/1‑93 (MQ=255)
                                                                              gCCGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCt            >  2:3339160/1‑93 (MQ=255)
                                                                                cGTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTAc          <  2:1579204/93‑1 (MQ=255)
                                                                                 gTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTAcc         >  1:6457322/1‑93 (MQ=255)
                                                                                 gTTGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTAcc         >  2:3758271/1‑93 (MQ=255)
                                                                                   tGTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCa       <  2:3881509/93‑1 (MQ=255)
                                                                                    gTCCAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCaa      <  2:2738277/93‑1 (MQ=255)
                                                                                     tCCAGTGCGGTCAGCGCCGGCGCGCCATGAATACCGAAGAACCAGAGCCGTGGAACAAAGAACACATAGGCCCAGCCCACCACGCTACCCAAc     <  1:6912099/93‑1 (MQ=255)
                                                                                     tCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAAc     >  1:6708899/1‑93 (MQ=255)
                                                                                       cAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGa   <  2:6165466/93‑1 (MQ=255)
                                                                                       cAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGa   <  1:1381572/93‑1 (MQ=255)
                                                                                       cAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGa   <  1:3948677/93‑1 (MQ=255)
                                                                                       cAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGa   <  1:896146/93‑1 (MQ=255)
                                                                                       cAGCGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGa   <  2:4819999/93‑1 (MQ=255)
                                                                                        aGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGAt  <  1:4118656/93‑1 (MQ=255)
                                                                                          |                                                                                          
ACCGGCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGAT  >  NC_000913/1818378‑1818558

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