Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,526,532 T→A 7.2% K471* (AAA→TAA)  dinG_1 ← putative ATP‑dependent helicase DinG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,526,5320TA 7.2% 795.3 / 36.7 236K471* (AAA→TAA) dinG_1putative ATP‑dependent helicase DinG
Reads supporting (aligned to +/- strand):  ref base T (104/115);  new base A (9/8);  total (113/123)
Fisher's exact test for biased strand distribution p-value = 8.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

Only 100 of 252 total aligned reads displayed.

ATTTTTCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAATCATCAAACAA  >  USA300TCH1516_ALE/1526442‑1526622
                                                                                          |                                                                                          
aTTTTTCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCTTTTTTGTTTTTACCAATCAAACCAAGCTGATAttt                                                                                            >  2:2293985/1‑91 (MQ=255)
aTTTTTCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATAttt                                                                                            >  1:3809942/1‑91 (MQ=255)
aTTTTTCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATAttt                                                                                            <  1:46485/91‑1 (MQ=255)
aTTTTTCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATAttt                                                                                            >  1:87365/1‑91 (MQ=255)
aTTTTTCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATAttt                                                                                            <  2:2325800/91‑1 (MQ=255)
aTTTTTCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATAttt                                                                                            <  1:2452739/91‑1 (MQ=255)
aTTTTTCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATAttt                                                                                            <  1:1795815/91‑1 (MQ=255)
 tttttCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTACATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTa                                                                                           <  1:2524605/91‑1 (MQ=255)
   tttCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTAAat                                                                                         <  2:1652723/91‑1 (MQ=255)
    tggaaaaaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAata                                                                                        >  1:2548108/5‑91 (MQ=255)
     tCCAAAATTCGCTGCTTTTCCGGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATAGTAAatat                                                                                       <  2:1316411/91‑1 (MQ=255)
     tCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAatat                                                                                       >  1:3756094/1‑91 (MQ=255)
     tCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTAAatat                                                                                       <  2:26880/91‑1 (MQ=255)
     ggaaaaaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAatat                                                                                       >  1:1533602/4‑91 (MQ=255)
      ccAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATc                                                                                      <  2:1248421/91‑1 (MQ=255)
      ccAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATc                                                                                      >  2:3522311/1‑91 (MQ=255)
       cAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCa                                                                                     >  1:778924/1‑91 (MQ=255)
       cAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTAAATATCa                                                                                     <  1:919880/91‑1 (MQ=255)
        aaaaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAg                                                                                    >  2:845940/1‑91 (MQ=255)
        aaaaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAg                                                                                    >  1:1123471/1‑91 (MQ=255)
          aaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCa                                                                                  <  2:3045588/91‑1 (MQ=255)
           aTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCAt                                                                                 <  2:706662/91‑1 (MQ=255)
              cGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaac                                                                              <  1:2850833/91‑1 (MQ=255)
              cGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTGTCGTTTTTGCCAATCAAACCAAGCTGATATTTAATATCAGCATaac                                                                              <  1:3091999/91‑1 (MQ=255)
               gctgctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaact                                                                             >  1:1979761/1‑91 (MQ=255)
                ctgctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAGTATCAGCATaactc                                                                            >  1:2266489/1‑91 (MQ=255)
                ctgctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactc                                                                            >  1:3475107/1‑91 (MQ=255)
                ctgctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactc                                                                            >  2:778162/1‑91 (MQ=255)
                ctgctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactc                                                                            >  1:1303276/1‑91 (MQ=255)
                 tgctTTTCCAGTTGATCAATTGCTTTTAACAATTTTTCATTTTCGTTTTTACTAATCAAACCAAGCTGCTATTTACTATCAGCATaactaa                                                                           >  2:1187413/1‑89 (MQ=255)
                 tgctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactca                                                                           <  2:1673894/91‑1 (MQ=255)
                  gctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactcaa                                                                          <  1:2011256/91‑1 (MQ=255)
                       tCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCAt                                                                     <  2:3182063/91‑1 (MQ=255)
                        ccAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCAtt                                                                    >  2:2375846/1‑91 (MQ=255)
                               aTCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACt                                                             <  2:1075324/91‑1 (MQ=255)
                                tCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACtt                                                            <  1:2417071/91‑1 (MQ=255)
                                 cAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTg                                                           <  2:2708109/91‑1 (MQ=255)
                                 cAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTg                                                           <  2:913574/91‑1 (MQ=255)
                                   aTTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGAt                                                         >  1:1780967/1‑91 (MQ=255)
                                    ttGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGAtt                                                        >  1:146358/1‑91 (MQ=255)
                                      gCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTa                                                      <  1:564030/91‑1 (MQ=255)
                                        cTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAg                                                    >  1:1419678/1‑91 (MQ=255)
                                         tttAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAgg                                                   <  1:3165045/91‑1 (MQ=255)
                                           tAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCa                                                 >  2:1017569/1‑91 (MQ=255)
                                           tAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCa                                                 >  2:1984633/1‑91 (MQ=255)
                                             aCAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATa                                               <  1:1270724/91‑1 (MQ=255)
                                              cAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAg                                              <  1:547997/91‑1 (MQ=255)
                                              cAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAg                                              <  2:2897285/91‑1 (MQ=255)
                                               aaTTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGt                                             >  1:1212884/1‑91 (MQ=255)
                                                aTTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTc                                            >  1:2898713/1‑91 (MQ=255)
                                                aTTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTc                                            >  1:2899591/1‑91 (MQ=255)
                                                aTTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTc                                            >  2:640264/1‑91 (MQ=255)
                                                aTTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTAACTTGATTTAAGGCATAGTc                                            >  1:1722237/1‑91 (MQ=255)
                                                 tttttCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCt                                           >  2:1874687/1‑91 (MQ=255)
                                                  ttttCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTg                                          >  1:39663/1‑91 (MQ=255)
                                                       aTTTTCGTTTTTACCAATCAAACCATGCTGATATTTAATATCAGCATAACGCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAg                                     <  1:3620122/91‑1 (MQ=255)
                                                        ttttCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGc                                    >  2:3877096/1‑91 (MQ=255)
                                                        ttttCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGc                                    <  1:48126/91‑1 (MQ=255)
                                                         tttCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCg                                   <  2:1742258/91‑1 (MQ=255)
                                                           tCGTTTTTACCAATCAAACCAAGCTGATATTTTATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGAt                                 <  1:3311392/91‑1 (MQ=255)
                                                            cGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGAtg                                <  1:1083111/91‑1 (MQ=255)
                                                            cGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGAtg                                <  2:1608962/91‑1 (MQ=255)
                                                              tttttACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGAtgtg                              <  1:3640510/91‑1 (MQ=255)
                                                                   aCCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCa                         <  1:3535745/91‑1 (MQ=255)
                                                                    ccAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCAt                        >  1:3453085/1‑91 (MQ=255)
                                                                     cAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATc                       <  2:1483206/91‑1 (MQ=255)
                                                                     cAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATc                       <  2:2890992/91‑1 (MQ=255)
                                                                      aaTCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATcc                      >  1:2719665/1‑91 (MQ=255)
                                                                      aaTCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATcc                      <  1:330500/91‑1 (MQ=255)
                                                                      aaTCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATcc                      <  2:1803147/91‑1 (MQ=255)
                                                                        tCAAACCAAGCTGATATTTAATATCAGCATAACTCGACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCAc                    <  2:2603331/91‑1 (MQ=255)
                                                                          aaaCCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACta                  <  2:2809006/91‑1 (MQ=255)
                                                                           aaCCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACtat                 <  1:2229867/91‑1 (MQ=255)
                                                                           aaCCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACtat                 <  1:2457782/91‑1 (MQ=255)
                                                                            aCCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACtata                <  2:1395304/91‑1 (MQ=255)
                                                                             ccAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATAc               <  1:1639185/91‑1 (MQ=255)
                                                                             ccAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATAc               <  2:1212884/91‑1 (MQ=255)
                                                                              cAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACa              >  2:481043/1‑91 (MQ=255)
                                                                              cAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACa              <  2:2868736/91‑1 (MQ=255)
                                                                               aaGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACaa             <  1:2142858/91‑1 (MQ=255)
                                                                               aaGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACaa             >  2:428893/1‑91 (MQ=255)
                                                                                aGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTg                                          >  2:610317/1‑61 (MQ=255)
                                                                                aGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAat            >  1:961768/1‑91 (MQ=255)
                                                                                aGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAat            >  1:2574304/1‑91 (MQ=255)
                                                                                  cTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATACTCTGGTAAGCGATGTGCTTCATCCACTATACAatca          >  2:3065351/1‑91 (MQ=255)
                                                                                  cTGATATTAAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAatca          >  1:3726779/1‑91 (MQ=255)
                                                                                   tGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAatcat         <  1:1511823/91‑1 (MQ=255)
                                                                                   tGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAatcat         <  2:1177691/91‑1 (MQ=255)
                                                                                   tGATATTAAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAatcat         >  1:618002/1‑91 (MQ=255)
                                                                                     atatTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTCAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAAtcatca       <  1:3868488/91‑1 (MQ=255)
                                                                                     atatTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAAtcatca       >  2:3248103/1‑91 (MQ=255)
                                                                                     atatTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAAtcatca       <  2:3702717/91‑1 (MQ=255)
                                                                                     atatTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAAtcatca       >  2:3823336/1‑91 (MQ=255)
                                                                                     atatTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAAtcatca       >  2:248583/1‑91 (MQ=255)
                                                                                     atatTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAAtcatca       >  1:2168021/1‑91 (MQ=255)
                                                                                     atatTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAAtcatca       <  2:2055096/91‑1 (MQ=255)
                                                                                      tatTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAATCATCaa      >  2:152310/1‑91 (MQ=255)
                                                                                      tatTAAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAATCATCaa      >  1:2331255/1‑91 (MQ=255)
                                                                                      tatTAAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAATCATCaa      >  1:1082081/1‑91 (MQ=255)
                                                                                        tttAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAATCATCAAAc    >  2:87006/1‑91 (MQ=255)
                                                                                          |                                                                                          
ATTTTTCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAATCATCAAACAA  >  USA300TCH1516_ALE/1526442‑1526622

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: