Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 1,025,623 C→A 5.4% intergenic (+392/‑184) SAS0954 → / → SAS0955 amidase/putative protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX5718571,025,6230CA 5.4% 573.6 / 5.2 165intergenic (+392/‑184)SAS0954/SAS0955amidase/putative protease
Reads supporting (aligned to +/- strand):  ref base C (82/71);  new base A (2/8);  total (86/79)
Fisher's exact test for biased strand distribution p-value = 5.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 167 total aligned reads displayed.

TTGCTGGGTGGTTGTTTTTTATGTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGTTTTTAAACG  >  BX571857/1025529‑1025719
                                                                                              |                                                                                                
ttGCTGGGTGGTTGTTTTTTATGTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCaa                                                                                                <  1:1399162/97‑1 (MQ=255)
 tGCTGGGTGGTTGTTTTTTATGTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCaaa                                                                                               >  2:172823/1‑97 (MQ=255)
  gCTGGGTGGTTGTTTTTTATGTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTTAGGGCaaaa                                                                                              >  2:1710100/1‑97 (MQ=255)
  gCTGGGTGGTTGTTTTTTATGTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCaaaa                                                                                              >  2:1609511/1‑97 (MQ=255)
    tGGGTGGTTGGTTTTTATGTTATATTATAAATGAGCAAACCACACCACCTATTTATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCGGGGCaaaaaa                                                                                            <  2:954602/97‑1 (MQ=255)
      ggtggtTGTTTTTTATGTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAg                                                                                          >  1:872827/1‑97 (MQ=255)
        tggtTGTTTTTTATGTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTCTTCACGGCAACACAAggg                                                                                        >  2:400750/1‑97 (MQ=255)
        tggtTGTTTTTTATGTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGCGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAggg                                                                                        >  2:1514506/1‑97 (MQ=255)
          gttgttTTTTATGTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCa                                                                                      <  1:2221879/97‑1 (MQ=255)
          gttgttTTTTATGTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCa                                                                                      >  1:2126799/1‑97 (MQ=255)
           ttgttTTTTATGTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAg                                                                                     <  1:1241403/97‑1 (MQ=255)
            tgttTTTTATGTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGa                                                                                    >  1:1411730/1‑97 (MQ=255)
            tgttTTTTATGTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGa                                                                                    >  2:727669/1‑97 (MQ=255)
            tgttTTTTATGTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGa                                                                                    >  1:2206639/1‑97 (MQ=255)
            tgttTTTTATGTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGa                                                                                    >  1:699379/1‑97 (MQ=255)
                ttttATGTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGattat                                                                                <  2:648973/97‑1 (MQ=255)
                  ttATGTTATATTATAAATGCTCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTAttt                                                                              <  1:1363480/97‑1 (MQ=255)
                    aTGTTATATTATAAATGATCAAACCACGCCTCCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTaa                                                                            <  1:3129/97‑1 (MQ=255)
                      gTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAAt                                                                          >  1:1704434/1‑97 (MQ=255)
                      gTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAAt                                                                          >  2:787720/1‑97 (MQ=255)
                      gTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAAt                                                                          >  1:1211407/1‑97 (MQ=255)
                         atattataAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTTAGGGCAAAAAAAGGGCAGATTATTTAAATAag                                                                       >  2:473599/1‑97 (MQ=255)
                         atattataAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGTAAAAAAAAGGTCAGATTATTTAAATAag                                                                       >  2:1016353/1‑97 (MQ=255)
                         atattataAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAag                                                                       <  1:845169/97‑1 (MQ=255)
                         atattataAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAag                                                                       >  2:64984/1‑97 (MQ=255)
                          tattataAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCATGTTAAAATATTTTTTGATTATTTAAATAaga                                                                      >  2:1453479/1‑97 (MQ=255)
                          tattataAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAaga                                                                      <  2:1850574/97‑1 (MQ=255)
                           attataAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGTCAAAAAAAGGGCAGATTATTTAAATAagag                                                                     >  2:2038841/1‑97 (MQ=255)
                           attataAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAagag                                                                     <  1:782996/97‑1 (MQ=255)
                              ataAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTAAAAAAAAAAAAAGGGCAGATTATTTAAATAAGAGCaa                                                                  <  2:1053117/97‑1 (MQ=255)
                               taAATGATCAAACCGCACCACCTAGGAATTTGGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCaaa                                                                 <  2:969066/97‑1 (MQ=255)
                               taAATGATCAAACCCCACCACCTATTAATTTAGGAGTGTGTGTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCaaa                                                                 <  2:2182658/97‑1 (MQ=255)
                                aaaTGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAac                                                                >  1:2263476/1‑97 (MQ=255)
                                aaaTGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAac                                                                <  1:1193925/97‑1 (MQ=255)
                                 aaTGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGTGCAAaca                                                               >  2:52702/1‑97 (MQ=255)
                                 aaTGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAaca                                                               >  1:2330707/1‑97 (MQ=255)
                                 aaTGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGCCCAAAAAAGGGCAGATTATTTAAATAAGAGCAAaca                                                               >  2:58419/1‑97 (MQ=255)
                                   tGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACg                                                             <  1:1114956/97‑1 (MQ=255)
                                   tGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACg                                                             <  1:132477/97‑1 (MQ=255)
                                   tGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACg                                                             >  1:1479882/1‑97 (MQ=255)
                                    gATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAATGGGCAGATTATTTAAATAAGAGCAAACACGa                                                            >  2:1810081/1‑97 (MQ=255)
                                     aTCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAg                                                           >  2:538558/1‑97 (MQ=255)
                                     aTCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAg                                                           >  2:1956762/1‑97 (MQ=255)
                                       cAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTg                                                         >  2:2153937/1‑97 (MQ=255)
                                        aaaCCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTTATGGGCAAAAAAGGGCAGGTTATTTAAAGAAGAGGAAATAGGTGTgg                                                        >  2:7749/1‑97 (MQ=25)
                                         aaCCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGa                                                       >  2:1445190/1‑97 (MQ=255)
                                          accacaccacCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGTCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACAAGAGTGGaa                                                      >  2:136542/1‑97 (MQ=255)
                                          accacaccacCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGaa                                                      <  2:292848/97‑1 (MQ=255)
                                          accacaccacCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTACAAAGAAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGaa                                                      <  2:178891/97‑1 (MQ=255)
                                           ccacgccacCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGaaa                                                     <  2:842208/97‑1 (MQ=255)
                                           ccacaccacCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGaaa                                                     >  1:1131482/1‑97 (MQ=255)
                                           ccacaccacCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGaaa                                                     <  2:1718853/97‑1 (MQ=255)
                                            cacaccacCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGaaaa                                                    <  1:930332/97‑1 (MQ=255)
                                             acaccacCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAat                                                   <  1:2500463/97‑1 (MQ=255)
                                             acaccacCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAat                                                   >  2:1215678/1‑97 (MQ=255)
                                               accaccTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAatat                                                 >  1:1748161/1‑97 (MQ=255)
                                                ccaccTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATAtt                                                <  1:2222138/97‑1 (MQ=255)
                                                  accTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATtta                                              >  2:1715106/1‑97 (MQ=255)
                                                   ccTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATttaa                                             <  2:2136337/97‑1 (MQ=255)
                                                   ccTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATttaa                                             >  2:2355548/1‑97 (MQ=255)
                                                    cTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTTAGGGCAAATAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATttaaa                                            >  2:1146810/1‑97 (MQ=255)
                                                    cTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATttaaa                                            >  1:1667782/1‑97 (MQ=255)
                                                      aTTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATttaaagt                                          <  1:1101781/97‑1 (MQ=255)
                                                          aTTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATttaaagttaaa                                      <  2:1109663/97‑1 (MQ=255)
                                                            ttAGGAGTGTGGTTATTTTTGTTTTTTTTTTTGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAATCACGAGTGGAAAATATTtaaagttaaagt                                    >  2:797100/1‑97 (MQ=255)
                                                             tAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTg                                   >  1:130872/1‑97 (MQ=255)
                                                                gAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATa                                >  1:932065/1‑97 (MQ=255)
                                                                gAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATa                                <  1:973772/97‑1 (MQ=255)
                                                                 aGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATaa                               <  2:1902709/97‑1 (MQ=255)
                                                                   tgtgGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGTCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATaaaa                             <  1:9241/97‑1 (MQ=255)
                                                                    gtgGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATaaaaa                            <  1:1885919/97‑1 (MQ=255)
                                                                    gtgGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATaaaaa                            >  1:2225560/1‑97 (MQ=255)
                                                                    gtgGTTATTTTTGTTTTATTTCAGGGAAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATaaaaa                            <  1:2062488/97‑1 (MQ=255)
                                                                     tgGTTATTTTTGTTTTTTTTCTGGGCAAATAAAGGGCTGATTATTTAAATAAGAGTAAACACGAGTGGAAAATATTTACAGTTAAAGTGATcaaaat                           >  2:1323352/1‑96 (MQ=255)
                                                                     tgGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAAc                           >  1:213833/1‑97 (MQ=255)
                                                                     tgGTTATTTTTGTTTTTTTTCAGGACAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAAc                           <  2:1111184/97‑1 (MQ=255)
                                                                      ggTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACt                          <  1:348261/97‑1 (MQ=255)
                                                                      gattATTTTTGTTTTTAATAAGAAAAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACt                          <  2:2061972/95‑1 (MQ=255)
                                                                       gTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACtt                         >  1:638815/1‑97 (MQ=255)
                                                                       gTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGACCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAAAtt                         >  1:703413/1‑97 (MQ=255)
                                                                         tATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTg                       >  1:474380/1‑97 (MQ=255)
                                                                           tttttGTTTTTTTTTAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGat                     >  2:351517/1‑97 (MQ=255)
                                                                           tttttGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGat                     <  2:1944135/97‑1 (MQ=255)
                                                                           tttttGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGat                     >  2:307755/1‑97 (MQ=255)
                                                                            ttttGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGata                    <  1:2131916/97‑1 (MQ=255)
                                                                              ttGTTTTTTTTCAGGGCAAAAAAAGGGAAGATTATATAAATAAGAGCAAAAACGAGTGGAAAATATTTAAAGTTAAAGTGATACAAACTTTGATAtt                  >  2:2541276/1‑97 (MQ=255)
                                                                                 ttttttttCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTAcc               >  1:142142/1‑96 (MQ=255)
                                                                                  tttttttCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAg              <  1:1300377/97‑1 (MQ=255)
                                                                                   ttttttCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGc             <  2:1615330/97‑1 (MQ=255)
                                                                                     ttttCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGt           <  1:2245842/97‑1 (MQ=255)
                                                                                      tttCATGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATGTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGtt          <  1:2048513/97‑1 (MQ=255)
                                                                                       ttCAGGGAAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGttt         <  2:842743/97‑1 (MQ=255)
                                                                                        tCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGtttt        >  1:2392275/1‑97 (MQ=255)
                                                                                            ggCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGTTTTTaaa    <  2:1303846/97‑1 (MQ=255)
                                                                                             gCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGTTTTTAAAc   <  1:2250435/97‑1 (MQ=255)
                                                                                             gCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGTTTTTAAAc   <  1:1645764/97‑1 (MQ=255)
                                                                                             gCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGTTTTTAAAc   >  1:1802731/1‑97 (MQ=255)
                                                                                             gCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGTTTTTAAAc   >  1:1831755/1‑97 (MQ=255)
                                                                                             gCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGTTTTTAAAc   <  1:543213/97‑1 (MQ=255)
                                                                                             gCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGTTTTTAAAc   <  1:543191/97‑1 (MQ=255)
                                                                                              |                                                                                                
TTGCTGGGTGGTTGTTTTTTATGTTATATTATAAATGATCAAACCACACCACCTATTAATTTAGGAGTGTGGTTATTTTTGTTTTTTTTCAGGGCAAAAAAAGGGCAGATTATTTAAATAAGAGCAAACACGAGTGGAAAATATTTAAAGTTAAAGTGATAAAAACTTTGATATTACAGCGTTTTTAAACG  >  BX571857/1025529‑1025719

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