Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 627,445 A→T 13.6% intergenic (‑133/‑35) SAS0546 ← / → SAS0547 conserved hypothetical protein/putative phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX571857627,4450AT13.6% 955.0 / 125.8 251intergenic (‑133/‑35)SAS0546/SAS0547conserved hypothetical protein/putative phosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (116/99);  new base T (21/13);  total (139/112)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 256 total aligned reads displayed.

TAAATTTTATCATAAAGCAAAGGGGTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTAAAAACGTTAAAATCGTATT  >  BX571857/627349‑627541
                                                                                                |                                                                                                
    ttttATCATAAAGCAAAGGGGTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGTATGa                                                                                              <  1:1050599/97‑1 (MQ=255)
    ttttATCATAAAGCAAAGGGGTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGa                                                                                              <  2:2219844/97‑1 (MQ=255)
    ttttATCATAAAGCAAAGGGGTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGa                                                                                              >  2:3175762/1‑97 (MQ=255)
    ttttATCATAAAGCAAAGGGGTAACGATAAAATTATGATTACAAATTGGTGACGCGGCATCATGCAATAAAATGGCGTCTAATTATACCGTGAATgg                                                                                              >  2:1596233/1‑96 (MQ=255)
     tttATCATAAAGCAACGGGGTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGAt                                                                                             >  2:36346/1‑97 (MQ=255)
         tCATAAAGCAAAGGGGTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATtaat                                                                                         <  2:311609/97‑1 (MQ=255)
           aTAAAGCAAAGGGGTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATtaataa                                                                                       <  1:1256874/97‑1 (MQ=255)
           aTAAAGCAAAGGGGTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATtaataa                                                                                       >  1:1116844/1‑97 (MQ=255)
           aTAAAGCAAAGGGGTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATtaataa                                                                                       <  1:208644/97‑1 (MQ=255)
           aTAAAGCAAAGGGGTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATtaataa                                                                                       >  2:1224089/1‑97 (MQ=255)
              aaGCAAAGGGGTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATATCGTGAATGATTAATAAGag                                                                                    >  2:693190/1‑96 (MQ=255)
              aaGCAAAGGGGTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGAt                                                                                    >  1:2991676/1‑97 (MQ=255)
               aGCAAAGGGGTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGAtt                                                                                   <  2:2030267/97‑1 (MQ=255)
                gCAAAGGGGTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGTATGATTAATAAGAttt                                                                                  >  1:1522641/1‑97 (MQ=255)
                gCAAAGGGGTAACGATAAAATTATGATTACAAATTGGTAACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGAttt                                                                                  >  1:2285532/1‑97 (MQ=255)
                 cAAAGGGGTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTta                                                                                 <  1:3291882/97‑1 (MQ=255)
                   aaGGGGTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGTATGATTAATAAGATTtata                                                                               >  2:1721696/1‑97 (MQ=255)
                   aaGGGGTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTtata                                                                               <  2:1628871/97‑1 (MQ=255)
                   aaGGGGTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTtata                                                                               <  1:2896248/97‑1 (MQ=255)
                     ggggTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATAtt                                                                             >  2:144281/1‑97 (MQ=255)
                     ggggTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATAtt                                                                             >  2:1865841/1‑97 (MQ=255)
                     ggggTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATAtt                                                                             >  1:262467/1‑97 (MQ=255)
                       ggTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGTATGATTAATAAGATTTATATTAc                                                                           <  1:2681519/97‑1 (MQ=255)
                       ggTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTAc                                                                           >  1:1949828/1‑97 (MQ=255)
                          aaCGATAAAATTCAGATTACATAATGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACagg                                                                        >  2:1737105/1‑97 (MQ=255)
                             gATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGTATGATTAATAAGATTTATATTACaggagg                                                                     >  2:869046/1‑97 (MQ=255)
                             gATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACaggagg                                                                     <  1:758230/97‑1 (MQ=255)
                              aTAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAggagga                                                                    >  1:3236001/1‑97 (MQ=255)
                               tAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGAc                                                                   >  2:924964/1‑97 (MQ=255)
                               tAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGAc                                                                   <  1:1713055/97‑1 (MQ=255)
                                aaaaTTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACa                                                                  >  2:2304344/1‑97 (MQ=255)
                                aaaaTTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACa                                                                  >  1:1830392/1‑97 (MQ=255)
                                aaaaTTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACa                                                                  <  2:618016/97‑1 (MQ=255)
                                aaaaTTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACa                                                                  >  2:2878380/1‑97 (MQ=255)
                                aaaaTTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACa                                                                  >  2:2965677/1‑97 (MQ=255)
                                 aaaTTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACAt                                                                 <  1:8057/97‑1 (MQ=255)
                                     tATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGTATGATTAATAAGATTTATATTACAGGAGGACATTATg                                                             >  1:1013080/1‑97 (MQ=255)
                                     tATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATg                                                             <  1:1982894/97‑1 (MQ=255)
                                        gATTACAAAATGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCt                                                          <  2:1677590/97‑1 (MQ=255)
                                         aTTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTg                                                         >  1:1754967/1‑97 (MQ=255)
                                         aTTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTg                                                         >  2:2832865/1‑97 (MQ=255)
                                         aTTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTg                                                         <  1:361201/97‑1 (MQ=255)
                                          ttACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGa                                                        <  1:858803/97‑1 (MQ=255)
                                           tACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGAt                                                       >  2:72343/1‑97 (MQ=255)
                                           tACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGAt                                                       >  2:3336969/1‑97 (MQ=255)
                                                aTTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATttatt                                                  >  2:574372/1‑97 (MQ=255)
                                                 ttGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACATGAGGACATTATGGCTGATttatta                                                 <  2:807836/97‑1 (MQ=255)
                                                 ttGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATttatta                                                 >  2:2496933/1‑97 (MQ=255)
                                                 ttGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATttatta                                                 <  2:1329220/97‑1 (MQ=255)
                                                  tGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTaa                                                >  1:1554062/1‑97 (MQ=255)
                                                  tGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTAAAGGAGGACATTATGGCTGATTTATTaa                                                >  2:2391203/1‑97 (MQ=255)
                                                     tGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATg                                             <  1:2164367/97‑1 (MQ=255)
                                                      gACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGTATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGt                                            <  1:1566749/97‑1 (MQ=255)
                                                        cGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTAt                                          >  2:3089507/1‑97 (MQ=255)
                                                         gTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTAtt                                         >  1:1409599/1‑97 (MQ=255)
                                                          tGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTa                                        <  1:2109762/97‑1 (MQ=255)
                                                          tGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTa                                        <  1:237437/97‑1 (MQ=255)
                                                           ggCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTaa                                       <  1:1132518/97‑1 (MQ=255)
                                                            gCATTATGAAATAAAATGGCGTATAATTATACCGTGTATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTaaa                                      >  1:2822459/1‑97 (MQ=255)
                                                             cATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTaaaa                                     <  1:694960/97‑1 (MQ=255)
                                                             cATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTaaaa                                     <  2:2248588/97‑1 (MQ=255)
                                                               ttATGAAATAAAATGGCGTATAATTATACCGTGAAAGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTAATAAAAGa                                   >  1:132121/1‑97 (MQ=255)
                                                                 aTGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACa                                 >  2:2373579/1‑97 (MQ=255)
                                                                  tGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACaa                                >  2:2451116/1‑97 (MQ=255)
                                                                  tGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACaa                                <  2:1988370/97‑1 (MQ=255)
                                                                   gAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACaaa                               >  2:1433097/1‑97 (MQ=255)
                                                                   gAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACaaa                               <  2:3051958/97‑1 (MQ=255)
                                                                    aaaTAAAATGGCGTATAATTATACCGTGTATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAAc                              <  2:1878830/97‑1 (MQ=255)
                                                                      aTAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACtt                            >  1:1497159/1‑97 (MQ=255)
                                                                      aTAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACtt                            <  1:47133/97‑1 (MQ=255)
                                                                      aTAAAATGGAGTATAATTATACCGTGAATGATTAAAAAGATTTATATTACAGGAGGACAATATGGATGATTTATTAAATGTATTAAAAGAAAAACtt                            >  2:380505/1‑97 (MQ=255)
                                                                        aaaaTGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTc                          >  1:359607/1‑97 (MQ=255)
                                                                         aaaTGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCt                         >  1:3131526/1‑97 (MQ=255)
                                                                         aaaTGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCt                         <  1:2463673/97‑1 (MQ=255)
                                                                            tGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGt                      >  2:2772969/1‑97 (MQ=255)
                                                                            tGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGt                      >  2:3338093/1‑97 (MQ=255)
                                                                             ggCGTATAATTATACCGTGTATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTa                     >  1:901741/1‑97 (MQ=255)
                                                                             ggCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTa                     <  2:1345124/97‑1 (MQ=255)
                                                                             ggCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTa                     >  1:2203235/1‑97 (MQ=255)
                                                                             ggCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTa                     >  2:1163852/1‑97 (MQ=255)
                                                                             ggCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTa                     >  2:3239846/1‑97 (MQ=255)
                                                                              gCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTaa                    >  2:2680687/1‑97 (MQ=255)
                                                                               cGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTaaa                   <  1:2831684/97‑1 (MQ=255)
                                                                                gTATAATTATACCGTGTATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTaaaa                  >  1:550614/1‑97 (MQ=255)
                                                                                   taATTATACCGTGTATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTAAAAACg               >  2:1050777/1‑97 (MQ=255)
                                                                                   taATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTAAAAACg               >  2:2267284/1‑97 (MQ=255)
                                                                                     aTTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCCGATTTATTAAATGTATTAAAAGACAACCTTTCTGGTAAAAACGtt             >  2:3270044/1‑97 (MQ=255)
                                                                                       tataCCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTAAAAACGTTaa           <  1:1970299/97‑1 (MQ=255)
                                                                                       tataCCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTAAAAACGTTaa           >  1:1199523/1‑97 (MQ=255)
                                                                                       tataCCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTAAAAACGTTaa           <  1:2440878/97‑1 (MQ=255)
                                                                                        ataCCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTAAAAACGTTaaa          >  1:735077/1‑97 (MQ=255)
                                                                                        ataCCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTAAAAACGTTaaa          >  1:3114644/1‑97 (MQ=255)
                                                                                         taCCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTAAAAACGTTaaaa         <  1:2965677/97‑1 (MQ=255)
                                                                                         taCCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTAAAAACGTTaaaa         <  1:2791989/97‑1 (MQ=255)
                                                                                         taCCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTAAAAACGTTaaaa         >  2:2106492/1‑97 (MQ=255)
                                                                                          aCCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTAAAAACGTTAAAAt        >  1:135494/1‑97 (MQ=255)
                                                                                           ccGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTAAAAACGTTAAAATc       >  1:1856552/1‑97 (MQ=255)
                                                                                            cGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTAAAAACGTTAAAATCg      <  2:558619/97‑1 (MQ=255)
                                                                                            cGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTAAAAACGTTAAAATCg      >  1:964295/1‑97 (MQ=255)
                                                                                                aaTGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTAAAAACGTTAAAATCGTAtt  >  2:2484857/1‑97 (MQ=255)
                                                                                                |                                                                                                
TAAATTTTATCATAAAGCAAAGGGGTAACGATAAAATTATGATTACAAATTGGTGACGTGGCATTATGAAATAAAATGGCGTATAATTATACCGTGAATGATTAATAAGATTTATATTACAGGAGGACATTATGGCTGATTTATTAAATGTATTAAAAGACAAACTTTCTGGTAAAAACGTTAAAATCGTATT  >  BX571857/627349‑627541

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