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.8% 1980.2 / ‑2.8 436intergenic (+23/+19)b4542/oraAhypothetical protein/OraA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (249/185);  total (250/186)
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.11e-01

Only 100 of 468 total aligned reads displayed.

GCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGC  >  NC_000913/2820051‑2820231
                                                                                            |                                                                                        
gCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAg                                                                                          <  1:8223034/93‑1 (MQ=255)
  tAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAgta                                                                                        <  2:2184392/93‑3 (MQ=255)
  tAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAgta                                                                                        >  2:3942193/1‑91 (MQ=255)
   aaCGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAgtaa                                                                                       >  2:7233300/1‑90 (MQ=255)
    aCGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAgtaaa                                                                                      >  1:9465844/1‑89 (MQ=255)
    aCGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAgtaaa                                                                                      <  1:8760110/93‑5 (MQ=255)
     cGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTaaaa                                                                                    >  1:4248184/1‑93 (MQ=255)
      gAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAAt                                                                                   <  2:8129176/93‑1 (MQ=255)
        aaTATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATAc                                                                                 >  1:129989/1‑93 (MQ=255)
        aaTATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATAc                                                                                 <  1:8526302/93‑1 (MQ=255)
         atatCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATAcc                                                                                >  1:5268043/1‑93 (MQ=255)
          tatCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCg                                                                               <  2:7570549/93‑1 (MQ=255)
            tCGCCATCAGCTCCAGTGTGTCTTAACCGCCGGGCTGGTAGCTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTa                                                                             <  1:7670725/93‑1 (MQ=255)
            tCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTa                                                                             <  2:1537274/93‑1 (MQ=255)
            tCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTa                                                                             <  2:2638569/93‑1 (MQ=255)
            tCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTa                                                                             >  2:3445414/1‑93 (MQ=255)
            tCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTa                                                                             <  2:4642889/93‑1 (MQ=255)
             cGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTAt                                                                            <  1:3111521/93‑1 (MQ=255)
              gCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATg                                                                           <  2:5310531/93‑1 (MQ=255)
              gCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATg                                                                           <  2:1622839/93‑1 (MQ=255)
              gCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATg                                                                           <  2:9376997/93‑1 (MQ=255)
               ccaccaGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGc                                                                          <  1:4603285/93‑1 (MQ=255)
                 accaGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGt                                                                        >  2:3619089/1‑93 (MQ=255)
                 accaGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGt                                                                        >  1:1981627/1‑93 (MQ=255)
                 accaGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGt                                                                        >  2:5765765/1‑93 (MQ=255)
                  ccaGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGtt                                                                       >  1:6906812/1‑93 (MQ=255)
                  ccaGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGtt                                                                       >  1:1851073/1‑93 (MQ=255)
                  ccaGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGtt                                                                       >  1:1615696/1‑93 (MQ=255)
                    aGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCa                                                                     >  2:8770805/1‑93 (MQ=255)
                    aGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCa                                                                     <  2:1457898/93‑1 (MQ=255)
                    aGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCa                                                                     <  1:4397128/93‑1 (MQ=255)
                     gCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAg                                                                    >  1:4924542/1‑93 (MQ=255)
                     gCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAg                                                                    >  1:919715/1‑93 (MQ=255)
                      cTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGt                                                                   >  2:3353681/1‑93 (MQ=255)
                         cAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCgg                                                                >  1:3983092/1‑93 (MQ=255)
                         cAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCgg                                                                <  1:2335542/93‑1 (MQ=255)
                           gCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCa                                                              >  1:8006281/1‑93 (MQ=255)
                              tgtCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCaaaa                                                           >  2:3824002/1‑93 (MQ=255)
                                tCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAAtt                                                         >  2:1733071/1‑93 (MQ=255)
                                 cTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAAttt                                                        <  1:326647/93‑1 (MQ=255)
                                  ttAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTc                                                       >  2:7323855/1‑93 (MQ=255)
                                   tAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAACATTTCg                                                      >  2:6832243/1‑93 (MQ=255)
                                   tAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCg                                                      >  2:9481015/1‑93 (MQ=255)
                                    aaCCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGc                                                     >  2:9263120/1‑93 (MQ=255)
                                        gccgGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCaaa                                                 >  1:2651506/1‑93 (MQ=255)
                                        gccgGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCaaa                                                 >  2:3097641/1‑93 (MQ=255)
                                         ccgGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAt                                                <  2:6323781/93‑1 (MQ=255)
                                          cgGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAtc                                               <  2:7462034/93‑1 (MQ=255)
                                           gggCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAtct                                              >  1:10001465/1‑93 (MQ=255)
                                           gggCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGAATGCGTTCAGCCGGCAAAATTTCGCCAAAtcc                                              >  2:1038088/1‑92 (MQ=255)
                                            ggCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAtctc                                             >  2:6481856/1‑93 (MQ=255)
                                             gCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTcc                                            >  1:630987/1‑93 (MQ=255)
                                             gCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTcc                                            >  1:9651304/1‑93 (MQ=255)
                                             gCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTcc                                            <  2:8902638/93‑1 (MQ=255)
                                               tGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTg                                          <  1:5118360/93‑1 (MQ=255)
                                                ggTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGTCAAATCTCCTgg                                         >  2:1185784/1‑93 (MQ=255)
                                                 gTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGa                                        <  2:4743439/93‑1 (MQ=255)
                                                 gTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGa                                        <  1:7410808/93‑1 (MQ=255)
                                                  tAACTGAAACGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGat                                       >  2:7076349/1‑93 (MQ=255)
                                                  tAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGat                                       <  2:717228/93‑1 (MQ=255)
                                                  tAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGat                                       <  1:4524676/93‑1 (MQ=255)
                                                   aaCTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGata                                      <  1:9324787/93‑1 (MQ=255)
                                                   aaCTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGata                                      <  2:5079848/93‑1 (MQ=255)
                                                     cTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATc                                    >  1:6986306/1‑93 (MQ=255)
                                                     cTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATc                                    >  2:4102504/1‑93 (MQ=255)
                                                       gAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCtt                                  >  2:6635830/1‑93 (MQ=255)
                                                       gAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCtt                                  <  2:5812903/93‑1 (MQ=255)
                                                       gAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCtt                                  >  1:9089588/1‑93 (MQ=255)
                                                                gAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGAt                         >  1:143070/1‑93 (MQ=255)
                                                                gAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGAt                         >  2:9753676/1‑93 (MQ=255)
                                                                gAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGAt                         <  1:6614581/93‑1 (MQ=255)
                                                                gAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGAt                         <  1:3419666/93‑1 (MQ=255)
                                                                gAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGAt                         >  2:8511738/1‑93 (MQ=255)
                                                                 aaTAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATa                        <  2:6074227/93‑1 (MQ=255)
                                                                  ataagataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAg                       >  1:3730890/1‑93 (MQ=255)
                                                                   taagataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGc                      >  1:5397463/1‑93 (MQ=255)
                                                                    aagataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGcc                     >  1:4757009/1‑93 (MQ=255)
                                                                     agataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCa                    <  1:7921819/93‑1 (MQ=255)
                                                                     agataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCa                    <  2:1820980/93‑1 (MQ=255)
                                                                      gataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCACATCTCCTGGATATCTTCCATCAGATAGcccc                   >  1:1099643/1‑91 (MQ=255)
                                                                      gataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCAc                   >  1:2017736/1‑93 (MQ=255)
                                                                      gataagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCAc                   <  1:2106837/93‑1 (MQ=255)
                                                                        taagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGa                 >  2:7739877/1‑93 (MQ=255)
                                                                        taagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGa                 >  1:6110792/1‑93 (MQ=255)
                                                                        taagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGa                 <  1:6468968/93‑1 (MQ=255)
                                                                        taagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGa                 >  1:3002418/1‑93 (MQ=255)
                                                                         aagTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGAt                <  2:9891014/93‑1 (MQ=255)
                                                                           gTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATag              <  2:6925867/93‑1 (MQ=255)
                                                                             tttCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATagag            >  2:2069916/1‑93 (MQ=255)
                                                                             tttCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATagag            <  1:5685186/93‑1 (MQ=255)
                                                                              ttCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGc           >  2:8254142/1‑93 (MQ=255)
                                                                              ttCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGc           >  1:7130369/1‑93 (MQ=255)
                                                                               tCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCa          >  2:241058/1‑93 (MQ=255)
                                                                                cTTGACTGGGAAGTAAAATACCGTATGCGGTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGTTAGCCACGATAGAGCAg         >  1:8965081/1‑93 (MQ=255)
                                                                                  tGACTGGTAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGaa       <  1:4120517/93‑1 (MQ=255)
                                                                                  tGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCGGGATATCTTCCATCAGATAGCCACGATAGAGCAGaa       <  2:291118/93‑1 (MQ=255)
                                                                                    aCTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGaaaa     <  1:9951045/93‑1 (MQ=255)
                                                                                      tGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACg   >  1:5857203/1‑93 (MQ=255)
                                                                                      tGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACg   >  1:1674088/1‑93 (MQ=255)
                                                                                       gggAAGTAAAATACCGTATGCGTTAAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGc  >  2:4306071/1‑93 (MQ=255)
                                                                                            |                                                                                        
GCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGC  >  NC_000913/2820051‑2820231

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