Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 318,690 A→G 54.2% S438P (TCA→CCA)  rclA ← reactive chlorine stress species (RCS) resistance protein; pyridine nucleotide‑dependent disulfide oxidoreductase family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913318,6900AG54.2% 49.4 / 212.8 168S438P (TCA→CCA) rclAreactive chlorine stress species (RCS) resistance protein; pyridine nucleotide‑dependent disulfide oxidoreductase family
Reads supporting (aligned to +/- strand):  ref base A (40/37);  new base G (47/44);  total (87/81)
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.73e-01

Only 100 of 177 total aligned reads displayed.

TAAGAAGTTATAGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCGGCATCCATCACCATTTTCACT  >  NC_000913/318592‑318787
                                                                                                  |                                                                                                 
tAAGAAGTTATAGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGaa                                                                                                  <  1:1386455/100‑1 (MQ=255)
 aagaagTTATAGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGgaa                                                                                                 <  1:2465352/100‑1 (MQ=255)
 aagaagTTATAGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGaaa                                                                                                 >  1:897258/1‑100 (MQ=255)
 aagaagTTATAGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTATGGCGTCTGATTTTGAGTTTATTTGACTAATGgaa                                                                                                 >  1:1753899/1‑100 (MQ=255)
  agaagTTATAGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCAAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGaaaa                                                                                                >  2:505030/1‑100 (MQ=255)
  agaagTTATAGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGACAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGgaaa                                                                                                >  2:1584284/1‑100 (MQ=255)
   tctgTTATAGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAAt                                                                                               <  1:3386499/97‑1 (MQ=255)
   gaagTTATAGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAAt                                                                                               <  2:460261/100‑1 (MQ=255)
   gaagTTATAGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAAt                                                                                               <  2:1502764/100‑1 (MQ=255)
    aagTTATAGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATa                                                                                              <  1:135653/100‑1 (MQ=255)
      gTTATAGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGa                                                                                            <  2:1961694/100‑1 (MQ=255)
      gTTATAGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGACTTTGAGTTTATTTGACTAATGGAAATAGa                                                                                            >  1:763952/1‑100 (MQ=255)
      gTTATAGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGAGTCTGATTTTGAGTTTATTTGACTAATGAAAATAg                                                                                             <  2:3961012/99‑1 (MQ=255)
       ttATAGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGa                                                                                            >  2:3415607/1‑99 (MQ=255)
       ttATAGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGAt                                                                                           >  1:916086/1‑100 (MQ=255)
       ttATAGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGAt                                                                                           <  2:3062976/100‑1 (MQ=255)
        tataGAGTAACAGAGAATTAATAAGATTCTTTTTTCTGTGACTCCTGAATATTTTTTCTGGAGTCTGATTTTGAGTTTATTTTACTAATGGAAATAGATc                                                                                          <  2:696574/100‑1 (MQ=255)
        tataGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATc                                                                                          <  2:1765590/100‑1 (MQ=255)
        tataGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATc                                                                                          <  2:834243/100‑1 (MQ=255)
        tataGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATc                                                                                          >  2:253138/1‑100 (MQ=255)
          taGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCan                                                                                        <  2:3155050/100‑2 (MQ=255)
           agagTAACAGGGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCAtt                                                                                       >  1:609222/1‑100 (MQ=255)
           agagTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCAtt                                                                                       >  2:3391443/1‑100 (MQ=255)
            gagTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTg                                                                                      >  2:2104136/1‑100 (MQ=255)
              gTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAg                                                                                    >  2:1114976/1‑100 (MQ=255)
              gTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAg                                                                                    <  1:819878/100‑1 (MQ=255)
                  cAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGAt                                                                                >  2:609122/1‑100 (MQ=255)
                  cAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGAt                                                                                >  1:2676164/1‑100 (MQ=255)
                   agagAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGAtt                                                                               <  1:1077880/100‑1 (MQ=255)
                    gagaATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTTTGATTTTGGGTTTATTTGACTAATGAAGATAGATCATTGAGTGATTc                                                                              <  2:1501044/100‑1 (MQ=255)
                    gagaATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGa                                                                                 >  2:2593252/1‑97 (MQ=255)
                    gagaATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTc                                                                              >  2:1338892/1‑100 (MQ=255)
                      gaATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGc                                                                            <  1:2143178/100‑1 (MQ=255)
                      gaATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGc                                                                            >  1:2264340/1‑100 (MQ=255)
                            ataaGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCg                                                                      <  2:2926017/100‑1 (MQ=255)
                             taaGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGCGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCg                                                                      >  2:3132301/1‑99 (MQ=255)
                             taaGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGa                                                                     >  1:1553490/1‑100 (MQ=255)
                                gATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACgg                                                                  >  1:3380860/1‑100 (MQ=255)
                                gATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACgg                                                                  >  1:3625612/1‑100 (MQ=255)
                                  ttCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGAt                                                                >  2:36488/1‑100 (MQ=255)
                                   tCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATg                                                               <  1:2543072/100‑1 (MQ=255)
                                   tCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATg                                                               <  2:2949613/100‑1 (MQ=255)
                                    cTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGa                                                              >  1:170443/1‑100 (MQ=255)
                                    cTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGa                                                              >  2:1948424/1‑100 (MQ=255)
                                     ttttCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAg                                                             <  2:2281719/100‑1 (MQ=255)
                                     ctcaaaatcagaCGCCAGAACAAATATTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAg                                                             >  2:1024460/10‑100 (MQ=21)
                                      tcaaaatcagaCGCCAGAACAAATATTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGt                                                            >  1:5026/9‑100 (MQ=21)
                                       ttCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTa                                                           <  1:2028923/100‑1 (MQ=255)
                                         ctctGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGAACATTGAGTGATTCGCTCATCGACGGATGAGTaaa                                                         >  1:2422504/1‑100 (MQ=255)
                                         ctctGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTaaa                                                         >  2:742962/1‑100 (MQ=255)
                                         aaatcagaCGCCAGAACAAATATTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTaaa                                                         >  2:3567564/6‑100 (MQ=25)
                                           atcagccGCCAGAACAAATATTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAata                                                       >  2:2819329/7‑100 (MQ=37)
                                           atcagaCGCCAGAACAAATATTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAata                                                       >  1:3312131/4‑100 (MQ=37)
                                            tGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGCCGGATGAGTAAatat                                                      >  1:4068765/1‑100 (MQ=255)
                                             cagaCGCCAGAACAAATATTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATc                                                     >  1:2438297/2‑100 (MQ=37)
                                              agaCGCCAGAATATTTGTTCTAGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCt                                                    <  2:3796075/100‑1 (MQ=255)
                                              agaCGCCAGAACAAATATTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCCTTGAGTGATTCGCTCATCGACGGATGCGTAAATATCt                                                    >  2:3496210/1‑100 (MQ=38)
                                                 cGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGAt                                                 >  1:1171558/1‑100 (MQ=255)
                                                     aGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATcgcg                                             <  1:1962875/100‑1 (MQ=255)
                                                       aaTATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTa                                           <  1:58764/100‑1 (MQ=255)
                                                        atatTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTaa                                          <  2:3370776/100‑1 (MQ=255)
                                                          atTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAata                                        >  2:2235581/1‑100 (MQ=255)
                                                             tGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGCTCGCGTAatatcc                                     <  1:2800606/100‑3 (MQ=255)
                                                             tGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGc                                     <  2:3565895/100‑1 (MQ=255)
                                                             tGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGc                                     >  2:2262013/1‑100 (MQ=255)
                                                             tGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGc                                     <  2:327469/100‑1 (MQ=255)
                                                              gTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCt                                    <  2:1092942/100‑1 (MQ=255)
                                                                 cTGGCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCtat                                  <  2:2700707/99‑1 (MQ=255)
                                                                   ggCGTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAg                               <  1:2992702/100‑1 (MQ=255)
                                                                    gCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTTGCTCTTCGACGGATGAGGAAATATCTGATCGTGTAATATGCTATAAgg                              >  2:2591540/1‑100 (MQ=255)
                                                                    gCGTCTGATTTTGAGTATGTTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGAGGAATGAGTAAATATCTGATCGCGTAAGATGCTATAAgg                              <  2:988162/100‑1 (MQ=255)
                                                                     cGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGc                             <  1:1109176/100‑1 (MQ=255)
                                                                      gTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCa                            >  1:1951054/1‑100 (MQ=255)
                                                                      gTCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCa                            <  1:609128/100‑1 (MQ=255)
                                                                      gTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCa                            >  1:2252430/1‑100 (MQ=255)
                                                                       tCTGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAg                           <  1:1794118/100‑1 (MQ=255)
                                                                        cTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCtagatcggaag                          >  2:3391523/1‑91 (MQ=255)
                                                                        cTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATAGCGTAATAAGCTATAAGGc                             >  2:2460365/1‑97 (MQ=255)
                                                                         tGATTTTGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGcc                         >  2:2916370/1‑100 (MQ=255)
                                                                            ttttGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCgg                      >  1:699937/1‑100 (MQ=255)
                                                                              ttGAGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCGGca                    <  1:3900821/100‑1 (MQ=255)
                                                                                gAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCGGcatc                  <  1:1260351/100‑1 (MQ=255)
                                                                                gAGTTTAGTTGACTAATGGAACTAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTTATCGCGTAATAGTCTATAAGGCAGCCCGGcatc                  <  2:422061/100‑1 (MQ=255)
                                                                                 aGTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCGGcatcc                 >  1:2409640/1‑100 (MQ=255)
                                                                                  gTTTATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCGGcatcca                >  1:1507043/1‑100 (MQ=255)
                                                                                  gTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCGGcatcca                >  1:3433883/1‑100 (MQ=255)
                                                                                  gTTTATGTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCGGcatcca                <  1:753020/100‑1 (MQ=255)
                                                                                   tttATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCGGcatccat               >  1:2379378/1‑100 (MQ=255)
                                                                                    ttATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCGGcatccatc              <  1:3958525/100‑1 (MQ=255)
                                                                                    ttATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCGGcatccatc              <  1:313063/100‑1 (MQ=255)
                                                                                    ttATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCGGcatccatc              >  2:909427/1‑100 (MQ=255)
                                                                                     tATTTGACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCGGCATCCATCa             <  1:1467787/100‑1 (MQ=255)
                                                                                          gACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCGGCATCCATCACCAtt        >  1:2495007/1‑100 (MQ=255)
                                                                                          gACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCGGCATCCATCACCAtt        >  1:1604357/1‑100 (MQ=255)
                                                                                          gACTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGAGTAATATGCTATAAGGCAGCCCGGCATCCATCACCAtt        >  2:2249235/1‑100 (MQ=255)
                                                                                          gACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATCAGGCAGCCCGGCATCCATCACCAtt        >  2:2671688/1‑100 (MQ=255)
                                                                                            cTAATGGAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCGGCATCCATCACCAtttt      >  1:704232/1‑100 (MQ=255)
                                                                                               aTGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCGGCATCCATCACCATTTTCAc   <  1:2748593/100‑1 (MQ=255)
                                                                                                tGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCGGCATCCATCACCATTTTCAc   <  1:3125526/99‑1 (MQ=255)
                                                                                                tGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCGGCATCCATCACCATTTTCACt  >  1:3360805/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
TAAGAAGTTATAGAGTAACAGAGAATTAATAAGATTCTTTTCTCTGAGACGCCAGAATATTTGTTCTGGCGTCTGATTTTGAGTTTATTTGACTAATGAAAATAGATCATTGAGTGATTCGCTCATCGACGGATGAGTAAATATCTGATCGCGTAATATGCTATAAGGCAGCCCGGCATCCATCACCATTTTCACT  >  NC_000913/318592‑318787

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