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% 901.1 / 6.8 254S77P (TCA→CCA) SAS0387exotoxin
Reads supporting (aligned to +/- strand):  ref base T (104/120);  new base C (13/11);  total (122/132)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 259 total aligned reads displayed.

AGCAAATGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGACACATC  >  BX571857/431408‑431602
                                                                                                 |                                                                                                 
aGCAAATGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAAt                                                                                                   <  1:890722/98‑1 (MQ=255)
 gCAAATGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATc                                                                                                  <  1:3169394/98‑1 (MQ=255)
  cAAATGCAACAACACCATCTACAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAACACCGCAACCa                                                                                                 >  2:3269763/1‑98 (MQ=255)
    aaTGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCaaa                                                                                               <  2:3418120/98‑1 (MQ=255)
     acgCAACAACACCATCTTCAACTAAAGTGGCAGCACCGCCACAAGCAGCAAACGCGACAACACCACCTTCACCTACAGTAGACGCCCCGCCACCaaaa                                                                                              >  2:1813898/3‑98 (MQ=11)
     aTGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCGACTAAAGTAGACGCAccacaaccacaa                                                                                              >  2:2532579/1‑98 (MQ=14)
     aTGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCaaaa                                                                                              <  2:1989428/98‑1 (MQ=255)
         aacaacACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCaa                                                                                          <  1:534623/98‑1 (MQ=255)
         aacaacACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAACCAAAACCac                                                                                          >  2:2990824/1‑97 (MQ=255)
          acaacaCCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCaaa                                                                                         <  1:626069/98‑1 (MQ=255)
           caacaCCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAAc                                                                                        >  1:1568903/1‑98 (MQ=255)
           caacaCCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAACCAAAACCAAAc                                                                                        >  2:2013635/1‑98 (MQ=255)
             acacCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCCACCAAAACCAAAcgc                                                                                      >  2:2178231/1‑98 (MQ=255)
              ctccATCTTCAACTAAAGTGGAGGCACCGCCACAAGCAGCAAACGCGACTACACCTGCTTCAACTAAAGGAGAAGCACCGCAATCAAAACCGAAcgcg                                                                                     <  2:2945880/96‑1 (MQ=21)
              cacCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAAcgcg                                                                                     <  1:2160436/98‑1 (MQ=255)
              cacCATCTCCAACTGTAGTGGAAGCACCGCAACAACCCGCAACCGCGACAACACCATCTTCAACTAGCGTAGAAGCTCCGCAATCAATCCCAAAcgcg                                                                                     <  2:1669365/98‑1 (MQ=18)
               acCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGa                                                                                    <  2:683835/98‑1 (MQ=255)
                ccATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGac                                                                                   >  1:1471263/1‑98 (MQ=255)
                 cATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCCACCAAAACCAAACGCGaca                                                                                  >  2:948546/1‑98 (MQ=255)
                  aTCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGacaa                                                                                 >  1:2117099/1‑98 (MQ=255)
                  aTCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGacaa                                                                                 >  1:3595833/1‑98 (MQ=255)
                   tCTTCAACTAAAGTGGAACCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCAACTAAAGTCGAAGCACCACAAGCAAAACCAACCACaacaac                                                                                >  2:1018619/1‑98 (MQ=17)
                     ttCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAacac                                                                              <  2:2332732/98‑1 (MQ=255)
                       cAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCa                                                                            >  2:702489/1‑98 (MQ=255)
                         aCTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATc                                                                          <  2:1755300/98‑1 (MQ=255)
                          cTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGCCAACACCATCTTCAAATAAACAAGAAGCACCGCAATCAAAACCAAACGAGACAACACCATCt                                                                         <  1:1116/98‑1 (MQ=21)
                          cTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCt                                                                         <  1:3246515/98‑1 (MQ=255)
                          cTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCt                                                                         <  2:1131116/98‑1 (MQ=255)
                           tAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCtt                                                                        >  2:1178889/1‑98 (MQ=255)
                             aaGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCa                                                                      <  2:795703/98‑1 (MQ=255)
                              aGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCaa                                                                     >  1:1668403/1‑98 (MQ=255)
                              aGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAAAACCATCTTCaa                                                                     >  1:1692200/1‑98 (MQ=255)
                                tGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACt                                                                   >  1:2893526/1‑98 (MQ=255)
                                   aaGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTaaa                                                                >  2:2635225/1‑98 (MQ=255)
                                   aaGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTaaa                                                                >  2:3244853/1‑98 (MQ=255)
                                   aaGCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTaaa                                                                >  1:2186001/1‑98 (MQ=255)
                                     gCACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGt                                                              <  1:1999054/98‑1 (MQ=255)
                                      cACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTg                                                             >  1:1740696/1‑98 (MQ=255)
                                      cACCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTg                                                             >  1:123883/1‑98 (MQ=255)
                                       aCCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTgg                                                            <  2:1480155/98‑1 (MQ=255)
                                       aCCGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTgg                                                            <  1:3344747/98‑1 (MQ=255)
                                        ccGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGa                                                           <  2:1545438/98‑1 (MQ=255)
                                         cGCAACAAGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGaa                                                          >  1:1879088/1‑98 (MQ=255)
                                           caacaaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGc                                                        <  2:592197/98‑1 (MQ=255)
                                            aacaaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCa                                                       <  2:1577953/98‑1 (MQ=255)
                                              caaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCAcc                                                     >  1:176822/1‑98 (MQ=255)
                                              caaGCAGCAAACGCGACAACACCACCTTCAACTAAAGGAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCAcc                                                     <  2:2233132/98‑1 (MQ=255)
                                               aaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCg                                                    >  2:1236394/1‑98 (MQ=255)
                                               aaGCAGCAAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCg                                                    >  1:506992/1‑98 (MQ=255)
                                                 gcagcaAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGca                                                  <  2:1291784/98‑1 (MQ=255)
                                                  cagcaAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGcaa                                                 >  2:3331503/1‑98 (MQ=255)
                                                  cagcaAACGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGcaa                                                 >  1:2555387/1‑98 (MQ=255)
                                                      aaaCGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAag                                             <  1:381029/98‑1 (MQ=255)
                                                        aCGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAagca                                           <  1:2267060/98‑1 (MQ=255)
                                                        aCGCGACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAagca                                           <  1:3155634/98‑1 (MQ=255)
                                                            gACAACACCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCaa                                       <  1:818163/98‑1 (MQ=255)
                                                             acaacaCCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCaaa                                      >  1:1812809/1‑98 (MQ=255)
                                                                acacCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAAcgc                                   >  1:2368402/1‑98 (MQ=255)
                                                                acacCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAAcgc                                   >  1:2794700/1‑98 (MQ=255)
                                                                 cacCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAAcgcg                                  <  1:3098587/98‑1 (MQ=255)
                                                                 cacCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAAcgcg                                  <  1:3041007/98‑1 (MQ=255)
                                                                 cacCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAAcgcg                                  >  1:270346/1‑98 (MQ=255)
                                                                  acCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGa                                 >  2:504257/1‑98 (MQ=255)
                                                                  acCATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGa                                 >  1:2674810/1‑98 (MQ=255)
                                                                    cATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGaca                               >  2:358315/1‑98 (MQ=255)
                                                                    cATCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGaca                               >  2:513941/1‑98 (MQ=255)
                                                                     aTCTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGacaa                              >  1:3547289/1‑98 (MQ=255)
                                                                       cTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGacaaca                            <  1:285058/98‑1 (MQ=255)
                                                                       cTTCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGacaaca                            >  1:1537970/1‑98 (MQ=255)
                                                                        ttCAACTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAacac                           >  1:1199361/1‑98 (MQ=255)
                                                                           aaCTAACGTAGAAGCACCGCAAACAAAACAAAACGCGCCCACACCACCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAAcaccac                        <  2:3384960/98‑1 (MQ=18)
                                                                           aaCTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAAcaccac                        >  1:1655023/1‑98 (MQ=255)
                                                                            aCTAAAGTGGAAACACCGCAATCAACAGCAAACGCAACAACCCCATCTTCACCTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACaccacc                       <  2:2853140/98‑1 (MQ=11)
                                                                            aCTAAAGTAGAAACACCGCAACCAGAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACaccacc                       <  2:1061062/98‑1 (MQ=255)
                                                                             cTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCt                      <  1:3241134/98‑1 (MQ=255)
                                                                             cTAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCt                      <  1:3370653/98‑1 (MQ=255)
                                                                             cTAAAGTAGAAGCACCGCAACCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCt                      <  2:2479445/98‑1 (MQ=255)
                                                                              tAAAGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCtt                     <  1:52675/98‑1 (MQ=255)
                                                                               aaaGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTc                    <  1:1217165/98‑1 (MQ=255)
                                                                               aaaGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTAAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCtac                    >  2:2801151/1‑96 (MQ=255)
                                                                                aaGTAGAAGCACCGCAACCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCg                   <  2:888042/98‑1 (MQ=255)
                                                                                 aGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGt                  <  1:2577098/98‑1 (MQ=255)
                                                                                 aGTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCGTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGt                  <  1:1268329/98‑1 (MQ=255)
                                                                                  gTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTc                 >  1:2804069/1‑98 (MQ=255)
                                                                                  gTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTc                 >  2:2716958/1‑98 (MQ=255)
                                                                                  gTAGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTc                 >  1:2528187/1‑98 (MQ=255)
                                                                                    aGAAGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTa               >  1:1505903/1‑98 (MQ=255)
                                                                                    aGAAGCACCGCAACCACAACCAAACGCAACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTa               <  2:2393160/98‑1 (MQ=255)
                                                                                      aaGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAAt             <  1:1046300/98‑1 (MQ=255)
                                                                                       aGCACCGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATg            >  1:682828/1‑98 (MQ=255)
                                                                                            cGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGac       >  2:3168595/1‑98 (MQ=255)
                                                                                            cGCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGac       <  1:754016/98‑1 (MQ=255)
                                                                                             gCAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGaca      >  2:3113929/1‑98 (MQ=255)
                                                                                              cAATCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGacac     <  1:2070623/98‑1 (MQ=255)
                                                                                               aaTCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGAcaca    <  2:369951/98‑1 (MQ=255)
                                                                                               aaTCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGAcaca    >  1:1940936/1‑98 (MQ=255)
                                                                                                aTCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACCAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGACACAt   >  1:1351590/1‑98 (MQ=255)
                                                                                                aTCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGACACAt   >  2:2101482/1‑98 (MQ=255)
                                                                                                aTCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGACACAt   >  1:3188849/1‑98 (MQ=255)
                                                                                                aTCAAAACCAAACGCGACAACACCATCTTCAACTAAAGTGGAAGCACCGCAACAAGCAGCAAACGCGACAACACCACCTTCGTCTAATGTAGACACAt   <  1:2975489/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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: