Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,818,550 (C)7→8 intergenic (‑50/+55) chbF ← / ← chbR phospho‑chitobiase; general 6‑phospho‑beta‑glucosidase activity/repressor of chb operon for N,N'‑diacetylchitobiose utilization

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,818,5431.C97.7% 454.7 / 3.2 129intergenic (‑43/+62)chbF/chbRphospho‑chitobiase; general 6‑phospho‑beta‑glucosidase activity/repressor of chb operon for N,N'‑diacetylchitobiose utilization
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (61/65);  total (62/67)
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.82e-01

Only 100 of 148 total aligned reads displayed.

CTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTG‑‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATA  >  NC_000913/1818444‑1818637
                                                                                                     |                                                                                               
cTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTg                                                                                                  >  1:2004989/1‑100 (MQ=255)
cTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTg                                                                                                  >  2:505900/1‑100 (MQ=255)
cTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTg                                                                                                  >  1:761479/1‑100 (MQ=255)
  ccGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTG‑‑cc                                                                                              <  1:250627/100‑1 (MQ=255)
  ccGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTG‑‑cc                                                                                              <  1:2090910/100‑1 (MQ=255)
  ccGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTG‑‑cc                                                                                              <  1:2918505/100‑1 (MQ=255)
  ccGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTG‑‑cc                                                                                              <  2:2264645/100‑1 (MQ=255)
  ccGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTG‑‑cc                                                                                              <  2:2936942/100‑1 (MQ=255)
      ggTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTG‑‑cccccc                                                                                          <  1:2236008/100‑1 (MQ=255)
        tataGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTG‑‑cccccccc                                                                                        <  1:1163630/100‑2 (MQ=255)
         ataGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTG‑‑ccccccccg                                                                                       >  2:1007159/1‑98 (MQ=255)
         ataGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTG‑‑cccccccc                                                                                        >  2:1976845/1‑98 (MQ=255)
          taGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTg‑‑cccgcccc                                                                                        >  2:2179245/1‑97 (MQ=255)
          taGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTG‑‑ccccccccgg                                                                                      <  1:2771054/100‑4 (MQ=255)
              tgctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACat                                                                                   >  1:1777984/1‑100 (MQ=255)
               gctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTG‑‑ccccccccgga                                                                                     >  2:1949932/1‑92 (MQ=255)
                 tCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACataaa                                                                                <  2:2004977/100‑1 (MQ=255)
                  cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACataaat                                                                               >  2:2059261/1‑100 (MQ=255)
                    cGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACataaataa                                                                             >  2:445734/1‑100 (MQ=255)
                     gCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAaataat                                                                            <  2:1884297/100‑1 (MQ=255)
                     gCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAaataat                                                                            >  2:103269/1‑100 (MQ=255)
                       ctccaATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATcc                                                                          <  1:1406403/98‑1 (MQ=255)
                       caccaATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATcc                                                                          >  1:1391399/1‑100 (MQ=255)
                       caccaATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATcc                                                                          >  2:280621/1‑100 (MQ=255)
                       caccaATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGAGATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATcc                                                                          <  2:1604866/100‑1 (MQ=255)
                        accaATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCa                                                                         >  1:757454/1‑100 (MQ=255)
                        accaATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCa                                                                         <  1:1220546/100‑1 (MQ=255)
                         ccaATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAg                                                                        >  2:1406395/1‑100 (MQ=255)
                          caATAGTGACGACTTTTAATTTCTGTATCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGc                                                                       <  2:2702727/100‑1 (MQ=255)
                          caATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGc                                                                       >  2:1709917/1‑100 (MQ=255)
                          caATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGc                                                                       <  1:1478872/100‑1 (MQ=255)
                          caATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGc                                                                       <  2:222528/100‑1 (MQ=255)
                          caATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGc                                                                       <  1:2156698/100‑1 (MQ=255)
                          caATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGc                                                                       <  1:113906/100‑1 (MQ=255)
                          caATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTCCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCCGACATAAATAACCCAGc                                                                       <  2:1343929/100‑1 (MQ=255)
                           aaTAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCa                                                                      <  1:2125519/100‑1 (MQ=255)
                           aaTAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCa                                                                      >  1:1644134/1‑100 (MQ=255)
                              aGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACa                                                                   <  1:1076643/100‑1 (MQ=255)
                               gTGACGACTTTTAATTTCTGGCTCATCATTTACCCCCTCAGTACAGAA‑TACTGATATCTGGCATATCTGCCCCCCCCCGGACATAAATACTCCAGCAACa                                                                   <  1:2990097/100‑1 (MQ=255)
                                  acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                               <  1:662579/100‑1 (MQ=255)
                                  acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                               <  1:196346/100‑1 (MQ=255)
                                  acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                               <  1:103271/100‑1 (MQ=255)
                                  acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                               <  2:1411867/100‑1 (MQ=255)
                                   cgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACa                                                              >  2:44777/1‑100 (MQ=255)
                                   cgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACa                                                              <  1:1851703/100‑1 (MQ=255)
                                    gacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAg                                                             <  2:2937580/100‑1 (MQ=255)
                                       ttttAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGata                                                          <  1:380946/100‑1 (MQ=255)
                                       ttttAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGata                                                          >  2:2423368/1‑100 (MQ=255)
                                        tttAATTTCTGGCTCATAATTTCCCCCTTCAGTACAGAA‑TAATGAAATCTGGCAAAACTGCCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGata                                                          <  2:1159899/100‑1 (MQ=255)
                                            aTTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGa                                                     >  2:871519/1‑100 (MQ=255)
                                             tttCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGaa                                                    <  2:2626932/100‑1 (MQ=255)
                                               tCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAAtt                                                  <  2:499286/100‑1 (MQ=255)
                                                cTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTg                                                 >  1:2595788/1‑100 (MQ=255)
                                                  ggCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTc                                               >  1:1768025/1‑100 (MQ=255)
                                                  ggCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTc                                               <  2:1056902/100‑1 (MQ=255)
                                                  ggCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTc                                               >  2:783172/1‑100 (MQ=255)
                                                   gCTTATTATTTTTTCCTTCAGTACAAAATTATTGTAATC‑GGCATAACTGCCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTg                                                 <  2:1228899/98‑1 (MQ=255)
                                                   gCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCa                                              <  2:963957/100‑1 (MQ=255)
                                                     tCATAATTTCTCCGTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAgg                                            >  1:2004096/1‑100 (MQ=255)
                                                     tCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAgg                                            <  1:919707/100‑1 (MQ=255)
                                                      cATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGt                                           >  2:822446/1‑100 (MQ=255)
                                                       aTAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGta                                          >  1:1216697/1‑100 (MQ=255)
                                                       aTAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGta                                          <  2:761476/100‑1 (MQ=255)
                                                        tAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGtat                                         >  1:668906/1‑100 (MQ=255)
                                                        tAATTTCTCCCTTCAGTACAGAA‑TACTGATATCGGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAAGTGTCAGGtat                                         >  2:1962929/1‑100 (MQ=255)
                                                         aaTTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGtata                                        <  2:2835668/100‑1 (MQ=255)
                                                                   ttCTGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACt                              <  2:746470/100‑1 (MQ=255)
                                                                   ttCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACt                              <  1:2423383/100‑1 (MQ=255)
                                                                    tCAGTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGGCGTAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTg                             <  1:974109/100‑1 (MQ=255)
                                                                      aGTACAGAA‑TACTGATATCTGGCATATCTG‑‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGAtt                          >  1:1763902/1‑100 (MQ=255)
                                                                       gTACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGAtt                          <  1:384818/100‑1 (MQ=255)
                                                                        tACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGAtt                          <  1:2059273/99‑1 (MQ=255)
                                                                        tACAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTa                         <  1:871522/100‑1 (MQ=255)
                                                                         aCAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTa                         <  2:706362/99‑1 (MQ=255)
                                                                          cAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTaaa                       >  1:532135/1‑100 (MQ=255)
                                                                          cAGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTaaa                       <  2:2699456/100‑1 (MQ=255)
                                                                           aGAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAt                      >  1:1099391/1‑100 (MQ=255)
                                                                            gAA‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt                     <  1:1590757/100‑1 (MQ=255)
                                                                             aa‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTc                    <  2:705488/100‑1 (MQ=255)
                                                                              a‑TACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCa                   <  2:979548/100‑1 (MQ=255)
                                                                                tACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAg                  >  1:1656194/1‑100 (MQ=255)
                                                                                tACTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAg                  <  2:1627995/100‑1 (MQ=255)
                                                                                 aCTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGt                 <  1:2623344/100‑1 (MQ=255)
                                                                                 aCTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGt                 <  2:2110217/100‑1 (MQ=255)
                                                                                 aCTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGt                 <  1:929302/100‑1 (MQ=255)
                                                                                  cTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTaaagatcgga                >  2:916567/1‑93 (MQ=255)
                                                                                  cTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTc                >  2:2605823/1‑100 (MQ=255)
                                                                                  cTGATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTc                <  2:1768014/100‑1 (MQ=255)
                                                                                    gATATCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCaa              >  2:1166395/1‑100 (MQ=255)
                                                                                        tCTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTc          >  1:2253738/1‑100 (MQ=255)
                                                                                         cTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCt         >  2:2793095/1‑100 (MQ=255)
                                                                                         cTGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCt         >  2:2075324/1‑100 (MQ=255)
                                                                                          tGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCtt        >  1:1957976/1‑100 (MQ=255)
                                                                                          tGGCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCtt        <  2:2607578/100‑1 (MQ=255)
                                                                                           ggCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAATAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTa       >  2:1835505/1‑100 (MQ=255)
                                                                                           ggCATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTa       >  1:563931/1‑100 (MQ=255)
                                                                                             cATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATCTGTGTACTGACAGCTATAACCACTTACTGATTAAATTCATTCAATTTCTTACg     >  2:2542729/1‑100 (MQ=255)
                                                                                             cATATCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTAc      <  1:783175/99‑1 (MQ=255)
                                                                                              atatCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGa    >  2:2843548/1‑100 (MQ=255)
                                                                                                atCTGC‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATa  <  2:1101543/100‑1 (MQ=255)
                                                                                                     |                                                                                               
CTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAA‑TACTGATATCTGGCATATCTG‑‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATA  >  NC_000913/1818444‑1818637

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