Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BX571857 321,709 T→A 5.2% D162E (GAT→GAA SAS0274 → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BX571857321,7090TA 5.2% 672.9 / 5.5 165D162E (GAT→GAASAS0274conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (78/77);  new base A (5/5);  total (83/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 165 total aligned reads displayed.

TTATTATATGTACAAAAAATTTGGGGTTATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGCAGTATGATACCGG  >  BX571857/321616‑321804
                                                                                             |                                                                                               
ttattaTATGTACAAAAAATTTGGGGTTATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGaa                                                                                              >  1:2000041/1‑97 (MQ=39)
ttattaTATGTACAAAAAATTTGGGGTTATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGaa                                                                                              <  2:1570031/97‑1 (MQ=39)
ttattaTATGTACAAAAAATTTGGGGTTATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGaa                                                                                              <  1:1857656/97‑1 (MQ=38)
  attatATGTACAAAAAATTTGGGGTTATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGc                                                                                            >  2:1062838/1‑97 (MQ=39)
    tataTGTACAAAAAATTTGGGGTTATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTg                                                                                          >  2:1990551/1‑97 (MQ=39)
        tGTACAAAAAATTTGGGGTTATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGaact                                                                                      >  1:86991/1‑95 (MQ=25)
        tGTACAAAAAATTCGAGGTTATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGa                                                                                         <  1:1142623/94‑1 (MQ=255)
        tGTACAAAAAATTCGAGGTTATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGa                                                                                         >  2:1142622/1‑94 (MQ=9)
          tACAAAAAATTTGGGGTTATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAatat                                                                                    >  2:870288/1‑97 (MQ=255)
               aaaaTTTGGGGTTATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAgggg                                                                               >  1:1690523/1‑97 (MQ=255)
               aaaaTTTGGGGTTATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAgggg                                                                               >  1:1690529/1‑97 (MQ=255)
                aaaTTTGGGGTTATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGga                                                                              <  2:473496/97‑1 (MQ=255)
                aaaTTTGGGGTTATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGga                                                                              <  2:1072097/97‑1 (MQ=255)
                     tGGGGTTATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATaa                                                                         <  1:333184/97‑1 (MQ=255)
                     tGGGGTTATACCAGAAATGGAATATAAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATaa                                                                         <  2:2401430/97‑1 (MQ=255)
                     cgaggTTATACCAGAAATGGAATATGAAATGGAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATaa                                                                         >  2:2147832/4‑97 (MQ=12)
                      ggggTTATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAAc                                                                        >  1:2159982/1‑97 (MQ=255)
                         gTTATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATg                                                                     <  2:1736899/97‑1 (MQ=255)
                         gTTATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATg                                                                     <  1:954097/97‑1 (MQ=255)
                          ttATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATgg                                                                    <  1:933158/97‑1 (MQ=255)
                           tataCCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGc                                                                   >  1:1837892/1‑97 (MQ=255)
                           tataCCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGc                                                                   >  2:2564371/1‑97 (MQ=255)
                            ataCCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCt                                                                  >  1:588340/1‑97 (MQ=255)
                             taCCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCtt                                                                 <  1:173538/97‑1 (MQ=255)
                             aaCCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCtt                                                                 <  1:84087/96‑1 (MQ=255)
                                 aGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCga                                                             >  1:193021/1‑97 (MQ=255)
                                 aGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCga                                                             >  1:4998/1‑97 (MQ=255)
                                    aaTGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCgaaga                                                          <  2:508338/97‑1 (MQ=255)
                                     aTGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCgaagaa                                                         >  1:2074653/1‑97 (MQ=255)
                                          atatGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACt                                                    >  1:2444910/1‑97 (MQ=255)
                                           tatGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGc                                                                   <  1:2107553/81‑1 (MQ=255)
                                           tatGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTa                                                   <  2:957894/97‑1 (MQ=255)
                                           tatGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTa                                                   <  1:1332083/97‑1 (MQ=255)
                                           aaggaaatgaaaGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTa                                                   <  1:1900118/94‑1 (MQ=255)
                                             tgaaatgaaaGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTaaa                                                 >  2:1460428/1‑97 (MQ=255)
                                             tgaaatgaaaGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTaaa                                                 >  1:617060/1‑97 (MQ=255)
                                               aaatgaaaGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTaaatg                                               >  2:1171234/1‑97 (MQ=255)
                                               aaatgaaaGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTaaatg                                               >  2:531217/1‑97 (MQ=255)
                                               aaatgaaaGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTaaatg                                               >  1:1780334/1‑97 (MQ=255)
                                               aaatgaaaGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTaaatg                                               <  1:2197640/97‑1 (MQ=255)
                                               aaatgaaaGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTaaatg                                               >  2:1590537/1‑97 (MQ=255)
                                                aatgaaaGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTaaatga                                              <  1:326087/97‑1 (MQ=255)
                                                aatgaaaGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTaaatga                                              <  1:937619/97‑1 (MQ=255)
                                                  tggaagaaaTTAAAGAAATCGAACAATATATTAAAGAGCAAGAAGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGaaaa                                                       >  2:200215/4‑86 (MQ=0)
                                                  tgaaaGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGAAGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGaaaa                                                       <  1:200215/86‑1 (MQ=255)
                                                   gaaagaaaTTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTaaatgaaat                                           >  2:892904/1‑97 (MQ=255)
                                                    aaagaaaTTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTaaatgaaatg                                          <  2:791253/97‑1 (MQ=255)
                                                    aaagaaaTTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTaaatgaaatg                                          <  2:61774/97‑1 (MQ=255)
                                                    aaagaaaTTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTaaatgaaatg                                          >  1:1510010/1‑97 (MQ=255)
                                                    aaagaaaTTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTaaatgaaatg                                          >  2:1782327/1‑97 (MQ=255)
                                                     aagaaaTTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGt                                         >  1:1954600/1‑97 (MQ=255)
                                                      agaaaTTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTa                                        >  1:298436/1‑97 (MQ=40)
                                                       gaaaTTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCt                                                                  >  2:44352/1‑70 (MQ=255)
                                                       gaaaTTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCt                                                                  <  1:44352/70‑1 (MQ=255)
                                                       gaaaTTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTAc                                       <  1:1011490/97‑1 (MQ=255)
                                                         aaTTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACaa                                     >  2:266788/1‑97 (MQ=255)
                                                          aTTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACt                                                    <  2:1877199/81‑1 (MQ=255)
                                                          aTTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAAt                                    <  2:1923338/97‑1 (MQ=255)
                                                          aTTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAAt                                    >  2:1252610/1‑97 (MQ=255)
                                                          aTTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAAt                                    >  2:1465868/1‑97 (MQ=255)
                                                           ttAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATg                                   >  2:1559166/1‑97 (MQ=255)
                                                            tAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATga                                  >  1:595836/1‑97 (MQ=255)
                                                            tAAAGAAATCGAACAATATACTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTCTCGAAGAAAAACTAAATGAAATGTACCATga                                  >  1:1159059/1‑97 (MQ=255)
                                                               aGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGAt                               <  2:2261352/97‑1 (MQ=255)
                                                               aGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGAt                               <  2:2261337/97‑1 (MQ=255)
                                                               aGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGAt                               <  2:4998/97‑1 (MQ=255)
                                                                gAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGAtt                              <  2:1801891/97‑1 (MQ=255)
                                                                 aaaTCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGaaa                                                        <  2:744141/70‑1 (MQ=255)
                                                                  aaTCGAACAATATATTAAAGAGCAAGAAGAAGCTGAAATATAGGGGAGATAACATGGCt                                                                  <  1:563076/59‑1 (MQ=255)
                                                                    tCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGa                          >  2:1877347/1‑97 (MQ=255)
                                                                     cGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTATAGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGaa                         <  2:1121692/97‑1 (MQ=255)
                                                                      gAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAAt                        >  2:253277/1‑97 (MQ=255)
                                                                        aCAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATaa                      <  2:2398707/97‑1 (MQ=255)
                                                                         cAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATaaa                     >  1:2467459/1‑97 (MQ=255)
                                                                         cAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATaaa                     >  1:2632588/1‑97 (MQ=255)
                                                                          aaTATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATaaaa                    <  1:370750/97‑1 (MQ=255)
                                                                          aaTATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATaaaa                    >  2:141605/1‑97 (MQ=255)
                                                                            tatatTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAAtt                  <  1:1557711/97‑1 (MQ=255)
                                                                             atatTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAACGAGATTGCGAATAAGATTa                 >  2:1632751/1‑97 (MQ=255)
                                                                               atTAGAGAGCAAGATGAAGCTGAAAGATAGGGGAGATAACATGGCGTTCGAAGAAAAACTAAATGAAATGTATAATGAGATTGCGAATAAAATTAGc               <  2:1624966/97‑1 (MQ=17)
                                                                               atTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGc               >  1:712663/1‑97 (MQ=255)
                                                                                 tAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGCAg             <  1:887056/97‑1 (MQ=255)
                                                                                 tAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGCAg             <  1:2195145/97‑1 (MQ=255)
                                                                                  aaaGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGCAGt            >  2:190659/1‑97 (MQ=255)
                                                                                  aaaGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGCAGt            >  2:2040493/1‑97 (MQ=255)
                                                                                  aaaGAGCAAGAAGAAGCTGAAATATAGGGGAGATGACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGCAgc            >  2:1963706/1‑96 (MQ=255)
                                                                                  aaaGAGCAAGAAGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAATGAAATGAACTGTACCCTGAGATTGCGAATAAAACTAGCAGt            >  2:1247103/1‑97 (MQ=14)
                                                                                    agagCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGCAGTAt          <  2:1031413/97‑1 (MQ=255)
                                                                                     gagCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGCAGTATg         >  1:2633592/1‑97 (MQ=255)
                                                                                     gagCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGCAGTATg         >  1:1401254/1‑97 (MQ=255)
                                                                                     gagCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGCAGTATg         >  2:2193926/1‑97 (MQ=255)
                                                                                     gagCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGCAGTATg         >  1:2110192/1‑97 (MQ=255)
                                                                                     gagCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTNCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGCAGTATg         >  1:22087/1‑97 (MQ=255)
                                                                                      agCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGCAGTATGa        >  1:1053004/1‑97 (MQ=255)
                                                                                       gCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGCAGTATGAt       <  1:2147833/97‑1 (MQ=255)
                                                                                         aaGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGCAGTATGATAc     >  2:74369/1‑97 (MQ=255)
                                                                                          aGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAGTGAAATGTACAATGAGATTGCGAATAAGATTAGCAGTATGATAcc    <  2:481310/97‑1 (MQ=255)
                                                                                           gATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGCAGTATGATACCg   <  1:609304/97‑1 (MQ=255)
                                                                                            aTGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGCAGTATGATACCgg  <  1:2535093/97‑1 (MQ=255)
                                                                                            aTGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGCAGTATGATACCgg  <  1:1883065/97‑1 (MQ=255)
                                                                                             |                                                                                               
TTATTATATGTACAAAAAATTTGGGGTTATACCAGAAATGGAATATGAAATGAAAGAAATTAAAGAAATCGAACAATATATTAAAGAGCAAGATGAAGCTGAAATATAGGGGAGATAACATGGCTTTCGAAGAAAAACTAAATGAAATGTACAATGAGATTGCGAATAAAATTAGCAGTATGATACCGG  >  BX571857/321616‑321804

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