Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 52,934 G→A 5.5% intergenic (+544/‑129) USA300HOU_RS00215 → / → USA300HOU_RS00220 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE52,9340GA 5.5% 402.9 / 8.8 133intergenic (+544/‑129)USA300HOU_RS00215/USA300HOU_RS00220hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (66/59);  new base A (2/6);  total (68/65)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 136 total aligned reads displayed.

TCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGTTTGAGTAAATTA  >  USA300TCH1516_ALE/52854‑53024
                                                                                |                                                                                          
tcatcGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAaataa                                                                                  >  1:3071637/1‑91 (MQ=31)
 catcGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAaataat                                                                                 >  2:2018186/1‑91 (MQ=32)
  atcGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAaataata                                                                                >  1:2378727/1‑90 (MQ=17)
  atcGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAaataata                                                                                >  2:2638966/1‑90 (MQ=17)
  atcGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAaataata                                                                                >  2:5768427/1‑90 (MQ=17)
   tcGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCa                                                                               >  1:3635064/1‑91 (MQ=34)
   tcGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCa                                                                               >  1:2500748/1‑91 (MQ=34)
   tcGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCa                                                                               >  1:97214/1‑91 (MQ=34)
    cGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCag                                                                              <  1:2074495/91‑1 (MQ=33)
     gtagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCaga                                                                             >  1:1190008/1‑91 (MQ=34)
     gtagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCaga                                                                             >  1:226701/1‑91 (MQ=34)
     gtagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCaga                                                                             >  1:4188140/1‑91 (MQ=34)
      tagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCaggg                                                                            >  2:4202714/1‑89 (MQ=21)
      tagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCagag                                                                            >  1:1986732/1‑91 (MQ=33)
      tagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCagag                                                                            <  1:3448925/91‑1 (MQ=34)
      tagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGAAATCGCGCGAAAATCGCGTGAAATAATAagag                                                                            >  2:4202713/1‑91 (MQ=16)
       agtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAgg                                                                           >  2:4486556/1‑91 (MQ=34)
       agtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAgg                                                                           >  1:488750/1‑91 (MQ=34)
       agtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGAGAAAATCGCGTGAAATAATAAGAgg                                                                           >  2:47361/1‑91 (MQ=255)
          aaTCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAAt                                                                        >  1:43350/1‑91 (MQ=34)
          aaTCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAAt                                                                        <  1:2594066/91‑1 (MQ=34)
           agcgcgtgAAAATGGCGTGATATAATCAGAGGAATGGTTGGAAGTCATCGTATGAATGTCGCGATAATGGCGTGAGATAATCAGATGAAtg                                                                       >  2:3826567/3‑90 (MQ=2)
           aTCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATc                                                                       >  2:648955/1‑91 (MQ=34)
             cgcgcgAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGt                                                                     >  1:4441959/1‑91 (MQ=34)
               cgcgAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGttt                                                                   <  1:776464/91‑1 (MQ=33)
               cgcgAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGttt                                                                   >  1:4108584/1‑91 (MQ=34)
                gcgAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTg                                                                  >  1:122766/1‑91 (MQ=34)
                gcgAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGAGATGGTTTg                                                                  >  2:4543304/1‑91 (MQ=255)
                 ggAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAAGCGCGTGAAATAATCAGAGGAATCGTTTGa                                                                 <  2:4491623/90‑1 (MQ=18)
                 cgAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGa                                                                 <  1:5071293/91‑1 (MQ=33)
                  gAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGaa                                                                <  1:3184400/91‑1 (MQ=33)
                  gAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGaa                                                                <  1:5317450/91‑1 (MQ=33)
                   aaaaTCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGaaa                                                               <  1:5478837/91‑1 (MQ=34)
                    aaaTCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAat                                                              >  1:4530303/1‑91 (MQ=33)
                    aaaTCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAat                                                              >  1:5286761/1‑91 (MQ=33)
                     aaTCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAatc                                                             <  1:5621897/91‑1 (MQ=34)
                      aTCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAatca                                                            <  1:5046469/91‑1 (MQ=33)
                       tCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAatcat                                                           >  2:2157229/1‑91 (MQ=33)
                       tCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAatcat                                                           >  1:3243045/1‑91 (MQ=33)
                         gcgTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCg                                                         >  2:4449486/1‑91 (MQ=31)
                         gcgTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCg                                                         >  2:4433353/1‑91 (MQ=31)
                                                         atcGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACtc                         >  2:3168405/1‑90 (MQ=17)
                                                         atcGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTa                         >  2:782966/1‑91 (MQ=33)
                                                         atcGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTa                         <  2:741276/91‑1 (MQ=33)
                                                         atcGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTa                         <  2:2475596/91‑1 (MQ=33)
                                                         atcGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTa                         >  1:4683889/1‑91 (MQ=33)
                                                         atcGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTa                         <  1:5066905/91‑1 (MQ=33)
                                                         atcGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTa                         >  2:3400894/1‑91 (MQ=33)
                                                           cGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGTGTGAAGTACCATg                       >  2:4343345/1‑91 (MQ=255)
                                                           cGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATg                       >  2:1883510/1‑91 (MQ=35)
                                                           cGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATg                       >  1:1628881/1‑91 (MQ=35)
                                                           cGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGGGAAATACTATg                       >  2:4537304/1‑91 (MQ=25)
                                                           cGTAGTAATCGCGCGAAAATCACGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATg                       >  2:4097697/1‑91 (MQ=17)
                                                            gtagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATgg                      >  1:2359495/1‑91 (MQ=37)
                                                            gtagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATgg                      <  1:3931651/91‑1 (MQ=37)
                                                             tagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGt                     <  2:3455723/91‑1 (MQ=38)
                                                             tagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGt                     <  1:2463168/91‑1 (MQ=38)
                                                             tagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGt                     <  2:3017473/91‑1 (MQ=38)
                                                             tagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGt                     >  2:2790335/1‑91 (MQ=38)
                                                             tagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGt                     <  2:2187353/91‑1 (MQ=38)
                                                             tagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGt                     <  2:787253/91‑1 (MQ=38)
                                                             tagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGt                     <  2:2365901/91‑1 (MQ=38)
                                                             tagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGt                     <  2:1574509/91‑1 (MQ=38)
                                                             tagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGt                     <  1:4291727/91‑1 (MQ=38)
                                                             tagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGt                     <  1:3851816/91‑1 (MQ=38)
                                                             tagtaATCGCGCGAAAATCGCGTGAAATAATCAGAGAAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATggg                     >  2:5643765/1‑90 (MQ=255)
                                                              agtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAATTCGCGTGAAATACTATGGTa                    >  2:2583890/1‑91 (MQ=255)
                                                              agtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTa                    <  2:1031638/91‑1 (MQ=39)
                                                              agtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTa                    >  1:370765/1‑91 (MQ=39)
                                                              agtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTa                    <  1:2898430/91‑1 (MQ=39)
                                                               gtgaTCTCGCGAAAATCGTGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAg                   <  1:4904523/88‑1 (MQ=255)
                                                               gtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGATATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTACGGTAg                   <  2:5750962/91‑1 (MQ=255)
                                                               gtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAg                   >  1:5679527/1‑91 (MQ=40)
                                                               gtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAg                   <  1:4488732/91‑1 (MQ=40)
                                                               gtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAg                   <  2:1024283/91‑1 (MQ=40)
                                                               gtaATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATAGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAg                   <  1:4860656/91‑1 (MQ=255)
                                                               gtaATCGCGCGAAAATCACGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATgaggg                   >  1:1748825/1‑87 (MQ=9)
                                                                taATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGa                  >  1:3782987/1‑91 (MQ=255)
                                                                   tCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACgg               >  1:441842/1‑91 (MQ=255)
                                                                   tCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAAACGCGTGAAATACTATGGTAGACgg               >  2:729008/1‑91 (MQ=255)
                                                                    cgcgcgAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGt              <  2:375348/91‑1 (MQ=255)
                                                                     gcgcgAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGtt             <  1:4299908/91‑1 (MQ=255)
                                                                      cgcgAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGttt            >  1:123340/1‑91 (MQ=255)
                                                                      cgcgAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGttt            >  1:1799255/1‑91 (MQ=255)
                                                                      cgcgAAAATCACGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGttt            <  2:2369797/91‑1 (MQ=255)
                                                                        cgAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGTTTGa          <  1:5200961/91‑1 (MQ=255)
                                                                         gAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGTTTGAg         <  1:2944821/91‑1 (MQ=255)
                                                                         gAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGTTTGAg         <  1:2383878/91‑1 (MQ=255)
                                                                         gAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGTTTGAg         <  1:5416446/91‑1 (MQ=255)
                                                                         gAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGTTTGAg         <  1:5768427/91‑1 (MQ=255)
                                                                         gAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAAGCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGTTTGAg         <  1:2672149/91‑1 (MQ=255)
                                                                         caaaaTCGCGTGAAATAATCAGAGGAATCGTTTGAACTCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTAAGGTAGACGGTTTGAg         <  1:3826567/90‑1 (MQ=255)
                                                                          aaaaTCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGTTTGAGt        <  1:648955/91‑1 (MQ=255)
                                                                          aaaaTCACGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGTTTGAGt        <  2:4204858/91‑1 (MQ=255)
                                                                           aaaTCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGTTTGAGTa       >  1:3102533/1‑91 (MQ=255)
                                                                            aaTCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGTTTGAGTaa      <  1:5100678/91‑1 (MQ=255)
                                                                              tCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGTTTGAGTAAAt    >  1:2684411/1‑91 (MQ=255)
                                                                              tCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGTTTGAGTAAAt    >  2:3078411/1‑91 (MQ=255)
                                                                                gcgTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGTTTGAGTAAATTa  >  1:3751996/1‑91 (MQ=255)
                                                                                gcgTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGTTTGAGTAAATTa  >  1:4225623/1‑91 (MQ=255)
                                                                                |                                                                                          
TCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATAATCAGAGGAATCGTTTGAAATCATCGTAGTAATCGCGCGAAAATCGCGTGAAATACTATGGTAGACGGTTTGAGTAAATTA  >  USA300TCH1516_ALE/52854‑53024

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