Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 963,971 C→A 43.0% F48L (TTC→TTA ihfB → integration host factor (IHF), DNA‑binding protein, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913963,9710CA43.0% 84.8 / 199.4 151F48L (TTC→TTAihfBintegration host factor (IHF), DNA‑binding protein, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (36/50);  new base A (37/28);  total (73/78)
Fisher's exact test for biased strand distribution p-value = 7.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

Only 100 of 153 total aligned reads displayed.

CACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAA  >  NC_000913/963873‑964070
                                                                                                  |                                                                                                   
cacaTTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCg                                                                                                    <  2:1100854/100‑1 (MQ=255)
cacaTTCCCGCCAAGACGGTTGAAGAGGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCg                                                                                                    <  2:3671251/100‑1 (MQ=255)
 acaTTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCgg                                                                                                   <  1:3981681/100‑1 (MQ=255)
  caTTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGc                                                                                                  >  1:1677737/1‑100 (MQ=255)
      cccGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGtt                                                                                              <  2:3691102/100‑1 (MQ=255)
        cGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTtc                                                                                            >  1:1876319/1‑100 (MQ=255)
        cGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGAGTATTGAAATCCGCGGTTTAGGCAGttt                                                                                             >  1:1876494/1‑99 (MQ=255)
           cAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTctct                                                                                         >  2:510742/1‑100 (MQ=255)
            aaGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCtt                                                                                        <  1:4013482/100‑1 (MQ=255)
             aGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCttt                                                                                       <  2:386729/100‑1 (MQ=255)
               aCGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGc                                                                                     <  2:1338911/100‑1 (MQ=255)
               aCGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGc                                                                                     <  1:2476057/100‑1 (MQ=255)
                cGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTg                                                                                      <  2:126531/98‑1 (MQ=255)
                cGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGAGCAGGGCGAGCGTATTGAAATCAGCGGTTTCGGCAGTTTCTCTTTGCa                                                                                    >  1:2837123/1‑100 (MQ=255)
                 ggTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCAc                                                                                   >  1:146166/1‑100 (MQ=255)
                  gTTGAAGATGCAGTAAAAGAGCTGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCAc                                                                                   >  2:2696987/1‑99 (MQ=255)
                  gTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACt                                                                                  >  1:2657549/1‑100 (MQ=255)
                    tGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTAc                                                                                >  2:1771567/1‑100 (MQ=255)
                    tGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTAc                                                                                >  1:2551015/1‑100 (MQ=255)
                     gAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTAcc                                                                               >  2:2257087/1‑100 (MQ=255)
                     gAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACGCCTGCGCAGGGCGAGCGCATTGAACTCAGCGGCTTAGGCCGCTTCCCTTCACACTAcc                                                                               >  1:3839777/1‑100 (MQ=255)
                      aaGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACcg                                                                              >  2:1835772/1‑100 (MQ=255)
                       aGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACcgc                                                                             <  2:133811/100‑1 (MQ=255)
                        gATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACcgcg                                                                            >  2:3318421/1‑100 (MQ=255)
                        gATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACcgcg                                                                            >  2:2074872/1‑100 (MQ=255)
                         aTGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCgcgc                                                                           <  1:2403000/100‑1 (MQ=255)
                          tGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACcgcg                                                                            >  1:3275301/1‑98 (MQ=255)
                          tGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCgcgc                                                                           >  1:4114229/1‑99 (MQ=255)
                          tGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGca                                                                          >  1:340825/1‑100 (MQ=255)
                          tGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGca                                                                          <  2:2360632/100‑1 (MQ=255)
                          tGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGca                                                                          >  2:3903485/1‑100 (MQ=255)
                          tGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACCACCGCGca                                                                          >  2:3087391/1‑100 (MQ=255)
                           gCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGcac                                                                         <  1:2829682/100‑1 (MQ=255)
                            cAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGCAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACAGCGcacc                                                                        >  1:3386656/1‑100 (MQ=255)
                            cAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGcacc                                                                        >  2:51325/1‑100 (MQ=255)
                            cAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGcacc                                                                        >  2:917949/1‑100 (MQ=255)
                             aGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGcacca                                                                       >  2:1199306/1‑100 (MQ=255)
                               tAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACg                                                                     >  2:1866013/1‑100 (MQ=255)
                                aaaaGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGt                                                                    <  2:498466/100‑1 (MQ=255)
                                 aaaGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTa                                                                   >  1:1528988/1‑100 (MQ=255)
                                    gagaTGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTAcc                                                                 <  2:1705610/99‑1 (MQ=255)
                                    gagaTGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCg                                                                >  2:3989025/1‑100 (MQ=255)
                                     agaTGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTAcccg                                                               >  2:3363188/1‑98 (MQ=255)
                                       aTGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGAc                                                             <  2:1541800/100‑1 (MQ=255)
                                        tGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACg                                                            <  2:2759207/100‑1 (MQ=255)
                                        tGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACg                                                            >  2:1557659/1‑100 (MQ=255)
                                        tGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACg                                                            >  2:1469186/1‑100 (MQ=255)
                                        tGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACg                                                            >  2:3974325/1‑100 (MQ=255)
                                         gCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGt                                                           >  1:1037999/1‑100 (MQ=255)
                                         gCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGt                                                           <  1:95646/100‑1 (MQ=255)
                                          cTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTa                                                          <  2:4131610/100‑1 (MQ=255)
                                           tGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTaa                                                         >  1:2502512/1‑100 (MQ=255)
                                           tGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTaa                                                         >  1:2905845/1‑100 (MQ=255)
                                            ggAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAAt                                                        >  2:3433799/1‑100 (MQ=255)
                                                 atatGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGaa                                                   <  1:989283/100‑1 (MQ=255)
                                                 atatGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTCTCTCTTTGCACTACCGCGCACCACGTACTGGACGTAATccta                                                    <  2:2895021/99‑3 (MQ=255)
                                                  tatGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAg                                                  <  2:121647/100‑1 (MQ=255)
                                                  tatGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAg                                                  >  1:4169554/1‑100 (MQ=255)
                                                  tatGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAg                                                  >  2:621044/1‑100 (MQ=255)
                                                    tGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGAc                                                >  2:641979/1‑100 (MQ=255)
                                                     ggCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACt                                               <  1:1670549/100‑1 (MQ=255)
                                                      gCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACt                                               <  1:2696977/99‑1 (MQ=255)
                                                       ccTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTgg                                             <  1:1434361/100‑1 (MQ=255)
                                                        cTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGc                                            <  2:3778576/100‑1 (MQ=255)
                                                        cTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGc                                            >  1:1311102/1‑100 (MQ=255)
                                                               ctTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAg                                     <  2:2954831/100‑1 (MQ=255)
                                                               ctTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAg                                     >  1:904811/1‑100 (MQ=255)
                                                                ttGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGt                                    <  2:2657558/100‑1 (MQ=255)
                                                                  gcgcAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAg                                  >  1:43328/1‑100 (MQ=255)
                                                                  gcgcAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAg                                  <  1:4216906/100‑1 (MQ=255)
                                                                  gcgcAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAg                                  >  2:2568839/1‑100 (MQ=255)
                                                                   cgcAGGGCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGa                                 >  2:1052603/1‑100 (MQ=255)
                                                                     cAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAAc                               >  2:2123296/1‑100 (MQ=255)
                                                                     cAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAAc                               <  1:2123291/100‑1 (MQ=255)
                                                                       gggCGAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTg                             >  1:2305018/1‑100 (MQ=255)
                                                                        ggCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTgg                            <  1:477156/100‑1 (MQ=255)
                                                                         gCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTgga                           <  2:1893551/100‑1 (MQ=255)
                                                                           gAGCGTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTggaag                         <  1:1306977/100‑1 (MQ=255)
                                                                            aGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTggaagg                        <  1:802032/100‑1 (MQ=255)
                                                                               gTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGaaa                     >  2:2378724/1‑100 (MQ=255)
                                                                               gTATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGaaa                     <  1:726318/100‑1 (MQ=255)
                                                                                tATTGAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGaaaa                    >  1:139148/1‑100 (MQ=255)
                                                                                 aTTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAAt                   >  1:1844102/1‑100 (MQ=255)
                                                                                    gAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACg                >  2:2110942/1‑100 (MQ=255)
                                                                                    gAAATCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACg                >  2:1430422/1‑100 (MQ=255)
                                                                                     aaaTCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGt               <  1:1780755/100‑1 (MQ=255)
                                                                                      aaTCCGCGGTTTAGGCAGTTTCTCTTTGCAATACCGCGCACCACGTACCGGACGTAATACGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGtt              >  1:2018722/1‑100 (MQ=255)
                                                                                       aTCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTc             <  2:1144/100‑1 (MQ=255)
                                                                                        tCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTcc            <  1:2463676/100‑1 (MQ=255)
                                                                                        tCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTcc            <  2:2265956/100‑1 (MQ=255)
                                                                                        tCCGCGGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTcc            <  1:4224578/100‑1 (MQ=255)
                                                                                          cgcgGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTc          <  2:1038002/100‑1 (MQ=255)
                                                                                          cgcgGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTc          <  2:3791862/100‑1 (MQ=255)
                                                                                            cgGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCGGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCAc        >  1:2873184/1‑100 (MQ=255)
                                                                                            cgGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCAc        <  1:1012525/100‑1 (MQ=255)
                                                                                            cgGTTTAGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCAc        <  2:1170430/100‑1 (MQ=255)
                                                                                              gTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACtt      <  2:324811/100‑1 (MQ=255)
                                                                                                ttCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATCAAGTAGCACTGGAACGAAAATACGTTCCTCACTTTa    >  1:2395091/1‑100 (MQ=255)
                                                                                                 tCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTaa   <  2:172453/100‑1 (MQ=255)
                                                                                                  cGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTaaa  >  1:301343/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAA  >  NC_000913/963873‑964070

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