Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 697,673 G→T 92.5% A121S (GCT→TCT)  SAS0619 → putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX571857697,6730GT92.5% 583.9 / 37.4 158A121S (GCT→TCT) SAS0619putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (8/4);  new base T (82/64);  total (90/68)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

Only 100 of 158 total aligned reads displayed.

TTTTAATCTTCTTCTTATTTGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTGCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATCGTTACACCATTTTTAATT  >  BX571857/697578‑697765
                                                                                               |                                                                                            
ttttAATCTTCTTCTTATTTGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTGc                                                                                             <  1:2068563/97‑1 (MQ=255)
  ttAATCTTCTTCTTATTTGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCtt                                                                                           <  2:198715/97‑1 (MQ=255)
  ttAATCTTCTTCTTATTTGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCtt                                                                                           >  2:1987081/1‑97 (MQ=255)
    aaTCTTCTTCTTATTTGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTa                                                                                         <  1:563569/97‑1 (MQ=255)
     aTCTTCTTCTTATTTGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTaa                                                                                        <  2:1105574/97‑1 (MQ=255)
       cttcttcttATTTGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATCACATACCTACATGGCTAGGCTCTTTAAtt                                                                                      >  2:409688/1‑97 (MQ=255)
       cttcttcttATTTGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAAtt                                                                                      >  1:1180572/1‑97 (MQ=255)
       cttcttcttATTTGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAAtt                                                                                      >  1:2087846/1‑97 (MQ=255)
       cttcttcttATTTGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAAtt                                                                                      <  2:261696/97‑1 (MQ=255)
        ttcttcttATTTGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTa                                                                                     >  1:1396880/1‑97 (MQ=255)
        ttcttcttATTTGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTa                                                                                     >  1:310770/1‑97 (MQ=255)
         tcttcttATTTGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTAt                                                                                    <  2:446880/97‑1 (MQ=255)
           ttcttATTTGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTGCTTTAATTATGa                                                                                  <  1:820644/97‑1 (MQ=255)
            tcttATTTGGTATAGCTGTAACCATGATAGCTGGTGCATGCTCAACATTTGATGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGAc                                                                                 <  1:1654737/97‑1 (MQ=255)
             cttATTTGGTATAGCTGTAACCATGATAGCTGGTGCAGGTGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACa                                                                                <  2:1605213/97‑1 (MQ=255)
              ttATTTGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACAt                                                                               >  1:2256274/1‑97 (MQ=255)
                aTTTGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTa                                                                             <  2:1588433/97‑1 (MQ=255)
                aTTTGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTGCTTTAATTATGACATTa                                                                             >  1:389275/1‑97 (MQ=255)
                  ttGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGc                                                                           >  1:1539809/1‑97 (MQ=255)
                   tGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCg                                                                          >  1:1066330/1‑97 (MQ=255)
                   tGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCg                                                                          >  2:157762/1‑97 (MQ=255)
                   tGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCg                                                                          >  2:1786831/1‑97 (MQ=255)
                   tGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTGCTTTAATTATGACATTAGCg                                                                          >  2:906769/1‑97 (MQ=255)
                      tataGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGAtt                                                                       <  2:1336148/97‑1 (MQ=255)
                      tataGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGAtt                                                                       <  2:1417477/97‑1 (MQ=255)
                      tataGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGAtt                                                                       <  1:2089697/97‑1 (MQ=255)
                      tataGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGAtt                                                                       <  2:149451/97‑1 (MQ=255)
                      tataGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGAtt                                                                       <  1:1878511/97‑1 (MQ=255)
                        taGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGCCATTAGCGATTta                                                                     >  2:1651308/1‑97 (MQ=255)
                         aGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTtat                                                                    >  2:218524/1‑97 (MQ=255)
                         aGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTtat                                                                    <  1:1416132/97‑1 (MQ=255)
                          gCTGTAACCATGATAGCTGGTGCAGGCGCGATGTTTGAAGAAAGTTATAACATACCTACATGGGTAGGGTCTTTATTTGTGTCGTTAGGGGTTtata                                                                   >  2:1224999/1‑97 (MQ=255)
                          gCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTtata                                                                   >  2:1196536/1‑97 (MQ=255)
                           cTGTATCCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTatat                                                                  <  1:1908521/97‑1 (MQ=255)
                            tGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATAtt                                                                 >  2:576210/1‑97 (MQ=255)
                             gTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTa                                                                <  1:1518085/97‑1 (MQ=255)
                              tAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCCAGGTTCTTTAATTATGACATTAGCGATTTATATTAc                                                               >  2:2095253/1‑97 (MQ=255)
                               aaCCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACg                                                              >  2:1525037/1‑97 (MQ=255)
                                  cATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCCAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTg                                                           >  2:1709308/1‑97 (MQ=255)
                                   aTGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGc                                                          <  1:3273/97‑1 (MQ=255)
                                    tGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCt                                                         <  2:1813679/97‑1 (MQ=255)
                                    tGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCt                                                         <  1:35645/97‑1 (MQ=255)
                                    tGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCt                                                         >  2:1645175/1‑97 (MQ=255)
                                       tAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTAtt                                                      >  2:2077149/1‑97 (MQ=255)
                                         gCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAg                                                    <  2:1124035/97‑1 (MQ=255)
                                               gCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTta                                              >  1:478910/1‑97 (MQ=255)
                                               gCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTta                                              >  1:256180/1‑97 (MQ=255)
                                                cAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTtaa                                             >  2:539348/1‑97 (MQ=255)
                                                   gcgcAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTtaataa                                          >  1:1067935/1‑97 (MQ=255)
                                                    cgcAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATaaa                                         >  2:285647/1‑97 (MQ=255)
                                                    cgcAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATaaa                                         <  1:1196203/97‑1 (MQ=255)
                                                    cgcAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATaaa                                         <  1:1838393/97‑1 (MQ=255)
                                                    cgcAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTGCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATaaa                                         >  1:864818/1‑97 (MQ=255)
                                                    cgcAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTGCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATaaa                                         >  2:1528515/1‑97 (MQ=255)
                                                     gcAACATTGGAAGAAAGTTATAACAGCCCTGCATGGCTACGTGGATTAAGTATGACATAAGCGACTTATGTTATGTTGATATTAGAGTTtaataaga                                        >  2:863304/1‑95 (MQ=255)
                                                        aCATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTGCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAg                                     >  2:2235535/1‑97 (MQ=255)
                                                        aCATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTGCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAg                                     >  2:2235890/1‑97 (MQ=255)
                                                         cATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGt                                    >  2:1343758/1‑97 (MQ=255)
                                                         cATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGt                                    >  2:1436184/1‑97 (MQ=255)
                                                          aTTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTa                                   <  1:1561954/97‑1 (MQ=255)
                                                          aTTTGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTa                                   >  2:375002/1‑97 (MQ=255)
                                                           tttGCAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATACTACGTCGCTATTAGACTTTAATAAAATAGTAc                                  >  2:807374/1‑97 (MQ=255)
                                                            ttGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACg                                 >  1:430519/1‑97 (MQ=255)
                                                            ttGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACg                                 >  2:282461/1‑97 (MQ=255)
                                                             tGAAGAAAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGt                                >  2:1245922/1‑97 (MQ=255)
                                                              gaagaaAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTg                               <  1:1917656/97‑1 (MQ=255)
                                                                agaaAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCa                             >  1:989097/1‑97 (MQ=255)
                                                                 taaaGTTATAACATACCTACATGGCTAGGTTCTGTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCAc                            <  1:1018512/96‑1 (MQ=255)
                                                                 gagaGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCAc                            <  1:1018329/94‑1 (MQ=255)
                                                                 gaaAGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCAc                            <  1:1980960/97‑1 (MQ=255)
                                                                  aaaGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACt                           >  2:2009037/1‑97 (MQ=255)
                                                                   aaGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTa                          >  2:335135/1‑97 (MQ=255)
                                                                   aaGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTa                          >  1:2035869/1‑97 (MQ=255)
                                                                    aGTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAg                         <  2:1741568/97‑1 (MQ=255)
                                                                     gTTATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAgg                        <  2:1600690/97‑1 (MQ=255)
                                                                      ttATAACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGt                       <  2:1723430/97‑1 (MQ=255)
                                                                        ataACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTAt                     >  2:1559759/1‑97 (MQ=255)
                                                                        ataACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTAt                     <  2:1848477/97‑1 (MQ=255)
                                                                        ataACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTAt                     >  2:791364/1‑97 (MQ=255)
                                                                        ataACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTAt                     <  1:2162194/97‑1 (MQ=255)
                                                                         taACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATc                    <  2:1039575/97‑1 (MQ=255)
                                                                         taACATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATc                    >  1:2167753/1‑97 (MQ=255)
                                                                         taACATACCTACATGGCTAGGTGCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATc                    <  1:855885/97‑1 (MQ=255)
                                                                          aaCATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATCg                   <  2:490220/97‑1 (MQ=255)
                                                                          aaCATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATCg                   >  1:1770158/1‑97 (MQ=255)
                                                                          aaCATACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATCg                   >  1:1770190/1‑97 (MQ=255)
                                                                              tACCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATCGTTac               >  2:912600/1‑97 (MQ=255)
                                                                               aCCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATCGTTaca              >  1:1553231/1‑97 (MQ=255)
                                                                               aCCTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATCGTTaca              >  1:1187390/1‑97 (MQ=255)
                                                                                ccTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATCGTTacac             >  2:2169988/1‑97 (MQ=255)
                                                                                 cTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATCGTTACAcc            <  2:437944/97‑1 (MQ=255)
                                                                                 cTACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATCGTTACAcc            <  2:644783/97‑1 (MQ=255)
                                                                                  tACATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATCGTTACACCa           <  2:289791/97‑1 (MQ=255)
                                                                                   aCATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATCGTTACACCAt          >  1:2005731/1‑97 (MQ=255)
                                                                                   aCATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATCGTTACACCAt          <  1:713048/97‑1 (MQ=255)
                                                                                    cATGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATCGTTACACCAtt         >  1:488606/1‑97 (MQ=255)
                                                                                     aTGGCTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATCGTTACACCAttt        >  1:2205777/1‑97 (MQ=255)
                                                                                         cTAGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATCGTTACACCATTTTTaa    >  1:1691844/1‑97 (MQ=255)
                                                                                           aGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATCGTTACACCATTTTTAAtt  >  1:1599589/1‑97 (MQ=255)
                                                                                           aGGTTCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATCGTTACACCATTTTTAAtt  <  2:186158/97‑1 (MQ=255)
                                                                                               |                                                                                            
TTTTAATCTTCTTCTTATTTGGTATAGCTGTAACCATGATAGCTGGTGCAGGCGCAACATTTGAAGAAAGTTATAACATACCTACATGGCTAGGTGCTTTAATTATGACATTAGCGATTTATATTACGTTGCTATTAGACTTTAATAAAATAGTACGTGCACTAGGTATCGTTACACCATTTTTAATT  >  BX571857/697578‑697765

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: