Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 396,540 G→A 41.7% A208T (GCA→ACA)  efeM → putative iron uptake system component EfeM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE396,5400GA41.7% 48.6 / 123.3 99A208T (GCA→ACA) efeMputative iron uptake system component EfeM
Reads supporting (aligned to +/- strand):  ref base G (24/33);  new base A (20/21);  total (44/55)
Fisher's exact test for biased strand distribution p-value = 5.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-01

Only 100 of 102 total aligned reads displayed.

AATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTG  >  USA300TCH1516_ALE/396453‑396629
                                                                                       |                                                                                         
aaTGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCacac                                                                                        <  1:2111904/91‑1 (MQ=255)
aaTGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGcac                                                                                        <  1:2613630/91‑1 (MQ=255)
 aTGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGcacaca                                                                                       <  1:3170728/91‑1 (MQ=255)
 aTGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGcaca                                                                                       <  1:606371/91‑1 (MQ=255)
 aTGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGcaca                                                                                       <  1:606308/91‑1 (MQ=255)
    aaGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACaaaa                                                                                    >  1:2571571/1‑91 (MQ=255)
    aaGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACaaaa                                                                                    <  1:3189925/91‑1 (MQ=255)
     aGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACaaaaa                                                                                   >  1:1858129/1‑91 (MQ=255)
       ttGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAAt                                                                                 <  1:3246248/91‑1 (MQ=255)
         gCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAAttt                                                                               >  1:3147569/1‑91 (MQ=255)
         gCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAAttt                                                                               >  1:3423443/1‑91 (MQ=255)
          cAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTa                                                                              <  1:3788508/91‑1 (MQ=255)
            aCTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATg                                                                            >  1:294361/1‑91 (MQ=255)
              ttCTAAACTCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGAc                                                                          <  1:12206/91‑1 (MQ=255)
              ttCTAAAATCACAGGTGAAGAAGAAATTTATTTACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGAc                                                                          <  1:3524049/91‑1 (MQ=255)
              ttCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGAc                                                                          >  1:1051965/1‑91 (MQ=255)
              ttCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGAc                                                                          >  1:3090602/1‑91 (MQ=255)
                cTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACtt                                                                        >  1:3033010/1‑91 (MQ=255)
                cTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACtt                                                                        >  1:1980746/1‑91 (MQ=255)
                   aaaTCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTAtt                                                                     <  1:3829559/91‑1 (MQ=255)
                    aaTCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTAttt                                                                    <  1:4095502/91‑1 (MQ=255)
                     aTCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTa                                                                   <  1:890837/91‑1 (MQ=255)
                      tCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTaa                                                                  >  1:3433648/1‑91 (MQ=255)
                      tCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTaa                                                                  <  1:1496297/91‑1 (MQ=255)
                      tCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTaa                                                                  <  1:3148953/91‑1 (MQ=255)
                       cacaGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTaaa                                                                 <  1:811713/91‑1 (MQ=255)
                       cacaGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTaaa                                                                 <  1:1207633/91‑1 (MQ=255)
                       cacaGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTaaa                                                                 <  1:1207675/91‑1 (MQ=255)
                        acaGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAAc                                                                <  1:2283324/91‑1 (MQ=255)
                         caGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAAcc                                                               <  1:1707182/91‑1 (MQ=255)
                              gaagaagaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTAttt                                                          >  1:1442700/1‑91 (MQ=255)
                              gaagaagaaATTTATTAACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTAttt                                                          >  1:1235808/1‑91 (MQ=255)
                               aagaagaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTAtttt                                                         <  1:3335834/91‑1 (MQ=255)
                                agaagaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTa                                                        <  1:2447071/91‑1 (MQ=255)
                                agaagaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTa                                                        <  1:2166386/91‑1 (MQ=255)
                                  aagaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTaga                                                      >  1:246690/1‑91 (MQ=255)
                                   agaaATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTagag                                                     >  1:4101975/1‑91 (MQ=255)
                                     aaaTTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGCGaa                                                   >  1:3342942/1‑91 (MQ=255)
                                     aaaTTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaa                                                   >  1:3342682/1‑91 (MQ=255)
                                       aTTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaaaa                                                 <  1:2341395/91‑1 (MQ=255)
                                       aTTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaaaa                                                 >  1:3424026/1‑91 (MQ=255)
                                       aTTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaaaa                                                 <  1:169943/91‑1 (MQ=255)
                                        tttATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAATGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGaaaaa                                                <  1:3288312/91‑1 (MQ=255)
                                          tATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAg                                              >  1:3335911/1‑91 (MQ=255)
                                           aTTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGa                                             <  1:2272394/91‑1 (MQ=255)
                                           aTTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGa                                             <  1:2272429/91‑1 (MQ=255)
                                           aTTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGa                                             <  1:2272490/91‑1 (MQ=255)
                                           aTTCACATACAGATTTATATGATTTTAAAGCGAACGTGGAAGGCGCACAAATAATTTATGACTTATTTAAACCTATTTTAGAGAATAAAGa                                             >  1:3859640/1‑91 (MQ=255)
                                               acaTACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATaaa                                         <  1:438172/91‑1 (MQ=255)
                                                caTACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATaaaa                                        >  1:222094/1‑91 (MQ=255)
                                                 aTACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATaaaaa                                       >  1:3224850/1‑91 (MQ=255)
                                                 aTACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATaaaaa                                       <  1:3208944/91‑1 (MQ=255)
                                                    cAGATTTATATGATTTTAAAGCGAACGTTGAAGGCTCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAAtt                                    <  1:151547/91‑1 (MQ=255)
                                                    cAGATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAAtt                                    >  1:2607954/1‑91 (MQ=255)
                                                      gATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTaa                                  <  1:3086445/91‑1 (MQ=255)
                                                      gATTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTaa                                  <  1:3086461/91‑1 (MQ=255)
                                                       aTTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGCTAAAAAATTAAg                                 <  1:3176807/91‑1 (MQ=255)
                                                       aTTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAg                                 >  1:1472908/1‑91 (MQ=255)
                                                       aTTTATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAg                                 <  1:88491/91‑1 (MQ=255)
                                                        tttATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGt                                <  1:3717512/91‑1 (MQ=255)
                                                        tttATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGt                                <  1:3893431/91‑1 (MQ=255)
                                                        tttATATGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGt                                >  1:327508/1‑91 (MQ=255)
                                                         ttATATGATTTTAAAGCGAACGTTGCAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtg                               <  1:3849077/91‑1 (MQ=255)
                                                         ttATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtg                               <  1:3315981/91‑1 (MQ=255)
                                                         ttATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtg                               <  1:1479328/91‑1 (MQ=255)
                                                           atatGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtgat                             <  1:3690847/91‑1 (MQ=255)
                                                           atatGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtgat                             >  1:55532/1‑91 (MQ=255)
                                                            tatGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtgatg                            >  1:847750/1‑91 (MQ=255)
                                                            tatGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtgatg                            <  1:3179954/91‑1 (MQ=255)
                                                            tatGATTTTAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGtgatg                            >  1:848127/1‑91 (MQ=255)
                                                                aTTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTGAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGatat                        >  1:4136949/1‑91 (MQ=255)
                                                                 ttttAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATc                       >  1:62118/1‑91 (MQ=255)
                                                                  tttAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATcc                      <  1:835448/91‑1 (MQ=255)
                                                                   ttAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCa                     <  1:1524727/91‑1 (MQ=255)
                                                                    tAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCaa                    >  1:497332/1‑91 (MQ=255)
                                                                    tAAAGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCaa                    <  1:2414958/91‑1 (MQ=255)
                                                                     aaaGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCaaa                   <  1:1137987/91‑1 (MQ=255)
                                                                     aaaGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCaaa                   <  1:1138240/91‑1 (MQ=255)
                                                                     aaaGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCaaa                   <  1:3007957/91‑1 (MQ=255)
                                                                      aaGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAAt                  >  1:14049/1‑91 (MQ=255)
                                                                      aaGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAAt                  <  1:913372/91‑1 (MQ=255)
                                                                       aGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATg                 <  1:242138/91‑1 (MQ=255)
                                                                       aGCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATg                 >  1:1550551/1‑91 (MQ=255)
                                                                        gCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGa                <  1:1309737/91‑1 (MQ=255)
                                                                        gCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGa                >  1:1290854/1‑91 (MQ=255)
                                                                        gCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGa                >  1:3095283/1‑91 (MQ=255)
                                                                        gCGAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGa                >  1:416496/1‑91 (MQ=255)
                                                                         cGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGaa               >  1:2709643/1‑91 (MQ=255)
                                                                          gAACGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAAc              >  1:3476418/1‑91 (MQ=255)
                                                                           aaCGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACt             <  1:3580020/91‑1 (MQ=255)
                                                                             cGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTc           >  1:1048477/1‑91 (MQ=255)
                                                                             cGTTGAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTc           >  1:472749/1‑91 (MQ=255)
                                                                               ttGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGa         <  1:4084239/91‑1 (MQ=255)
                                                                                tGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGAt        >  1:1960652/1‑91 (MQ=255)
                                                                                 gAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATa       >  1:3415790/1‑91 (MQ=255)
                                                                                 gAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATa       >  1:2986012/1‑91 (MQ=255)
                                                                                 gAAGGCACACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATa       >  1:3798092/1‑91 (MQ=255)
                                                                                   aGGCGCACAAAAAATTTATGACTTATTTAAACGTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATaaa     <  1:3539007/91‑1 (MQ=255)
                                                                                   aGGCGCACAAAAAATTTATGACTTATTTAAACGTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATaaa     <  1:3539458/91‑1 (MQ=255)
                                                                                      cgcACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTg  >  1:2556994/1‑91 (MQ=255)
                                                                                       |                                                                                         
AATGAAGTTGCAACTTCTAAAATCACAGGTGAAGAAGAAATTTATTCACATACAGATTTATATGATTTTAAAGCGAACGTTGAAGGCGCACAAAAAATTTATGACTTATTTAAACCTATTTTAGAGAAAAAAGATAAAAAATTAAGTGATGATATCCAAATGAACTTCGATAAAGTG  >  USA300TCH1516_ALE/396453‑396629

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