Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 431,505 T→C 7.4% S77P (TCA→CCA)  SAS0387 → exotoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX571857431,5050TC 7.4% 730.3 / 5.6 203S77P (TCA→CCA) SAS0387exotoxin
Reads supporting (aligned to +/- strand):  ref base T (73/109);  new base C (7/12);  total (80/123)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 210 total aligned reads displayed.

AGCAAATGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGACACATC  >  BX571857/431408‑431602
                                                                                                 |                                                                                                 
aGCAAATGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAAt                                                                                                   >  1:2735388/1‑98 (MQ=255)
aGCAAATGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAAt                                                                                                   <  2:1415455/98‑1 (MQ=255)
aGCAAATGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAAt                                                                                                   <  2:494610/98‑1 (MQ=255)
 gCAAATGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATc                                                                                                  >  1:676516/1‑98 (MQ=255)
  cAAATGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCa                                                                                                 >  1:2946072/1‑98 (MQ=255)
  cAAATGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCa                                                                                                 <  2:921737/98‑1 (MQ=255)
  cAAATGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCa                                                                                                 <  1:3289576/98‑1 (MQ=255)
   aaaTGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCaa                                                                                                <  2:2227533/98‑1 (MQ=255)
   aaaTGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCaa                                                                                                <  1:1767119/98‑1 (MQ=255)
     aTGCAACAACACCATCTTCAACTAAAGTGTAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCaaaa                                                                                              <  2:2029039/98‑1 (MQ=255)
     aTGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGACGCACCGCAACCaaaa                                                                                              >  2:1575415/1‑98 (MQ=255)
     aTGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCaaaa                                                                                              <  2:1847272/98‑1 (MQ=255)
     aTGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCaaaa                                                                                              >  1:1550352/1‑98 (MQ=255)
       gCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCacaaccacaacc                                                                                            >  2:474345/1‑98 (MQ=255)
       gCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAAcc                                                                                            <  1:1719289/98‑1 (MQ=255)
       gCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAAcc                                                                                            <  2:1458661/98‑1 (MQ=255)
        caacaaCACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGACGCACCGCAACCACAAccc                                                                                           >  2:2263376/1‑97 (MQ=17)
        caacaaCACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCa                                                                                           <  2:2989848/98‑1 (MQ=255)
         aacaacACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCaa                                                                                          >  1:1581751/1‑98 (MQ=255)
         aacaacACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCaa                                                                                          <  2:2564951/98‑1 (MQ=255)
         aacaacACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCaa                                                                                          <  1:531721/98‑1 (MQ=255)
         aacaacACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCaa                                                                                          <  2:2697431/98‑1 (MQ=255)
          acaacaCCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCaaa                                                                                         <  1:3085818/98‑1 (MQ=255)
          acaacaCCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCaaa                                                                                         <  1:2836891/98‑1 (MQ=255)
             acacCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAAcgc                                                                                      >  1:1572941/1‑98 (MQ=255)
             acacCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAAcgc                                                                                      >  1:2994300/1‑98 (MQ=255)
              cacCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAAcgcg                                                                                     <  2:112420/98‑1 (MQ=255)
              cacCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAAcgcg                                                                                     <  2:2792975/98‑1 (MQ=255)
               acCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAATTCCATACGCGa                                                                                    <  2:1912024/98‑1 (MQ=255)
                 cATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGaca                                                                                  <  1:2378216/98‑1 (MQ=255)
                 cATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGaca                                                                                  <  2:1818109/98‑1 (MQ=255)
                  ctCTTCAGCTAAGGTGCAAGCACCCCAACAAGCATCAAACGCGTCAACACCATCTTCACCTAATGTAGAAGCGCCGCCATCAAAACCAAACGCGtcaa                                                                                 <  2:1392244/97‑1 (MQ=255)
                   tCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGacaac                                                                                >  1:1333065/1‑98 (MQ=255)
                   tCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGacaac                                                                                >  1:2315086/1‑98 (MQ=255)
                    cTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGacaaca                                                                               >  2:3077136/1‑98 (MQ=255)
                     ttCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAacac                                                                              <  1:1995288/98‑1 (MQ=255)
                      tCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACAcc                                                                             <  1:3179069/98‑1 (MQ=255)
                        aaCTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCAt                                                                           <  1:3018960/98‑1 (MQ=255)
                         aCTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATc                                                                          <  1:2918331/98‑1 (MQ=255)
                         aCTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATc                                                                          <  2:2481165/98‑1 (MQ=255)
                          cTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCt                                                                         <  1:598800/98‑1 (MQ=255)
                           tAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCtt                                                                        <  1:1300282/98‑1 (MQ=255)
                             aaGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCa                                                                      <  2:1317470/98‑1 (MQ=255)
                                tGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACt                                                                   <  2:204372/98‑1 (MQ=255)
                                 ggAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTa                                                                  >  1:363873/1‑98 (MQ=255)
                                 ggAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTa                                                                  >  1:776237/1‑98 (MQ=255)
                                  gAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTaa                                                                 <  2:839130/98‑1 (MQ=255)
                                    aGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACACGACAACACCATCTTCAACTAAAg                                                               >  2:119553/1‑98 (MQ=255)
                                         cGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGaa                                                          >  2:724609/1‑98 (MQ=255)
                                           caacaaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGc                                                        <  1:411160/98‑1 (MQ=255)
                                           caacaaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGc                                                        <  1:259836/98‑1 (MQ=255)
                                              caaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCAcc                                                     >  2:3286666/1‑98 (MQ=255)
                                                agcagcAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGc                                                   >  2:2598841/1‑98 (MQ=255)
                                                 gcagcaAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGca                                                  <  1:2490154/98‑1 (MQ=255)
                                                 gcagcaAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGca                                                  >  1:1317530/1‑98 (MQ=255)
                                                 gcagcaAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGca                                                  >  2:1683882/1‑98 (MQ=255)
                                                      aaaCGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAag                                             >  2:1897208/1‑98 (MQ=255)
                                                      aaaCGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAag                                             <  1:1000407/98‑1 (MQ=255)
                                                      aaaCGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAag                                             <  2:2178101/98‑1 (MQ=255)
                                                      aaaCGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAag                                             >  1:1276525/1‑98 (MQ=255)
                                                       aaCGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAagc                                            <  2:1650147/98‑1 (MQ=255)
                                                        aCGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAagca                                           <  1:632392/98‑1 (MQ=255)
                                                        aCGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAagca                                           >  1:462029/1‑98 (MQ=255)
                                                        aCGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAagca                                           >  1:219343/1‑98 (MQ=255)
                                                        aCGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAagca                                           >  1:2578076/1‑98 (MQ=255)
                                                              caacaCCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAAc                                     >  2:3151482/1‑98 (MQ=255)
                                                              caacaCCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAAc                                     >  2:1984637/1‑98 (MQ=255)
                                                              caacaCCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAAc                                     >  2:235600/1‑98 (MQ=255)
                                                                  acCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGa                                 <  1:145987/98‑1 (MQ=255)
                                                                   ccATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGac                                >  1:2012014/1‑98 (MQ=255)
                                                                    caccTTCAACTAAAGTAGCAGCACCGCAATCAGAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGaca                               <  2:175799/95‑1 (MQ=255)
                                                                      tCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGacaac                             <  1:2820417/98‑1 (MQ=255)
                                                                       cTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGacaaca                            >  2:428364/1‑98 (MQ=255)
                                                                       cTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGacaaca                            <  1:18483/98‑1 (MQ=255)
                                                                         tCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAAcacc                          <  1:1237210/98‑1 (MQ=255)
                                                                         tCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAAcacc                          >  1:89253/1‑98 (MQ=255)
                                                                          cAACTAAAGTAGAAGCACCGCAATCAGAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAAcacca                         <  1:21126/98‑1 (MQ=255)
                                                                           aaCTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAAcaccac                        >  2:1457323/1‑98 (MQ=255)
                                                                            aCTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACaccacc                       <  1:3130002/98‑1 (MQ=255)
                                                                            aCTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACaccacc                       <  1:1883276/98‑1 (MQ=255)
                                                                             cTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCt                      >  2:205717/1‑98 (MQ=255)
                                                                             cTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCt                      >  2:459528/1‑98 (MQ=255)
                                                                               aaaGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATTTTCAACGAAACTGGAAGCCCCGCCACAAGCAGCTAACGCGACCACACCACCtta                    >  1:1380446/1‑97 (MQ=25)
                                                                               aaaGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTc                    <  1:2148881/98‑1 (MQ=255)
                                                                                 aGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGt                  >  2:1715610/1‑98 (MQ=255)
                                                                                 aGTAGAAACACCGCAACCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGt                  <  2:2278991/98‑1 (MQ=255)
                                                                                  gTAGAAGCACCGCAACCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTc                 <  2:1877785/98‑1 (MQ=255)
                                                                                    aGAAGCACCGCAACCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTa               <  2:2997041/98‑1 (MQ=255)
                                                                                      aaGCACCGCAATCAAAACCAAACGCGACCACACCCTCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAAt             <  2:1627646/98‑1 (MQ=255)
                                                                                       aGCACCGCAACAAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATg            <  2:689638/98‑1 (MQ=255)
                                                                                        gCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGt           >  1:1008053/1‑98 (MQ=255)
                                                                                        gCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGt           >  2:364174/1‑98 (MQ=255)
                                                                                         cACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTa          <  1:2877619/98‑1 (MQ=255)
                                                                                         cACCGCAATCAAAACCAAACGCGACAACACCAGCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTa          <  2:82991/98‑1 (MQ=255)
                                                                                          aCCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAg         <  1:647449/98‑1 (MQ=255)
                                                                                           ccGCAACCAGAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGa        <  2:2622283/98‑1 (MQ=255)
                                                                                            cGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGac       >  1:997518/1‑98 (MQ=255)
                                                                                             gCAATTAGACCCAAACGCGACATCACCATCTTCAACTCAAGTGGAAGCACCGCAAAAAGCAGCGAACGAGACAGCACCACCTTCGTCGAACGTAGaca      <  2:2109594/98‑1 (MQ=17)
                                                                                              cAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGacac     <  1:732296/98‑1 (MQ=255)
                                                                                                 tCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGACACATc  <  1:1333703/98‑1 (MQ=255)
                                                                                                 |                                                                                                 
AGCAAATGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGACACATC  >  BX571857/431408‑431602

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