Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,957,957 C→T intergenic (‑121/+78) ppiC ← / ← yifN peptidyl‑prolyl cis‑trans isomerase C (rotamase C)/conserved protein (pseudogene)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,957,9570CT99.8% 1674.2 / NA 442intergenic (‑121/+78)ppiC/yifNpeptidyl‑prolyl cis‑trans isomerase C (rotamase C)/conserved protein (pseudogene)
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (143/298);  total (144/298)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 524 total aligned reads displayed.

TGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCCTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGG  >  NC_000913/3957859‑3958054
                                                                                                  |                                                                                                 
tctCAGGCGGCGATTTTAGCATGATCCGCAGATTCTAACGTAATTGTCAGTCAGGGCGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTc                                                                                                 <  2:7001111/99‑1 (MQ=255)
tGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGGCCTTTTTGGGAAAAGCCCGGAAAATCCCGGCTTc                                                                                                 >  1:12674238/1‑101 (MQ=255)
  tCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGATTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGc                                                                                               <  1:13248879/101‑1 (MQ=255)
    aGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGGCATATATGGCAAAAACCCGGAAAATCCCGGCTTCgnnn                                                                                             >  2:12248451/1‑98 (MQ=255)
     nnnggcgATTTGAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGGCCGGAAAATCCGGGCTTCGCAgg                                                                                            <  2:10612639/98‑1 (MQ=255)
       cggcgATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTc                                                                                          <  2:1979854/101‑1 (MQ=255)
       cggcgATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAAACCGGGCTTCGCAGGTc                                                                                          >  1:14231214/1‑101 (MQ=255)
         gcgATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGACCGGAAAAACCGGGCTTCGCAGGtnnn                                                                                        >  2:11108999/1‑98 (MQ=255)
          nnnttttAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCt                                                                                       <  2:13699021/98‑1 (MQ=255)
           gATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTCCATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTc                                                                                      <  1:10707531/101‑1 (MQ=255)
           gATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGGAAATCCGGGCTTCGCAGGTCCCTc                                                                                      >  1:7648984/1‑101 (MQ=255)
           gATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAAACCGGGCTTCGCAGGTCCCTc                                                                                      >  1:59993/1‑101 (MQ=255)
             ttttAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCaa                                                                                    <  2:3742343/101‑1 (MQ=255)
               ttAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgg                                                                                  <  1:11870526/101‑1 (MQ=255)
               nnngCATGATCCGCAGATTCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgg                                                                                  <  2:14167119/98‑1 (MQ=255)
               nnngCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgg                                                                                  <  2:10219655/98‑1 (MQ=255)
                 aGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgnnn                                                                                >  2:13928753/1‑98 (MQ=255)
                   gaTGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGgag                                                                              <  1:5642967/100‑1 (MQ=255)
                    aTGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGgaga                                                                             >  2:898267/1‑101 (MQ=255)
                    aTGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGgaga                                                                             <  1:7115568/101‑1 (MQ=255)
                    aTGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGgaga                                                                             <  2:1552408/101‑1 (MQ=255)
                      gATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCCGGCTTCGCAGGTCCCTCAAGGGGAGaaa                                                                           <  1:4760685/101‑1 (MQ=255)
                       aTCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAACCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaaaa                                                                          >  2:5614306/1‑101 (MQ=255)
                        tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTTGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAAc                                                                         <  1:12010837/101‑1 (MQ=255)
                        tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAACCCGGAAAATCCGGGCGTTGCAGGTCCCTCAAGGGGAGAAAAc                                                                         >  2:2874038/1‑101 (MQ=255)
                         ccGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACt                                                                        <  1:6178723/101‑1 (MQ=255)
                           gCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGGAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTa                                                                      >  2:5541297/1‑101 (MQ=255)
                            nnnaTGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                                     <  2:8814706/98‑1 (MQ=255)
                            cAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAACAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                                     <  1:2098459/101‑1 (MQ=255)
                            cAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                                                                     >  2:5980428/1‑101 (MQ=255)
                             nnntGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAgg                                                                    <  2:11545229/98‑1 (MQ=255)
                             aGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGGCCCGAAAATCCGGGGTTCGCAGGTCACTCAGAGGGAGGAAACTTAgg                                                                    >  2:6380318/1‑101 (MQ=255)
                              gATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTGCATATATGGGAAAAGCACGGAAAATCCGGGCTTCGGAGGTCCCTCAAGGGGAGAGAACTTannn                                                                   >  2:13167181/1‑98 (MQ=255)
                              gATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAggg                                                                   <  1:9406354/101‑1 (MQ=255)
                                 gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGACCGGAAAATCCGGGCTTCGCAGGGCCCTCAAGGGGAGAAAACTTAGGGGct                                                                >  1:677308/1‑101 (MQ=255)
                                  cTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcnnn                                                               >  2:12858399/1‑98 (MQ=255)
                                   tAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctct                                                              <  2:1792158/101‑1 (MQ=255)
                                   tAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACGTAGGGCctct                                                              >  2:6831458/1‑101 (MQ=255)
                                    nnngTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTa                                                             <  2:9532232/98‑1 (MQ=255)
                                    aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTa                                                             <  1:14601621/101‑1 (MQ=255)
                                     aCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAt                                                            <  1:6480023/101‑1 (MQ=255)
                                      catAATTGTCAGTCAGGGAGTGGCTGGTACATATCTGGGCAAAGCCCGGAAAATCCGGGGCTCGCAGGGCCCTCAAGGGGAGAAAACCTAGGGCCTCTATg                                                           >  2:3507790/3‑101 (MQ=255)
                                      cGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAACTCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATg                                                           <  2:3980052/101‑1 (MQ=255)
                                      cGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGGGAAAACTTAGGGCCTCTATg                                                           >  1:3854199/1‑101 (MQ=255)
                                      cGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctctnnn                                                           >  2:13316972/1‑98 (MQ=255)
                                        tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCCAGGGGGGGAAACTTAGGGCCTCTATggg                                                         >  1:2752339/1‑101 (MQ=255)
                                         aaTTGTCAGTCAGGGAGTGGCTGGTACCTATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGc                                                        >  1:4238804/1‑101 (MQ=255)
                                           ttGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTa                                                      <  1:7265487/101‑1 (MQ=255)
                                            tGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAACCCCGGAAAACCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                                     <  1:4973375/101‑1 (MQ=255)
                                             nnnagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                                    <  2:7811156/98‑1 (MQ=255)
                                             gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                                    >  1:7405864/1‑101 (MQ=255)
                                             gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                                    <  1:9827248/101‑1 (MQ=255)
                                               cagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                                                  <  1:1154639/101‑1 (MQ=255)
                                               cagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaaa                                                  <  1:1124953/101‑1 (MQ=255)
                                                 gtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGt                                                <  1:8010349/101‑1 (MQ=255)
                                                 gtcagGGAGTGGCTGGTACATATATGGCAAAAACCCGGAAAATCCCGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGt                                                >  1:12187153/1‑101 (MQ=255)
                                                   cagGGAGTGGCTGGTACATATATGTCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCa                                              <  2:4647828/101‑1 (MQ=255)
                                                   cagGGAGTGGCGGGTACATATATGGCAAAAGCCCGGAAAATCCGGGTTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCa                                              <  1:5377313/101‑1 (MQ=255)
                                                     gggAGGGGCTGGTACATATATGGCAAAACCCCGGGAAATACGGGGTTCTCAGGTCCCTCAAGGGGAGAAAACTTAAGTCCTCTATTGGCTACAAAGGCAtt                                            >  2:2249314/1‑101 (MQ=255)
                                                       gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTAt                                          >  1:3419192/1‑101 (MQ=255)
                                                        aGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATc                                         <  1:9691822/101‑1 (MQ=255)
                                                           nnntGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                                      <  2:8100691/98‑1 (MQ=255)
                                                           ggCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGCCATTATCggg                                      >  1:8698037/1‑101 (MQ=255)
                                                            gCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgggg                                     <  2:3422705/101‑1 (MQ=255)
                                                            gCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgggg                                     <  1:11108999/101‑1 (MQ=255)
                                                              tGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGAt                                   >  1:3550329/1‑101 (MQ=255)
                                                               ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATg                                  <  1:1441106/101‑1 (MQ=255)
                                                               ggTACATATATGGCAAAAGCCAGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATg                                  <  1:10203958/101‑1 (MQ=255)
                                                                gTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATgg                                 <  2:6164218/101‑1 (MQ=255)
                                                                     tatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttnnn                            >  2:12313856/1‑98 (MQ=255)
                                                                     tatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatt                            <  1:1999725/101‑1 (MQ=255)
                                                                       tatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttt                          <  1:6231982/101‑1 (MQ=255)
                                                                       tatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttt                          >  1:2260650/1‑101 (MQ=255)
                                                                       nnnggCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGCGGATGGttatttt                          <  2:11313411/98‑1 (MQ=255)
                                                                        nnngCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttca                         <  2:11821318/98‑3 (MQ=255)
                                                                        atGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttatttta                         <  2:5145319/101‑1 (MQ=255)
                                                                            nnnaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGt                     <  2:8530993/98‑1 (MQ=255)
                                                                             aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttattnnn                    >  2:9220399/1‑98 (MQ=255)
                                                                             aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTc                    <  1:11655983/101‑1 (MQ=255)
                                                                             aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTc                    >  1:14001086/1‑101 (MQ=255)
                                                                             aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTc                    <  2:3212411/101‑1 (MQ=255)
                                                                             aaaaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTc                    >  1:1604287/1‑101 (MQ=255)
                                                                               aaGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCaa                  <  1:898267/101‑1 (MQ=255)
                                                                                aGCCAGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAAt                 >  2:1211512/1‑101 (MQ=255)
                                                                                 gCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAAtt                <  1:12849390/101‑1 (MQ=255)
                                                                                   nnngAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTa              <  2:14305199/98‑1 (MQ=255)
                                                                                   nnngAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTa              <  2:7613786/98‑1 (MQ=255)
                                                                                    cGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCa                   <  2:2557587/95‑1 (MQ=255)
                                                                                     nnnaaaTCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGTATGGTTATTTTATTGTCAATTTATg            <  2:7329620/98‑1 (MQ=255)
                                                                                     ggAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATg            <  1:10858821/101‑1 (MQ=255)
                                                                                      gAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTGATTGTCAATTTATGa           <  1:13388517/101‑1 (MQ=255)
                                                                                       gaaaTCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAg          <  1:1371136/100‑1 (MQ=255)
                                                                                        aaaTCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGc         <  1:6905912/101‑1 (MQ=255)
                                                                                        aaaTCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGc         <  2:4327303/101‑1 (MQ=255)
                                                                                         aaTCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCg        <  1:10373453/101‑1 (MQ=255)
                                                                                             cGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGatat    <  1:826273/101‑1 (MQ=255)
                                                                                             cGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGatat    <  1:8489113/101‑1 (MQ=255)
                                                                                              gggCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATg   >  1:6638115/1‑101 (MQ=255)
                                                                                              gggCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATg   <  1:1299423/101‑1 (MQ=255)
                                                                                              gggCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGCGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATg   >  1:4518219/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
TGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCCTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGG  >  NC_000913/3957859‑3958054

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