Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 2,727,708 G→A 5.8% T430I (ACC→ATC)  SAS2534 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX5718572,727,7080GA 5.8% 893.2 / 35.1 225T430I (ACC→ATC) SAS2534hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (126/86);  new base A (8/5);  total (134/91)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 238 total aligned reads displayed.

ACAATGGCATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTGATGCAG  >  BX571857/2727612‑2727804
                                                                                                |                                                                                                
aCAATGGCATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAgg                                                                                                  >  2:829853/1‑97 (MQ=255)
aCAATGGCATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAgg                                                                                                  <  2:455644/97‑1 (MQ=255)
   aTGGCATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTg                                                                                               >  1:1728627/1‑97 (MQ=255)
     ggCATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTggg                                                                                             >  1:1041112/1‑97 (MQ=255)
      gCATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGt                                                                                            >  1:499915/1‑97 (MQ=255)
      gCATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGt                                                                                            >  2:1788172/1‑97 (MQ=255)
      gCATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGt                                                                                            >  1:1364983/1‑97 (MQ=255)
       cATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGtt                                                                                           >  1:382315/1‑97 (MQ=255)
       cATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGtt                                                                                           <  1:1776114/97‑1 (MQ=255)
         tGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTAc                                                                                         <  1:1565379/97‑1 (MQ=255)
          gAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTAcc                                                                                        >  2:1831060/1‑97 (MQ=255)
            gCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCAt                                                                                      <  2:2220439/97‑1 (MQ=255)
              tAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATaa                                                                                    <  2:1013852/97‑1 (MQ=255)
              tAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATaa                                                                                    <  2:581523/97‑1 (MQ=255)
              tAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGATTGGGTTACCATaa                                                                                    <  1:211142/97‑1 (MQ=255)
                  tttGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTca                                                                                <  2:154133/97‑1 (MQ=255)
                   ttGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTcac                                                                               >  2:2133039/1‑97 (MQ=255)
                       gATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACAttt                                                                           >  1:1825894/1‑97 (MQ=255)
                         tAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTaa                                                                         <  1:1395112/97‑1 (MQ=255)
                          aaCATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAAt                                                                        >  1:1337264/1‑97 (MQ=255)
                          aaCATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAAt                                                                        >  1:938955/1‑97 (MQ=255)
                          aaCATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAAt                                                                        >  1:761503/1‑97 (MQ=255)
                            catcatCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAACGGTTGGGTTACCATAAATCACATTTAATcc                                                                      >  2:1498319/1‑97 (MQ=255)
                               catcCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAAt                                                                   <  2:2076341/97‑1 (MQ=255)
                               catcCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAAt                                                                   >  1:2002194/1‑97 (MQ=255)
                                atcCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATc                                                                  >  1:561283/1‑97 (MQ=255)
                                 tcCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGATTGGGTTACCATAATTCACATTTAATCCAATCa                                                                 <  1:1526066/97‑1 (MQ=255)
                                  ccACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCAt                                                                <  1:425203/97‑1 (MQ=255)
                                  ccACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTCACCATAATTCACATTTAATCCAATCAt                                                                >  2:38828/1‑97 (MQ=255)
                                  ccACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATCTTCCCCATCTCCAATGAAGGGTGGGTTACCATAAATCACATTTAATCCAATCAt                                                                >  2:91023/1‑97 (MQ=255)
                                   cacaCTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATa                                                               >  2:617114/1‑97 (MQ=255)
                                   cacaCTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATa                                                               <  1:878834/97‑1 (MQ=255)
                                   cacaCTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATa                                                               >  1:2244473/1‑97 (MQ=255)
                                     cacTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATAcc                                                             >  1:464509/1‑97 (MQ=255)
                                       cTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTa                                                           >  2:1426049/1‑97 (MQ=255)
                                        tATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAg                                                          >  2:1579833/1‑97 (MQ=255)
                                        tATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAg                                                          <  1:372427/97‑1 (MQ=255)
                                        tATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAg                                                          <  1:160874/97‑1 (MQ=255)
                                        tATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAg                                                          <  2:2166301/97‑1 (MQ=255)
                                         aTCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGt                                                         <  1:1428383/97‑1 (MQ=255)
                                         aTCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGGTACCATAATTCACATTTAATCCAATCACACCTAGt                                                         >  1:22839/1‑97 (MQ=255)
                                            ttCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCt                                                      >  1:1683100/1‑97 (MQ=255)
                                            ttCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCt                                                      >  1:1811903/1‑97 (MQ=255)
                                              cATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGa                                                    <  2:925540/97‑1 (MQ=255)
                                              cATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGa                                                    >  1:825410/1‑97 (MQ=255)
                                              cATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGa                                                    >  1:403127/1‑97 (MQ=255)
                                               aTCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAg                                                   <  1:1306230/97‑1 (MQ=255)
                                                  ggTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCt                                                >  1:1701920/1‑97 (MQ=255)
                                                  ggTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCt                                                >  1:620197/1‑97 (MQ=255)
                                                  ggTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCt                                                <  2:33733/97‑1 (MQ=255)
                                                    tATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTAc                                              >  2:1106419/1‑97 (MQ=255)
                                                     aTCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTAcc                                             >  1:964034/1‑97 (MQ=255)
                                                     aTCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTAcc                                             >  2:955919/1‑97 (MQ=255)
                                                     aTCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTAcc                                             <  2:680028/97‑1 (MQ=255)
                                                       cAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGc                                           >  2:1287003/1‑97 (MQ=255)
                                                       cAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGc                                           <  2:657143/97‑1 (MQ=255)
                                                        aaaaGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCt                                          <  1:2061178/97‑1 (MQ=255)
                                                         aaaGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGATTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCtt                                         >  1:1254088/1‑97 (MQ=255)
                                                          aaGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTc                                        <  1:2190933/97‑1 (MQ=255)
                                                           aGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTcc                                       <  1:564138/97‑1 (MQ=255)
                                                           aGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGATTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTcc                                       >  1:1840923/1‑97 (MQ=255)
                                                            ggTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCa                                      <  1:1516881/97‑1 (MQ=255)
                                                             gTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCAt                                     >  1:1058167/1‑97 (MQ=255)
                                                             gTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCAt                                     >  2:1700212/1‑97 (MQ=255)
                                                                      aGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTaa                            >  1:984978/1‑97 (MQ=255)
                                                                      aGCCATTTTCCCCATCCCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCAAACATAGACCTGAGCCTACCGCTTCCATTAGTGTTaa                            >  2:1010199/1‑97 (MQ=255)
                                                                       gCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAAc                           >  2:488414/1‑97 (MQ=255)
                                                                         cctcttCCCCATCTCGCATAAAGGTTCGGCTACCATACTTCCCATTTAACCCCCTCATACCTAGTCCTCAGCCTACCGCTTCCATCAGTGTTAAccc                         >  2:354204/5‑97 (MQ=255)
                                                                         cATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAAccc                         >  2:2103003/1‑97 (MQ=255)
                                                                         cATCTTCCCCATCTCGAACAAAGGCGCGGCTACCATAATTCACATTTCCTCCCATCTTACCTAGTCCCGAGCCGCCCGCATCCACAAGACTTAAccc                         >  2:206011/1‑97 (MQ=255)
                                                                           ttttCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCaa                       >  2:896419/1‑97 (MQ=255)
                                                                           ttttCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCaa                       >  2:2155811/1‑97 (MQ=255)
                                                                            tttCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCaaa                      <  1:1177471/97‑1 (MQ=255)
                                                                             ttCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCaaaa                     >  2:970934/1‑97 (MQ=255)
                                                                             ttCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCaaaa                     <  2:1203133/97‑1 (MQ=255)
                                                                             ttCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCaaaa                     <  2:1044354/97‑1 (MQ=255)
                                                                                cccATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCt                  >  1:67801/1‑97 (MQ=255)
                                                                                cccATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCt                  <  1:2012081/97‑1 (MQ=255)
                                                                                 ccATCTCGAATAAAGGTTGGGTTACCATAATTCACCTTTAATCCAATCATACCTAGTCCTGACCCTACCGCTTCCATTAGTGTTAACCCAAACCCtt                 >  1:1698011/1‑97 (MQ=255)
                                                                                 ccATCTCGAATAAAGATTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCtt                 >  1:595809/1‑97 (MQ=255)
                                                                                  cATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTc                <  1:2147283/97‑1 (MQ=255)
                                                                                   aTCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCg               <  1:1662154/97‑1 (MQ=255)
                                                                                   aTCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCg               >  2:1709417/1‑97 (MQ=255)
                                                                                    tctcGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGc              >  2:1827894/1‑97 (MQ=255)
                                                                                      tcGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTa            <  2:1869164/97‑1 (MQ=255)
                                                                                       cGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAg           <  2:2247752/97‑1 (MQ=255)
                                                                                        gAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGt          >  1:510928/1‑97 (MQ=255)
                                                                                         aaTAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGtt         <  2:1823189/97‑1 (MQ=255)
                                                                                           tAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTGa       >  2:676372/1‑97 (MQ=255)
                                                                                           tAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTGa       >  1:1022998/1‑97 (MQ=255)
                                                                                           tAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTGa       >  2:2006718/1‑97 (MQ=255)
                                                                                            aaaGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTGAt      >  2:810754/1‑97 (MQ=255)
                                                                                            aaaGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTGAt      >  1:324054/1‑97 (MQ=255)
                                                                                             aaGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTGATg     >  1:1272966/1‑97 (MQ=255)
                                                                                             aaGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTGATg     >  2:370798/1‑97 (MQ=255)
                                                                                             aaGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTGATg     <  1:1527721/97‑1 (MQ=255)
                                                                                             aaGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTGATg     <  1:415910/97‑1 (MQ=255)
                                                                                             aaGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTGATg     <  1:2208743/97‑1 (MQ=255)
                                                                                               ggTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTGATGCa   >  2:1251440/1‑97 (MQ=255)
                                                                                                gTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTGATGCag  >  2:146437/1‑97 (MQ=255)
                                                                                                |                                                                                                
ACAATGGCATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTGATGCAG  >  BX571857/2727612‑2727804

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: ‑