Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 391,456 C→A 6.0% intergenic (+74/+192) SAS0343 → / ← SAS0344 30S ribosomal protein S18/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX571857391,4560CA 6.0% 970.8 / 18.3 245intergenic (+74/+192)SAS0343/SAS034430S ribosomal protein S18/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (101/127);  new base A (2/15);  total (103/142)
Fisher's exact test for biased strand distribution p-value = 9.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 247 total aligned reads displayed.

CATATGTTAAAGAAGAACAATAATATATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAAC  >  BX571857/391360‑391552
                                                                                                |                                                                                                
cATATGTTAAAGAAGAACAATAATATATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGc                                                                                                  >  1:2814059/1‑97 (MQ=255)
 atatGTTAAAGAAGAACAATAATATATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCa                                                                                                 >  1:501162/1‑97 (MQ=255)
  tatGTTAAAGAAGAACAATAATATATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCaa                                                                                                >  2:3272239/1‑97 (MQ=255)
  tatGTTAAAGAAGAACAATAATATATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCaa                                                                                                <  1:17831/97‑1 (MQ=255)
  tatGTTAAAGAAGAACAATAATATATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCaa                                                                                                >  1:731605/1‑97 (MQ=255)
       tAAAGAAGAACAATAATATATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAg                                                                                           <  1:2493530/97‑1 (MQ=255)
        aaaGAAGAACAATAATATATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAgg                                                                                          <  2:1966596/97‑1 (MQ=255)
         aagaagAACAATAATATATAATTTATTGTCAAACCCCGTAGGCATAGGCTTATGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGga                                                                                         >  1:789112/1‑97 (MQ=255)
          agaagaACAATAATATATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGgat                                                                                        <  2:601605/97‑1 (MQ=255)
           gaagaaCAATAATTTATAATTTCTTGTCAAACCTCGTAGGCATAGGCTTTCGGGGCTTTTTGTGTTTTGGGGTATTGATAAAGGGCAAAATAGgatg                                                                                       <  2:312073/97‑1 (MQ=255)
           gaagaaCAATAATATATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGgatg                                                                                       >  1:4645/1‑97 (MQ=255)
              gaaCAATAATATATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGatgatg                                                                                    >  1:2892234/1‑97 (MQ=255)
               aaCAATAATATATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGAtgt                                                                                   >  1:3261568/1‑97 (MQ=255)
                aCAATAATATATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGAtgtg                                                                                  <  1:1052125/97‑1 (MQ=255)
                 cAATAATATATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGa                                                                                 >  1:1586896/1‑97 (MQ=255)
                  aataatATATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGaa                                                                                >  1:2578883/1‑97 (MQ=255)
                  aataatATATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTGGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGaa                                                                                <  2:3234898/97‑1 (MQ=255)
                   ataataTATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAAt                                                                               <  1:415832/97‑1 (MQ=255)
                   ataataTATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAAt                                                                               <  1:3478296/97‑1 (MQ=255)
                   ataataTATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAAt                                                                               <  1:2609227/97‑1 (MQ=255)
                     aataTATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAACGGGCAAAAAAGGATGATGTGAATGt                                                                             <  1:2700421/97‑1 (MQ=255)
                     aaaatataATTTATTGTCAAACCCCGTAGGCATAAGCTTACGGGGCTTTTTGTGTTTTGGGGTATAAAAAAAAAGAAAAAAAGGATGATGTGAATGt                                                                             <  2:1612959/94‑1 (MQ=255)
                      atatatAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGtt                                                                            <  1:2848692/97‑1 (MQ=255)
                       tatataATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGttt                                                                           <  1:2105526/97‑1 (MQ=255)
                        atataATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGtttt                                                                          <  1:2390062/97‑1 (MQ=255)
                        atataATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGtttt                                                                          <  1:2295771/97‑1 (MQ=255)
                        atagaATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGAAAAAAAAAGGATGATGTGAATGtttt                                                                          <  2:140068/97‑1 (MQ=255)
                         tataATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGCGTTTTGGGGTATAGAAACAGGGCAAAAAAGGATGATGTGCATGTTTtg                                                                         >  1:1599559/1‑97 (MQ=255)
                          ataATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTtgt                                                                        <  2:2132809/97‑1 (MQ=255)
                             aTTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGCAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGtt                                                                     >  2:116291/1‑97 (MQ=255)
                               ttATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCg                                                                   >  1:2054737/1‑97 (MQ=255)
                                  ttGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAAAGAAAAAAAGGATGATGTGAATGTTTTGTGTTCGGaa                                                                <  2:1781019/97‑1 (MQ=255)
                                  ttGTCAAACCCCGTAGGCATAGGCTTAAGGGGCTTTTTGTGTTTTGGGGTATAAAAAAAAAAAAAAAAAGGATGATGTGAATGTTTTGTGTTCGGaa                                                                <  2:1480164/97‑1 (MQ=255)
                                   tggcAAACCCCGTAGGCGTAGGCTTACGGGGCTTTTGGTGGTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAAt                                                               <  2:3383057/94‑1 (MQ=255)
                                   cgccAACCCCCGTCGGCATAGGCGTACGTGGCGTTTTGTGGTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAAt                                                               <  2:2225755/94‑1 (MQ=255)
                                    gTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAAtt                                                              >  2:1954572/1‑97 (MQ=255)
                                    gTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAAtt                                                              >  1:529249/1‑97 (MQ=255)
                                     tCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAAttt                                                             >  2:1108867/1‑97 (MQ=255)
                                      cAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTg                                                            >  1:3342640/1‑97 (MQ=255)
                                       aaaCCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGc                                                           <  2:2092600/97‑1 (MQ=255)
                                       aaaCCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGc                                                           >  2:1671385/1‑97 (MQ=255)
                                       aaaCCCCGTAGGAAAAGGCTGACGGGGATGTTTGTGTTTTGGGGTATAGAAAAAGAAAAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGc                                                           <  2:2517592/97‑1 (MQ=255)
                                        aaCCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGca                                                          <  1:2420619/97‑1 (MQ=255)
                                         aCCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGcac                                                         >  1:3189132/1‑97 (MQ=255)
                                          ccccGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGcaca                                                        <  1:632069/97‑1 (MQ=255)
                                          ccccGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGcaca                                                        <  1:192893/97‑1 (MQ=255)
                                          ccccGTAGGCATAGGCTTACGGGGCTTTTTGTGTGTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGcaca                                                        >  1:3207431/1‑97 (MQ=255)
                                             cGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAg                                                     <  2:1841229/97‑1 (MQ=255)
                                                aGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGata                                                  <  2:2358469/97‑1 (MQ=255)
                                                aGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGata                                                  <  1:2013471/97‑1 (MQ=255)
                                                aGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGata                                                  >  2:3012755/1‑97 (MQ=255)
                                                 ggCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGatat                                                 >  1:817003/1‑97 (MQ=255)
                                                 ggCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGatat                                                 <  1:303058/97‑1 (MQ=255)
                                                 ggCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGatat                                                 <  1:2198003/97‑1 (MQ=255)
                                                 ggCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGatat                                                 <  1:2224324/97‑1 (MQ=255)
                                                 ggCATAGGATTACGGGGCTTTTTGTGTTTTGGTGGATAGAAAAAGAGAAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGatat                                                 <  2:25820/97‑1 (MQ=255)
                                                  gCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATg                                                <  1:995115/97‑1 (MQ=255)
                                                   cATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGt                                               >  1:2413896/1‑97 (MQ=255)
                                                    aTAGACTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAAAGAAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGtt                                              <  2:1338262/97‑1 (MQ=255)
                                                     tAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGttt                                             >  1:557343/1‑97 (MQ=255)
                                                     tAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGttt                                             >  1:1682100/1‑97 (MQ=255)
                                                     tAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGttt                                             >  2:2881298/1‑97 (MQ=255)
                                                      aGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTta                                            >  2:2272749/1‑97 (MQ=255)
                                                       ggCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTtat                                           >  2:217321/1‑97 (MQ=255)
                                                       ggCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTtat                                           <  1:2516290/97‑1 (MQ=255)
                                                         cTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTatat                                         >  1:1343952/1‑97 (MQ=255)
                                                         cTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTatat                                         <  1:846083/97‑1 (MQ=255)
                                                               gggCTTTTTTTGTTTTAGGGTATAGAAAAACGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCaaa                                   <  2:724533/97‑1 (MQ=255)
                                                                ggCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCaaaa                                  <  1:975916/97‑1 (MQ=255)
                                                                ggCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCaaaa                                  <  1:1445410/97‑1 (MQ=255)
                                                                 gCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCaaaaa                                 >  1:1644649/1‑97 (MQ=255)
                                                                    ttttGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAaataa                              <  2:3596407/97‑1 (MQ=255)
                                                                    ttttGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAaataa                              <  2:2088329/97‑1 (MQ=255)
                                                                    ttttGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAaataa                              <  1:1257893/97‑1 (MQ=255)
                                                                    ttttGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAaataa                              <  1:207972/97‑1 (MQ=255)
                                                                    ttttGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAaataa                              <  1:2932726/97‑1 (MQ=255)
                                                                    ttttGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAaataa                              >  2:2896535/1‑97 (MQ=255)
                                                                     tttGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAaataat                             <  1:3510021/97‑1 (MQ=255)
                                                                      ttGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTAAAAATGCAAAAataata                            <  2:3154297/97‑1 (MQ=255)
                                                                       tgtgTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAatat                           <  1:3667979/97‑1 (MQ=255)
                                                                          gtTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGaa                        <  2:534138/97‑1 (MQ=255)
                                                                          gtTTTGGGGTATAGAAAAAGGGAAAAAAAGTATGATGTGAATGTTTTGTGGTTGGAATTTGCACAAAGATATGTTTATATTGCTAAAATAATATGaa                        >  2:3567618/1‑97 (MQ=255)
                                                                            tttGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAAtt                      <  1:2755033/97‑1 (MQ=255)
                                                                              tGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTa                    <  2:959237/97‑1 (MQ=255)
                                                                               ggggTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAg                   <  2:2341901/97‑1 (MQ=255)
                                                                                gggTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTGGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGa                  <  2:563983/97‑1 (MQ=255)
                                                                                 ggTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGAt                 <  1:746798/97‑1 (MQ=255)
                                                                                  gTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATg                >  2:3687188/1‑97 (MQ=255)
                                                                                   tataGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGc               <  2:899854/97‑1 (MQ=255)
                                                                                    ataGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCa              >  2:3501884/1‑97 (MQ=255)
                                                                                    ataGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCa              <  2:1304693/97‑1 (MQ=255)
                                                                                       gAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATaa           >  1:358443/1‑97 (MQ=255)
                                                                                       gAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATaa           <  1:1061775/97‑1 (MQ=255)
                                                                                        aaaaaGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATaaa          <  2:3061921/97‑1 (MQ=255)
                                                                                          aaaGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATaaaaa        <  1:1720928/97‑1 (MQ=255)
                                                                                           aaGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATaaaaaa       >  2:315865/1‑97 (MQ=255)
                                                                                            aGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATaaaaaaa      >  2:3407843/1‑97 (MQ=255)
                                                                                              ggCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGa    >  1:1111332/1‑97 (MQ=255)
                                                                                               gCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGaa   >  2:1560834/1‑97 (MQ=255)
                                                                                               gCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGaa   >  1:1903578/1‑97 (MQ=255)
                                                                                                |                                                                                                
CATATGTTAAAGAAGAACAATAATATATAATTTATTGTCAAACCCCGTAGGCATAGGCTTACGGGGCTTTTTGTGTTTTGGGGTATAGAAAAAGGGCAAAAAAGGATGATGTGAATGTTTTGTGTTCGGAATTTGCACAAAGATATGTTTATATTGCAAAAATAATATGAATTTAGATGCATAAAAAAAGAAC  >  BX571857/391360‑391552

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: