Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,820,144 Δ1 bp intergenic (+24/+18) b4542 → / ← oraA hypothetical protein/OraA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,820,1430G.99.1% 2525.4 / 10.6 536intergenic (+23/+19)b4542/oraAhypothetical protein/OraA
Reads supporting (aligned to +/- strand):  ref base G (5/0);  new base . (301/229);  total (306/230)
Fisher's exact test for biased strand distribution p-value = 7.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 574 total aligned reads displayed.

CAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATC  >  NC_000913/2820048‑2820237
                                                                                               |                                                                                              
 aaGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAgt                                                                                               >  1:5386840/1‑95 (MQ=255)
  aGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAgta                                                                                              >  2:8128328/1‑94 (MQ=255)
   gCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAgtaa                                                                                             >  1:5852406/1‑93 (MQ=255)
     tAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTaaaa                                                                                          >  2:5884061/1‑96 (MQ=255)
     tAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTaaaa                                                                                          >  1:4434438/1‑96 (MQ=255)
       aCGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATa                                                                                        <  2:12433459/96‑1 (MQ=255)
       aCGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATa                                                                                        <  2:6075600/96‑1 (MQ=255)
       aCGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATa                                                                                        >  1:1774378/1‑96 (MQ=255)
         gAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATAcc                                                                                      >  1:6272715/1‑96 (MQ=255)
         gAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGCGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATAcc                                                                                      >  2:9364584/1‑96 (MQ=255)
           aaTATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGt                                                                                    >  1:13957096/1‑96 (MQ=255)
            atatCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTa                                                                                   <  2:6885524/96‑1 (MQ=255)
               tagCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGGAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGc                                                                                >  1:13022841/3‑96 (MQ=255)
                   caccaGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTc                                                                            >  1:1017109/1‑96 (MQ=255)
                    accaGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGCGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCa                                                                           >  2:14403663/1‑96 (MQ=255)
                     ccaGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAg                                                                          >  2:2380623/1‑96 (MQ=255)
                        gCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCg                                                                       >  2:6590297/1‑96 (MQ=255)
                        gCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCg                                                                       <  2:13457767/96‑1 (MQ=255)
                         cTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCgg                                                                      <  2:11702622/96‑1 (MQ=255)
                          tCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGc                                                                     <  2:1270650/96‑1 (MQ=255)
                          tCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGc                                                                     >  1:8911627/1‑96 (MQ=255)
                            cAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCaa                                                                   <  2:10027365/96‑1 (MQ=255)
                            cAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCaa                                                                   >  2:901627/1‑96 (MQ=255)
                            cAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCaa                                                                   >  2:3736473/1‑96 (MQ=255)
                            aaGCGTGTCTTAATCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCaa                                                                   >  1:11417622/2‑96 (MQ=255)
                              gCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCaaaa                                                                 >  2:7428435/1‑96 (MQ=255)
                                 tgtCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAAttt                                                              <  1:696539/96‑1 (MQ=255)
                                  gtCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTc                                                             >  1:6523864/1‑96 (MQ=255)
                                   tCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCg                                                            >  1:8547736/1‑96 (MQ=255)
                                      tAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCa                                                         <  1:4233695/96‑1 (MQ=255)
                                      tAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCa                                                         <  2:3061942/96‑1 (MQ=255)
                                       aaCCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAGATTTCGCCaa                                                        >  2:7810411/1‑96 (MQ=255)
                                        aCCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCaaa                                                       <  1:1510644/96‑1 (MQ=255)
                                           gccgGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAtct                                                    >  2:5631240/1‑96 (MQ=255)
                                           gccgGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAtct                                                    <  2:14157475/96‑1 (MQ=255)
                                           gccgGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAtct                                                    >  2:10223862/1‑96 (MQ=255)
                                              gggCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCt                                                 >  2:11564511/1‑96 (MQ=255)
                                              gggCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCt                                                 <  2:2079383/96‑1 (MQ=255)
                                              gggCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCt                                                 >  2:4360071/1‑96 (MQ=255)
                                               ggCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTg                                                >  1:9509226/1‑96 (MQ=255)
                                                 cTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGa                                              >  2:3206760/1‑96 (MQ=255)
                                                 cTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGa                                              <  1:1190777/96‑1 (MQ=255)
                                                  tGGTAACTGAAAAGAGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGat                                             >  1:2990639/1‑96 (MQ=255)
                                                   ggTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGata                                            <  2:4062734/96‑1 (MQ=255)
                                                    gTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGatat                                           >  2:8595295/1‑96 (MQ=255)
                                                    gTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGatat                                           >  1:2118358/1‑96 (MQ=255)
                                                    gTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGatat                                           >  2:11473639/1‑96 (MQ=255)
                                                     tgaCTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATc                                          >  1:10643364/3‑96 (MQ=255)
                                                     tAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATc                                          >  1:9181579/1‑96 (MQ=255)
                                                     tAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATc                                          >  1:765745/1‑96 (MQ=255)
                                                      aaCTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCt                                         <  1:4013131/96‑1 (MQ=255)
                                                      aaCTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCt                                         <  1:3454576/96‑1 (MQ=255)
                                                       aCTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCtt                                        >  2:1092588/1‑96 (MQ=255)
                                                       aCTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCtt                                        <  1:9728429/96‑1 (MQ=255)
                                                       aCTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCtt                                        <  1:11198094/96‑1 (MQ=255)
                                                        cTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTc                                       >  1:12433350/1‑96 (MQ=255)
                                                          gAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCa                                     <  2:14303510/96‑1 (MQ=255)
                                                            aaaGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATc                                   >  1:7289549/1‑96 (MQ=255)
                                                             aaGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCa                                  >  1:1208746/1‑96 (MQ=255)
                                                              aGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAg                                 >  2:5242064/1‑96 (MQ=255)
                                                              aGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAg                                 >  2:7251712/1‑96 (MQ=255)
                                                              aGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAg                                 <  2:13621115/96‑1 (MQ=255)
                                                                 gggAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATCCCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATa                              <  1:2803540/96‑1 (MQ=255)
                                                                   gAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGc                            >  1:1279973/1‑96 (MQ=255)
                                                                   gAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGc                            >  1:7484438/1‑96 (MQ=255)
                                                                    aaTAAGATAAGTTTTCTTGACTGGGAAGGAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGcc                           >  2:3639287/1‑96 (MQ=255)
                                                                    aaTAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGcc                           >  2:12382126/1‑96 (MQ=255)
                                                                    aaTAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGcc                           >  2:1770112/1‑96 (MQ=255)
                                                                       aagataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACg                        >  2:8482658/1‑96 (MQ=255)
                                                                        agataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGa                       >  1:8279208/1‑96 (MQ=255)
                                                                         gataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGAt                      >  1:13288123/1‑96 (MQ=255)
                                                                          ataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATa                     <  2:10624405/96‑1 (MQ=255)
                                                                           taagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATag                    <  1:2011237/96‑1 (MQ=255)
                                                                           taagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATag                    >  2:11929077/1‑96 (MQ=255)
                                                                           taagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATag                    <  1:5088660/96‑1 (MQ=255)
                                                                            aagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATaga                   <  2:5386840/96‑1 (MQ=255)
                                                                             agTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATagag                  <  2:5003076/96‑1 (MQ=255)
                                                                             agTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATagag                  <  1:4484040/96‑1 (MQ=255)
                                                                               ttttCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCa                >  2:9719321/1‑96 (MQ=255)
                                                                                tttCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAg               <  2:2752088/96‑1 (MQ=255)
                                                                                tttCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAg               <  1:10967942/96‑1 (MQ=255)
                                                                                 ttCTTGACTGGGAAGGAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGa              >  1:1164256/1‑96 (MQ=255)
                                                                                  tCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGaa             <  1:4693104/96‑1 (MQ=255)
                                                                                   cTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGaaa            >  2:3816677/1‑96 (MQ=255)
                                                                                   cTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGaaa            <  2:13868474/96‑1 (MQ=255)
                                                                                   cTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGaaa            <  1:14238050/96‑1 (MQ=255)
                                                                                    ttGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGaaaa           <  2:476419/96‑1 (MQ=255)
                                                                                    gtGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGaaaa           >  1:13811148/2‑96 (MQ=255)
                                                                                     tGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAATTCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAAc          <  2:12898568/96‑1 (MQ=255)
                                                                                     tGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAAc          >  1:10439495/1‑96 (MQ=255)
                                                                                     tGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAAc          <  2:6459556/96‑1 (MQ=255)
                                                                                     tGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAAc          >  2:13336019/1‑96 (MQ=255)
                                                                                     tGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAAc          >  2:1906770/1‑96 (MQ=255)
                                                                                      gACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACg         >  1:13516421/1‑96 (MQ=255)
                                                                                       aCTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGc        <  1:6141244/96‑1 (MQ=255)
                                                                                        cTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCt       >  1:6723807/1‑96 (MQ=255)
                                                                                        cTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCt       >  1:7911462/1‑96 (MQ=255)
                                                                                          gggAAGTACAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTGCCATCAGATAGCCACGCTAGAGCAGAAAACGCTgg     <  1:4848668/96‑1 (MQ=255)
                                                                                          gggAAGTAAAATCCCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTgg     <  2:11079967/96‑1 (MQ=255)
                                                                                          gggAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTgg     >  2:6299816/1‑96 (MQ=255)
                                                                                               |                                                                                              
CAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATC  >  NC_000913/2820048‑2820237

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