Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,309,899 G→T 12.7% S322* (TCA→TAA)  yfmC_2 ← Fe(3+)‑citrate‑binding protein YfmC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,309,8990GT12.7% 143.1 / 2.6 62S322* (TCA→TAA) yfmC_2Fe(3+)‑citrate‑binding protein YfmC
Reads supporting (aligned to +/- strand):  ref base G (34/18);  new base T (7/2);  total (41/21)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTGCCCAAACATCACGGTCAACAATATCCACGCGATTATTTTTAACTGCATTCAACTTTTTCCATGTTGCATCTT  >  USA300TCH1516_ALE/2309760‑2310028
                                                                                                                                           |                                                                                                                                 
aaTATATACAAATAAAAAGAAAGCACAACCTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTg                                                                                                                                   <  1:628098/140‑1 (MQ=255)
aaTATATACAAATAAAAAGAAAGCACAAACCGACTATATACCTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTg                                                                                                                                   <  1:446283/140‑1 (MQ=255)
             aaaaaGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAg                                                                                                                      >  1:259794/1‑140 (MQ=255)
               aaagaacgCACAAACACACTATATAGGTGAGTGTTGTGCGTCTTTTTCTTTTATTATCTATGGCAGATTGTCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTTTTCTTTTTTTTATAATTCAACTAGtt                                                                                                                    >  1:334042/1‑140 (MQ=255)
               aaagaaagCACAAACTCACTCTATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGCTCTATGGCATATTGGCTACTTATTTCTTTTATAGCCATTTACTTTCATCTTACTTTTGTTGTTTTTTTNATAATTCAACAAAtt                                                                                                                    >  1:498220/1‑140 (MQ=255)
               aaagaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTCATTTTGTTCTTTTTTTGGTATTTCAACAAGtt                                                                                                                    >  2:384976/1‑140 (MQ=255)
                 agaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTCTGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTATTTCCACTTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCt                                                                                                                  >  1:152146/1‑140 (MQ=255)
                 agaaagCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCAACTTACTTTTGTTCTTTTTTTGATAATTCAACAGGTTcg                                                                                                                  >  2:51663/1‑139 (MQ=255)
                       cacaAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGCTTGGCTACTTCTTTCTTTTATAGCCATTTACTTCCCCCTTACTTTTTTTCTTTTTTTGTTATTTCAACAAGTTCTTTATcc                                                                                                            >  1:485889/1‑140 (MQ=255)
                       cacaAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGcc                                                                                                            <  1:300756/140‑1 (MQ=255)
                         caAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCCCCTTACTTTTGTTCTTTTTTTTATAATTCACCAAGTTTTTTAGCCAt                                                                                                          >  1:199373/1‑140 (MQ=255)
                           aaCTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTACATAATTTGATAATTCAACAAGTCCTTTAGCCAttt                                                                                                        >  1:285664/1‑140 (MQ=255)
                                aCTATATACGTGAGTGTTGTGCGTCTTTTTCATTTATGATATGTGTCTGTTTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCa                                                                                                   <  2:301677/140‑1 (MQ=255)
                                   atataCGTGGGTGTTGTGGGTCTTTTTCTTTTATGATCTATGGCAGGTTGGCTAGTTATTTCTTTTATAGCGATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTTCACGAGTTCTTTAGCCAGTTCTTTagaa                                                                                                >  2:249397/1‑140 (MQ=255)
                                   atataCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCagaa                                                                                                <  2:77273/140‑1 (MQ=255)
                                     ataCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCCTTTACTTCCACCTTCCTTTTGTTCTTTTTTTGTTAATTTAACAAGTTCTTTTGCCATTTCTTCataaga                                                                                              >  2:20573/1‑140 (MQ=255)
                                       aCGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGTCTCCTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGTTAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGaaa                                                                                            >  1:334435/1‑140 (MQ=255)
                                           gAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTTCTTATTTCTTTTTTAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTaa                                                                                        <  1:145921/140‑1 (MQ=255)
                                             gtgtTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCAATTTCTTCAGAAGAAATTAAGc                                                                                      >  1:508872/1‑140 (MQ=255)
                                              tgttgaGCGTCATTTACTTTTATAATCTCTGGCAGATTGGCTACGTATTTCTTTTAAAGCCATTTACTTCCACCTACCTTTTGTTCTTTTTTTGATAATTCAACAAGTTTTTTAGCCATTTCTACAGAAGAAAATacaca                                                                                     >  2:289454/1‑136 (MQ=255)
                                               gttgtGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCa                                                                                    >  2:588304/1‑140 (MQ=255)
                                                 tgtgCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTTTTTATTCTACAAGTTCTTTAGCCATTTCTTCAGAAGAAAATAAGCCACg                                                                                  >  2:416360/1‑140 (MQ=255)
                                                 tgtgCGTCTTTTTCTTTTAAGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACg                                                                                  >  1:399004/1‑140 (MQ=255)
                                                  gtgCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACga                                                                                 <  1:515057/140‑1 (MQ=255)
                                                   tgCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACgag                                                                                <  2:666547/140‑1 (MQ=255)
                                                    gCGTCTTTTTCATTTATGAACCATGGCAGATTGGATCCTTGTTTCTGTTATCGACATATACTTACCCCGTACTTTTGTTCTTTTTTTGATAATTCCACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACgcga                                                                               >  2:169637/1‑140 (MQ=255)
                                                      gTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTATTCTTTTTTTGATAATTCAGCAAGTTCTTTGTCCATTTCTTCAGAAGCAATTAAGCCACGAGATc                                                                             >  2:542722/1‑140 (MQ=255)
                                                       tctttttcttttATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCt                                                                            <  2:508472/140‑1 (MQ=255)
                                                       tccttttcttttATGATCTATGGCAGCTTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCt                                                                            <  2:444002/137‑1 (MQ=255)
                                                        ctttttcttttATGATCTATGGCAGATTAGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGTTAATTCAACAATTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACAAGATCtt                                                                           >  1:456161/1‑140 (MQ=255)
                                                         tttttcttttATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTg                                                                          >  2:287025/1‑140 (MQ=255)
                                                               ttttATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACCTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCATAAGCACTTAAGCCCAGAGATCTTGGCCaaa                                                                    >  2:260039/1‑140 (MQ=255)
                                                               ttttATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACCTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGCGATCTTGCCCaaa                                                                    >  2:168356/1‑140 (MQ=255)
                                                                  ttttctCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTTTACAATTTCTTTAGCCATTTCTTCAGAAGATATTAAGCCCCGAAATCTTGCGCAAACAt                                                                 >  2:461506/6‑140 (MQ=255)
                                                                   aTGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTGCCCAAACATc                                                                >  1:351530/1‑140 (MQ=255)
                                                                      aTCTATTGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTGCCCAAACATCACg                                                             <  2:336732/140‑1 (MQ=255)
                                                                         tATGGCAAATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCCTTTCTTCAGGAGAAATTAAGCAACGAGATCTTGCCCAAACATCACGGTc                                                          >  2:468852/1‑140 (MQ=255)
                                                                          aTGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTGCCCAAACATCACGGTca                                                         >  2:226652/1‑140 (MQ=255)
                                                                          aTGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCAATTTCTTCAGAAGAAATTAAGCCCCGAGATCTTGCCCAAACATCACgttat                                                         >  2:84529/1‑136 (MQ=255)
                                                                              cAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTAAGAAGAAATTAAGCCACGAGATCTTACCCAAACATCACGGTcaacaa                                                     >  2:524941/1‑140 (MQ=255)
                                                                                 aTTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTGCCCAAACATCACGGTCAACAattt                                                  >  1:507456/1‑138 (MQ=255)
                                                                                  tgggCTACTTATTTCTTTTGTAGCCATTTTCTTCCTCCTCTCTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTTAGAAGAAATTAAGCCACGAGATCTTGCCCAAACATCACGGTCAACAATATc                                                 <  2:424602/138‑1 (MQ=255)
                                                                                      cTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTGCCCAAACATCAAGGTCACCAAAATCCCcg                                             >  2:442637/1‑140 (MQ=255)
                                                                                       tACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTGCCCAAACATCACGGTCAACAATATCCAcgc                                            <  1:469846/140‑1 (MQ=255)
                                                                                                    tatcGCCATTTATTTCAACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGATGAAATTAAGCCACGAGATCTTGCCCAAACATCACGGTCAACAATATCCACGCGATTATTTTTAAc                               <  1:91473/140‑1 (MQ=255)
                                                                                                     ataGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTGCCCAAACATCACGGTCAACAATATCCACGCGATTATTTTTAACt                              <  2:663336/140‑1 (MQ=255)
                                                                                                       aGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTTATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTGCCCATACATCACGGTCAACAATCTCCATCCGACTATTTTTAACtcc                            >  1:321840/1‑138 (MQ=255)
                                                                                                       aGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTGCCCAAACATCACGGTCAACAATATCCACGCGATTATTTTTAACTGc                            >  2:33884/1‑140 (MQ=255)
                                                                                                         ccATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTGCCCAAACATCACGGTCAACAATATCCACGCGATTATTTTTAACTTCAt                          >  1:628084/1‑140 (MQ=255)
                                                                                                           aTTTACTTCCACCTTACTTTTGTTCGTTTTTTGATAATTCATCAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACTAGATCTTGCCCAAACATCACGGTCAACAATATCCACGGGATTATTTTTAACTGCATTc                        <  1:406445/140‑1 (MQ=255)
                                                                                                              tACTTCCACCTTACTTTTGTTCTTTTTTTTATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTTCCCAAACATCACGGTCAACAATATCCACGCGATTATTTTTAACTGCATTCAAc                     >  2:206380/1‑140 (MQ=255)
                                                                                                              tACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCGTACCCAAACATCACGGTCAACAATATCCACGCGATTATTTTTAACTGCATTCAAc                     >  2:401579/1‑140 (MQ=255)
                                                                                                                cgtCCAGCTTACTTTTGTGCTTTTTTTAATAGTTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTGCCCAAACATCACGGTCAACAATATCCACGCGATTATTTTTAACTGCATTCAACtt                   <  2:597747/138‑1 (MQ=255)
                                                                                                                cTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTGCCCAAACATCACGGTAAACAATATCCTCGCGATTTTTTTTTACTGCATTCAACtt                   >  1:436081/1‑140 (MQ=255)
                                                                                                                cTTCCACCTTACTTTTGCTCATTGCTTTAGAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTGCCCAAACATCACGGTCAACAATATCCACGCGATTATTTTTCACTGCACTCCACtt                   >  2:642970/1‑140 (MQ=255)
                                                                                                                    cACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGCAATTAAGCCACGAGATCTTGCCCAAACATCACGGTCAACAATATTCACGCGATTATTTTTAACTTCATTTAACTTTTTc               >  1:286937/1‑140 (MQ=255)
                                                                                                                         tACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTGCCCAAACATCACGGTCAACAATATCCACGCGATTATTTTTAACTGCATTCAACTTTTTCCATGt          <  1:91510/140‑1 (MQ=255)
                                                                                                                          aCTTTTGTTCTTTTTTTTTTTTTTTAACAAGATCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTGCCCAAACATCACGGTCAACAATATCCACGCGATTATTTTTAACTGCATTCAACTTTTTCCATGtt         >  2:613140/1‑140 (MQ=255)
                                                                                                                           cTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTGCCCAAACATCACGGTCAACAATATCCACGCGATTATTTTTAACTGCATTCAACTTTTTCCATGTTg        <  1:606219/140‑1 (MQ=255)
                                                                                                                             tttGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAGTTAAGCCACGAGATGTTGCCCAAACATCACGGTCAACAATATACACGCGATTATTTTTAACGGCATTCACCTTTTTCTATGTTGct      >  2:108438/1‑139 (MQ=255)
                                                                                                                               tGTTCTTTTTTTTATAATTCAACAAGTTCTTTAGCCATCTCTTCAGAATAAATTAAGCCACGAGATCTTGCCCAAACATCACGGTCCACTATATCCACGCGATTATTTTTAACTGCATTCAACTTTTTCCATGTTGCAtc    <  1:321409/140‑1 (MQ=255)
                                                                                                                               tGTTCTTTTTTTTATAATTCAACAAGTTCTTTAGCCATCTCTTCAGAAGAACTTAAGCCACGAGACCTTGCCCACACCTCACGGTCAACAATATCCACTCGATTATTTTTAACTGCATTCAACTTTTTCCATGTTGCAtc    <  2:355204/140‑1 (MQ=255)
                                                                                                                               tGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGATATCTTGCTCAATCATCATGGTTAACAATATCCACGCGATTATTTTTAACTGCATTCAACTTTTTCCATGTTGCAtc    <  1:61402/140‑1 (MQ=255)
                                                                                                                                gTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTGCCCAAACATCACGGTCAACAATATCCACGCGATTATTTTTAACTGCATTCAACTTTTTCCATGTTGCAtct   >  1:409659/1‑140 (MQ=255)
                                                                                                                                 ttCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGTAGAAATTAAGCCACGAGATCTTGCCCAAACATCACGGTCAACAATATCCACGCGATTATTTTTAACTGCATTCACCTTTTTCCATGTTGCAtctt  <  1:467986/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                 
AATATATACAAATAAAAAGAAAGCACAAACTCACTATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTTCTTTTTTTGATAATTCAACAAGTTCTTTAGCCATTTCTTCAGAAGAAATTAAGCCACGAGATCTTGCCCAAACATCACGGTCAACAATATCCACGCGATTATTTTTAACTGCATTCAACTTTTTCCATGTTGCATCTT  >  USA300TCH1516_ALE/2309760‑2310028

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