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.

CTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATA  >  NC_000913/1818444‑1818637
                                                                                                     |                                                                                               
cTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTg                                                                                                  >  1:2004989/1‑100 (MQ=255)
cTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTg                                                                                                  >  2:505900/1‑100 (MQ=255)
cTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTg                                                                                                  >  1:761479/1‑100 (MQ=255)
  ccGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑cc                                                                                              <  1:250627/100‑1 (MQ=255)
  ccGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑cc                                                                                              <  1:2090910/100‑1 (MQ=255)
  ccGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑cc                                                                                              <  1:2918505/100‑1 (MQ=255)
  ccGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑cc                                                                                              <  2:2264645/100‑1 (MQ=255)
  ccGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑cc                                                                                              <  2:2936942/100‑1 (MQ=255)
      ggTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑cccccc                                                                                          <  1:2236008/100‑1 (MQ=255)
        tataGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑cccccccc                                                                                        <  1:1163630/100‑2 (MQ=255)
         ataGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑ccccccccg                                                                                       >  2:1007159/1‑98 (MQ=255)
         ataGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑cccccccc                                                                                        >  2:1976845/1‑98 (MQ=255)
          taGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTg‑‑cccgcccc                                                                                        >  2:2179245/1‑97 (MQ=255)
          taGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑ccccccccgg                                                                                      <  1:2771054/100‑4 (MQ=255)
              tgctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACat                                                                                   >  1:1777984/1‑100 (MQ=255)
               gctCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑ccccccccgga                                                                                     >  2:1949932/1‑92 (MQ=255)
                 tCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACataaa                                                                                <  2:2004977/100‑1 (MQ=255)
                  cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACataaat                                                                               >  2:2059261/1‑100 (MQ=255)
                    cGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACataaataa                                                                             >  2:445734/1‑100 (MQ=255)
                     gCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAaataat                                                                            <  2:1884297/100‑1 (MQ=255)
                     gCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAaataat                                                                            >  2:103269/1‑100 (MQ=255)
                       ctccaATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATcc                                                                          <  1:1406403/98‑1 (MQ=255)
                       caccaATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATcc                                                                          >  1:1391399/1‑100 (MQ=255)
                       caccaATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATcc                                                                          >  2:280621/1‑100 (MQ=255)
                       caccaATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGAGATCTGGCATATCTGCCCCCCCCGGACATAAATAATcc                                                                          <  2:1604866/100‑1 (MQ=255)
                        accaATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCa                                                                         >  1:757454/1‑100 (MQ=255)
                        accaATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCa                                                                         <  1:1220546/100‑1 (MQ=255)
                         ccaATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAg                                                                        >  2:1406395/1‑100 (MQ=255)
                          caATAGTGACGACTTTTAATTTCTGTATCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGc                                                                       <  2:2702727/100‑1 (MQ=255)
                          caATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGc                                                                       >  2:1709917/1‑100 (MQ=255)
                          caATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGc                                                                       <  1:1478872/100‑1 (MQ=255)
                          caATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGc                                                                       <  2:222528/100‑1 (MQ=255)
                          caATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGc                                                                       <  1:2156698/100‑1 (MQ=255)
                          caATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGc                                                                       <  1:113906/100‑1 (MQ=255)
                          caATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTCCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCCGACATAAATAACCCAGc                                                                       <  2:1343929/100‑1 (MQ=255)
                           aaTAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCa                                                                      <  1:2125519/100‑1 (MQ=255)
                           aaTAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCa                                                                      >  1:1644134/1‑100 (MQ=255)
                              aGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACa                                                                   <  1:1076643/100‑1 (MQ=255)
                               gTGACGACTTTTAATTTCTGGCTCATCATTTACCCCCTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCCGGACATAAATACTCCAGCAACa                                                                   <  1:2990097/100‑1 (MQ=255)
                                  acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                               <  1:662579/100‑1 (MQ=255)
                                  acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                               <  1:196346/100‑1 (MQ=255)
                                  acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                               <  1:103271/100‑1 (MQ=255)
                                  acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                                               <  2:1411867/100‑1 (MQ=255)
                                   cgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACa                                                              >  2:44777/1‑100 (MQ=255)
                                   cgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACa                                                              <  1:1851703/100‑1 (MQ=255)
                                    gacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAg                                                             <  2:2937580/100‑1 (MQ=255)
                                       ttttAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGata                                                          <  1:380946/100‑1 (MQ=255)
                                       ttttAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGata                                                          >  2:2423368/1‑100 (MQ=255)
                                        tttAATTTCTGGCTCATAATTTCCCCCTTCAGTACAGAATAATGAAATCTGGCAAAACTGCCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGata                                                          <  2:1159899/100‑1 (MQ=255)
                                            aTTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGa                                                     >  2:871519/1‑100 (MQ=255)
                                             tttCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGaa                                                    <  2:2626932/100‑1 (MQ=255)
                                               tCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAAtt                                                  <  2:499286/100‑1 (MQ=255)
                                                cTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTg                                                 >  1:2595788/1‑100 (MQ=255)
                                                  ggCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTc                                               >  1:1768025/1‑100 (MQ=255)
                                                  ggCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTc                                               <  2:1056902/100‑1 (MQ=255)
                                                  ggCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTc                                               >  2:783172/1‑100 (MQ=255)
                                                   gCTTATTATTTTTTCCTTCAGTACAAAATTATTGTAATCGGCATAACTGCCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTg                                                 <  2:1228899/98‑1 (MQ=255)
                                                   gCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCa                                              <  2:963957/100‑1 (MQ=255)
                                                     tCATAATTTCTCCGTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAgg                                            >  1:2004096/1‑100 (MQ=255)
                                                     tCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAgg                                            <  1:919707/100‑1 (MQ=255)
                                                      cATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGt                                           >  2:822446/1‑100 (MQ=255)
                                                       aTAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGta                                          >  1:1216697/1‑100 (MQ=255)
                                                       aTAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGta                                          <  2:761476/100‑1 (MQ=255)
                                                        tAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGtat                                         >  1:668906/1‑100 (MQ=255)
                                                        tAATTTCTCCCTTCAGTACAGAATACTGATATCGGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAAGTGTCAGGtat                                         >  2:1962929/1‑100 (MQ=255)
                                                         aaTTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGtata                                        <  2:2835668/100‑1 (MQ=255)
                                                                   ttCTGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACt                              <  2:746470/100‑1 (MQ=255)
                                                                   ttCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACt                              <  1:2423383/100‑1 (MQ=255)
                                                                    tCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGGCGTAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTg                             <  1:974109/100‑1 (MQ=255)
                                                                      aGTACAGAATACTGATATCTGGCATATCTG‑‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGAtt                          >  1:1763902/1‑100 (MQ=255)
                                                                       gTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGAtt                          <  1:384818/100‑1 (MQ=255)
                                                                        tACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGAtt                          <  1:2059273/99‑1 (MQ=255)
                                                                        tACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTa                         <  1:871522/100‑1 (MQ=255)
                                                                         aCAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTa                         <  2:706362/99‑1 (MQ=255)
                                                                          cAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTaaa                       >  1:532135/1‑100 (MQ=255)
                                                                          cAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTaaa                       <  2:2699456/100‑1 (MQ=255)
                                                                           aGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAt                      >  1:1099391/1‑100 (MQ=255)
                                                                            gAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt                     <  1:1590757/100‑1 (MQ=255)
                                                                             aaTACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTc                    <  2:705488/100‑1 (MQ=255)
                                                                              aTACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCa                   <  2:979548/100‑1 (MQ=255)
                                                                                tACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAg                  >  1:1656194/1‑100 (MQ=255)
                                                                                tACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAg                  <  2:1627995/100‑1 (MQ=255)
                                                                                 aCTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGt                 <  1:2623344/100‑1 (MQ=255)
                                                                                 aCTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGt                 <  2:2110217/100‑1 (MQ=255)
                                                                                 aCTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGt                 <  1:929302/100‑1 (MQ=255)
                                                                                  cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTaaagatcgga                >  2:916567/1‑93 (MQ=255)
                                                                                  cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTc                >  2:2605823/1‑100 (MQ=255)
                                                                                  cTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTc                <  2:1768014/100‑1 (MQ=255)
                                                                                    gATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCaa              >  2:1166395/1‑100 (MQ=255)
                                                                                        tCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTc          >  1:2253738/1‑100 (MQ=255)
                                                                                         cTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCt         >  2:2793095/1‑100 (MQ=255)
                                                                                         cTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCt         >  2:2075324/1‑100 (MQ=255)
                                                                                          tGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCtt        >  1:1957976/1‑100 (MQ=255)
                                                                                          tGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCtt        <  2:2607578/100‑1 (MQ=255)
                                                                                           ggCATATCTGCCCCCCCCGGACATAAATAATCCAGCAATAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTa       >  2:1835505/1‑100 (MQ=255)
                                                                                           ggCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTa       >  1:563931/1‑100 (MQ=255)
                                                                                             cATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATCTGTGTACTGACAGCTATAACCACTTACTGATTAAATTCATTCAATTTCTTACg     >  2:2542729/1‑100 (MQ=255)
                                                                                             cATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTAc      <  1:783175/99‑1 (MQ=255)
                                                                                              atatCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGa    >  2:2843548/1‑100 (MQ=255)
                                                                                                atCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCAGTCAATTTCTTACGATa  <  2:1101543/100‑1 (MQ=255)
                                                                                                     |                                                                                               
CTCCGGGGTATAGCTGCTCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑‑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: