Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 1,408,208 G→C 5.3% G240A (GGT→GCT)  SAS1312 → tryptophan synthase beta chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX5718571,408,2080GC 5.3% 1157.1 / 48.9 281G240A (GGT→GCT) SAS1312tryptophan synthase beta chain
Reads supporting (aligned to +/- strand):  ref base G (152/114);  new base C (7/8);  total (159/122)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-01

Only 100 of 303 total aligned reads displayed.

GTTAGAGATTTTCAGAGTGTGATTGGTAAAGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAGGTCAAGGCGATGATACTGATAA  >  BX571857/1408111‑1408304
                                                                                                 |                                                                                                
  tAGAGATTTTCAGAGTGTGATTGGTAAAGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGgtggtg                                                                                                <  1:865199/98‑1 (MQ=255)
  tAGAGATTTTCAGAGTGTGATTGGTAAAGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGgtggtg                                                                                                >  2:1553170/1‑98 (MQ=255)
  tAGAGATTTTCAGAGTGTGATTGGTAAAGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGgtggtg                                                                                                >  2:652662/1‑98 (MQ=255)
  tAGAGATTTTCAGAGTGTGATTGGTAAAGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGgtggtg                                                                                                <  1:2909636/98‑1 (MQ=255)
  tAGAGATTTTCAGAGTGTGATTGGTAAAGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGgtggtg                                                                                                >  2:2538033/1‑98 (MQ=255)
     agaTTTTCAGAGTGTGATTGGTAAAGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGcc                                                                                             >  2:1865654/1‑97 (MQ=255)
       aTTTTCAGAGTGTGATTGGTAAAGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCa                                                                                           <  1:970773/98‑1 (MQ=255)
            cAGAGTGGGATGGGTAAAGAAATAAGAGGGCGGTTGGTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGc                                                                                      <  2:2000016/98‑1 (MQ=255)
              gagTGTGATTGGTAAAGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCaa                                                                                    <  1:3254425/98‑1 (MQ=255)
               agTGTGATTGGTAAAGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAAt                                                                                   >  1:3212912/1‑98 (MQ=255)
                 tgtgATTGGTAAAGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCg                                                                                 >  2:2453420/1‑98 (MQ=255)
                    taTTGGTAAAGAAATGAAATTACAGATATTGAAGAAAGAAGGGAGAATACCGGATGCAATTGTAGAATGTATCGGTGGTGGCTCAAATGCATTCGGTa                                                                              <  2:376137/97‑1 (MQ=255)
                     aTTGGTAAAGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTAc                                                                             >  2:1240002/1‑98 (MQ=255)
                     aTTGGTAAAGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTAc                                                                             <  1:1262682/98‑1 (MQ=255)
                      ttGGTAAAGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACa                                                                            >  1:2560575/1‑98 (MQ=255)
                      ttGGTAAAGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACa                                                                            >  1:1630504/1‑98 (MQ=255)
                       tcgTAAAGATATAAAAGCACAGATAGTGAGTAAAGAAGGGAGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACAt                                                                           <  2:1485144/96‑1 (MQ=255)
                          tAAAGAAATAAAATCACGGATATTGGAGAGAGAAGGTCAAATACCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACAtttt                                                                        <  2:49953/98‑1 (MQ=255)
                          tAAAGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACAtttt                                                                        >  1:569079/1‑98 (MQ=255)
                            aaGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTAt                                                                      >  1:1202468/1‑98 (MQ=255)
                            aaGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGCTGGCTCAAATGCAATCGGTACATTTTAt                                                                      <  1:790366/98‑1 (MQ=255)
                             aGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATc                                                                     >  1:1788101/1‑98 (MQ=255)
                              gAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATcc                                                                    <  1:2359914/98‑1 (MQ=255)
                              gAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATcc                                                                    >  1:3002368/1‑98 (MQ=255)
                              gAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATcc                                                                    >  2:3158699/1‑98 (MQ=255)
                               aaaTAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCa                                                                   <  1:2135358/98‑1 (MQ=255)
                               aaaTAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCa                                                                   <  1:1977475/98‑1 (MQ=255)
                                aaTAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCAt                                                                  >  1:388503/1‑98 (MQ=255)
                                aaTAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCAt                                                                  >  1:2164132/1‑98 (MQ=255)
                                aaTAAAAGCACAGATATTGAAGAGAGAAGGTCGTCTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCAt                                                                  <  2:2746080/98‑1 (MQ=255)
                                 aTAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCAtt                                                                 >  1:1045714/1‑98 (MQ=255)
                                  tAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCAttt                                                                >  1:1339839/1‑98 (MQ=255)
                                  tAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGCTGGCTCAAATGCAATCGGTACATTTTATCCAttt                                                                <  2:1323601/98‑1 (MQ=255)
                                  tAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGCTGGCTCAAATGCAATCGGTACATTTTATCCAttt                                                                >  2:33607/1‑98 (MQ=255)
                                   aaaaTCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATtta                                                               >  1:1083390/1‑98 (MQ=255)
                                   aaaaTCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATtta                                                               >  2:2028665/1‑98 (MQ=255)
                                   aaaaTCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATtta                                                               <  1:2179340/98‑1 (MQ=255)
                                     aaTCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCGCAAATGCAATCGGTACATTTTATCCATttatt                                                             >  2:3134423/1‑98 (MQ=255)
                                       tCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTaa                                                           <  1:491707/98‑1 (MQ=255)
                                          caGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAga                                                        >  2:2331500/1‑98 (MQ=255)
                                          caGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAga                                                        >  2:3080935/1‑98 (MQ=255)
                                          caGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAga                                                        >  1:2810292/1‑98 (MQ=255)
                                           aGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAgat                                                       >  1:1724233/1‑98 (MQ=255)
                                             atatTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAACgatga                                                     >  2:616307/1‑98 (MQ=255)
                                              tatTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAgatgat                                                    >  1:1323713/1‑98 (MQ=255)
                                               atTGAAGAATGAAGGTCGATTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGatgatg                                                   <  2:2190470/98‑1 (MQ=255)
                                               atTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGatgatg                                                   >  2:215429/1‑98 (MQ=255)
                                               atTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGatgatg                                                   >  1:2731025/1‑98 (MQ=255)
                                                 tGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGtt                                                 >  2:1890265/1‑98 (MQ=255)
                                                  gaagaaAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTg                                                <  2:456660/98‑1 (MQ=255)
                                                  gaagaaAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTg                                                <  2:2080708/98‑1 (MQ=255)
                                                    agaaagaaGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCa                                              >  2:181332/1‑98 (MQ=255)
                                                    agaaagaaGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTAGCCATTTATTAAAGATGATGTTGCa                                              <  2:3024358/98‑1 (MQ=255)
                                                     gaaagaaGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCAt                                             >  2:2387170/1‑98 (MQ=255)
                                                     gaaagaaGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCAt                                             >  2:2385401/1‑98 (MQ=255)
                                                     gaaagaaGGTCGACTTCCGGATGCAATTGTAGCATGTACCGGAGGTGGCTCAAATGCAAACGGTACATCTTATCCATTTACTAAAGATGACGTTCCAt                                             >  2:922662/1‑98 (MQ=255)
                                                       aagaagGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATta                                           >  1:631141/1‑98 (MQ=255)
                                                       aagaagGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGCTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATta                                           <  2:1944798/98‑1 (MQ=255)
                                                        agaagGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATtat                                          <  1:2615598/98‑1 (MQ=255)
                                                        agaagGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATtat                                          >  2:716915/1‑98 (MQ=255)
                                                        agaagGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATtat                                          >  1:941015/1‑98 (MQ=255)
                                                         gaagGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATtata                                         <  1:1851398/98‑1 (MQ=255)
                                                         gaagGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATtata                                         >  2:1739751/1‑98 (MQ=255)
                                                            ggTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACgg                                      <  1:3068778/98‑1 (MQ=255)
                                                            ggTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGCTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACgg                                      >  1:2266208/1‑98 (MQ=255)
                                                             gTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGCTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGgt                                     >  1:341437/1‑98 (MQ=255)
                                                              tccaCTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGgtg                                    <  2:2222975/95‑1 (MQ=255)
                                                              tCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGgtg                                    <  1:1548982/98‑1 (MQ=255)
                                                              tCGACTTCCGGATGCAATTGTAGCATGTATCGGTGCTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGgtg                                    <  1:209308/98‑1 (MQ=255)
                                                                 aCTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTg                                 >  1:3171699/1‑98 (MQ=255)
                                                                 aCTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCACTTATTAAAGATGATGTTGCATTATACGGTGTTg                                 >  1:1484227/1‑98 (MQ=255)
                                                                  ttgccGGATGCGATTGTAGTGTGTATCGGTGGTGGCGCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGa                                <  2:874719/95‑1 (MQ=255)
                                                                  cTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGa                                >  1:660864/1‑98 (MQ=255)
                                                                   ttCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGaa                               <  1:2293176/98‑1 (MQ=255)
                                                                   ttCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGaa                               <  2:809597/98‑1 (MQ=255)
                                                                   ttCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGaa                               >  2:1037728/1‑98 (MQ=255)
                                                                   ttCCAGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAAAGAAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGaa                               <  2:1694246/98‑1 (MQ=255)
                                                                   tcccGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGaa                               <  2:85032/96‑1 (MQ=255)
                                                                    tCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTCTATCCATTTATTAAAGATGATGTTGCATTATACGGGGTTGAAg                              >  2:1185724/1‑98 (MQ=255)
                                                                      cGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGcc                            >  1:837230/1‑98 (MQ=255)
                                                                         aTGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCa                         >  1:577670/1‑98 (MQ=255)
                                                                          tGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAg                        <  2:2927273/98‑1 (MQ=255)
                                                                           gCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAgg                       <  1:294424/98‑1 (MQ=255)
                                                                           gCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAgg                       <  2:2241670/98‑1 (MQ=255)
                                                                            cAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAGGt                      <  2:279008/98‑1 (MQ=255)
                                                                             aaTTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAGGTc                     >  2:1398461/1‑98 (MQ=255)
                                                                              aTTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAGGTCa                    >  2:1387075/1‑98 (MQ=255)
                                                                               ttGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAGGTCaa                   >  1:444593/1‑98 (MQ=255)
                                                                               ttGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAGGTCaa                   <  1:2453383/98‑1 (MQ=255)
                                                                                    gCATGTATCGGTGCTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAGGTCAAGGCga              >  2:517463/1‑98 (MQ=255)
                                                                                      aTGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGATGCCGCAGGTCATGGTgatg            <  2:3326202/98‑1 (MQ=255)
                                                                                       tGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAGGTCAAGGCgatga           <  2:1649574/98‑1 (MQ=255)
                                                                                         tATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAGGTCAAGGCGATGATa         >  1:2295797/1‑98 (MQ=255)
                                                                                           tCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAGGTCAAGGCGATGATACt       >  2:1703549/1‑98 (MQ=255)
                                                                                            cGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAGGTCAAGGCGATGATACTg      >  2:2949205/1‑98 (MQ=255)
                                                                                            cGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAGGTCAAGGCGATGATACTg      >  1:2015218/1‑98 (MQ=255)
                                                                                               tggtggCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAGGTCAAGGCGATGATACTGATa   >  1:1532222/1‑98 (MQ=255)
                                                                                                ggtggCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAGGTCAAGGCGATGATACTGATaa  <  1:2391862/98‑1 (MQ=255)
                                                                                                ggtggCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAGGTCAAGGCGATGATACTGATaa  <  1:1924439/98‑1 (MQ=255)
                                                                                                ggtggCTCAAAGGCAATCGGTACGTTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAGGTCAAGGCGATGATACTGATaa  <  1:40466/98‑1 (MQ=255)
                                                                                                 |                                                                                                
GTTAGAGATTTTCAGAGTGTGATTGGTAAAGAAATAAAATCACAGATATTGAAGAAAGAAGGTCGACTTCCGGATGCAATTGTAGCATGTATCGGTGGTGGCTCAAATGCAATCGGTACATTTTATCCATTTATTAAAGATGATGTTGCATTATACGGTGTTGAAGCCGCAGGTCAAGGCGATGATACTGATAA  >  BX571857/1408111‑1408304

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