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% 801.6 / 36.2 235K471* (AAA→TAA) dinG_1putative ATP‑dependent helicase DinG
Reads supporting (aligned to +/- strand):  ref base T (100/118);  new base A (6/11);  total (106/129)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 250 total aligned reads displayed.

ATTTTTCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAATCATCAAACAA  >  USA300TCH1516_ALE/1526442‑1526622
                                                                                          |                                                                                          
aTTTTTCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATAttt                                                                                            >  1:693519/1‑91 (MQ=255)
  ttttGCAAAATTCGCTCCCTTTCCAGTTGAACAATTGCCTTGAACAAATTTTCCTATTCGCTTTTACCAATCAAACCAATCTAATCTCTaa                                                                                          >  1:1731796/1‑91 (MQ=255)
  ttttCCAAAAGTCGCTGCTTTTCCAGTTGATCAATTGCCGTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTaaa                                                                                          <  2:283855/91‑1 (MQ=255)
   tttCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTAAat                                                                                         <  2:2519930/91‑1 (MQ=255)
    ttCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAata                                                                                        >  1:968247/1‑91 (MQ=255)
    ttCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTGTTTACCAATCAAACCAAGCTGATATTTAata                                                                                        <  2:1078291/91‑1 (MQ=255)
       cAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCa                                                                                     >  1:2720080/1‑91 (MQ=255)
       cAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATACTTAATATCa                                                                                     >  1:3779420/1‑91 (MQ=255)
        aaaaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAg                                                                                    >  1:3612099/1‑91 (MQ=255)
        aaaaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAg                                                                                    >  2:2368636/1‑91 (MQ=255)
        aaaaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAg                                                                                    >  2:1501942/1‑91 (MQ=255)
         aaaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGc                                                                                   >  1:3529798/1‑91 (MQ=255)
         aaaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGc                                                                                   >  1:2064269/1‑91 (MQ=255)
         aaaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGc                                                                                   >  1:2057267/1‑91 (MQ=255)
         aaaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGc                                                                                   <  1:1199047/91‑1 (MQ=255)
          aaTTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCa                                                                                  <  2:311633/91‑1 (MQ=255)
             tCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTAAATATCAGcttgg                                                                               <  2:3724077/91‑5 (MQ=255)
              cGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaac                                                                              >  2:115369/1‑91 (MQ=255)
               gctgctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaact                                                                             >  1:1816125/1‑91 (MQ=255)
                ctgctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactc                                                                            <  1:749129/91‑1 (MQ=255)
                ctgctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactc                                                                            <  1:221291/91‑1 (MQ=255)
                ctgctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactc                                                                            >  2:3465610/1‑91 (MQ=255)
                 tgctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactca                                                                           >  2:2330093/1‑91 (MQ=255)
                 tgctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactca                                                                           >  1:555554/1‑91 (MQ=255)
                  gctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactcaa                                                                          <  1:1501439/91‑1 (MQ=255)
                  gctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactcaa                                                                          <  1:1470039/91‑1 (MQ=255)
                  gctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactcaa                                                                          <  1:2807180/91‑1 (MQ=255)
                  gctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCACaactcaa                                                                          >  1:657082/1‑91 (MQ=255)
                   ctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactcaac                                                                         <  1:2334008/91‑1 (MQ=255)
                   ctTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactcaac                                                                         <  2:33103/91‑1 (MQ=255)
                    ttttCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactcaact                                                                        >  1:1024748/1‑91 (MQ=255)
                     tttCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATaactcaactc                                                                       <  1:836373/91‑1 (MQ=255)
                         cAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCAttt                                                                   >  2:3404542/1‑91 (MQ=255)
                         cAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCAttt                                                                   >  1:421750/1‑91 (MQ=255)
                          aGTTGATCAATTGCCTTTAACAATTTTTCATTTGCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTg                                                                  <  1:3245174/91‑1 (MQ=255)
                            ttGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTc                                                                >  2:1368347/1‑91 (MQ=255)
                            ttGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTc                                                                >  1:2140191/1‑91 (MQ=255)
                             tGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTcc                                                               >  2:1519418/1‑90 (MQ=255)
                             tGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCa                                                               >  2:3875218/1‑91 (MQ=255)
                              gATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCAc                                                              <  1:1461247/91‑1 (MQ=255)
                              gATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCAc                                                              <  2:1359027/91‑1 (MQ=255)
                              gATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCAc                                                              <  2:3227793/91‑1 (MQ=255)
                               aTCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACt                                                             <  1:3839266/91‑1 (MQ=255)
                                tCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATTAGCATAACTCAACTCAATTGTTACtt                                                            >  2:1694107/1‑91 (MQ=255)
                                tCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACtt                                                            <  2:2987625/91‑1 (MQ=255)
                                    ttGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGAtt                                                        <  2:567007/91‑1 (MQ=255)
                                    ttGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGAtt                                                        >  2:1638675/1‑91 (MQ=255)
                                    ttGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGAtt                                                        >  1:3799643/1‑91 (MQ=255)
                                     tGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGAttt                                                       <  1:2535491/91‑1 (MQ=255)
                                       ccTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTaa                                                     <  2:1971512/91‑1 (MQ=255)
                                         tttAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAgg                                                   <  1:1240070/91‑1 (MQ=255)
                                         tttAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAgg                                                   <  2:422057/91‑1 (MQ=255)
                                         tttAACAAGTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAgg                                                   <  2:792520/91‑1 (MQ=255)
                                           tAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCa                                                 <  2:2573416/91‑1 (MQ=255)
                                           tAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCa                                                 <  1:1635544/91‑1 (MQ=255)
                                              cAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAg                                              <  1:3067181/91‑1 (MQ=255)
                                                 tttttCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCt                                           <  1:3296730/91‑1 (MQ=255)
                                                   tttCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTgg                                         >  2:2925390/1‑91 (MQ=255)
                                                   tttCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTgg                                         <  1:3718774/91‑1 (MQ=255)
                                                     tCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTa                                       <  1:100417/91‑1 (MQ=255)
                                                       aTTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAg                                     <  2:855176/91‑1 (MQ=255)
                                                       aTTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAg                                     <  2:2180889/91‑1 (MQ=255)
                                                        ttttCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGc                                    >  2:2594581/1‑91 (MQ=255)
                                                         tttcttttctACAACTCAATCCAAGCTGATATTTAATATCGGCATAACTCAACTCAATTGTGACTAGATTTAAGGCATAGTCTGGTAAGCg                                   <  2:3442188/91‑1 (MQ=255)
                                                         tttCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCg                                   <  2:3554349/91‑1 (MQ=255)
                                                         tttCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCg                                   <  1:1316842/91‑1 (MQ=255)
                                                          ttCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGa                                  >  1:1761622/1‑91 (MQ=255)
                                                            cGTTTTTACCAATCAAACCAAGCTGCTATCTAACATCAGCATAACCAAACTCATTTGTCCCTTCATTCAAGGCATCGTCTGGTCAGCGAtg                                >  2:3715786/1‑91 (MQ=255)
                                                            cGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGAtg                                <  2:2876833/91‑1 (MQ=255)
                                                            cGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGAtg                                >  2:971430/1‑91 (MQ=255)
                                                            cGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGAtg                                >  2:1875418/1‑91 (MQ=255)
                                                              tttttACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGAtgtg                              <  1:497952/91‑1 (MQ=255)
                                                              tttttACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGAtgtg                              <  1:633637/91‑1 (MQ=255)
                                                               ttttACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGc                             <  2:678768/91‑1 (MQ=255)
                                                                   aCCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCa                         <  1:1570172/91‑1 (MQ=255)
                                                                    ccAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCAt                        >  1:2863222/1‑91 (MQ=255)
                                                                       aTCAAACCAAGCTGATATTTAATATCAGCATAACTCATCTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCa                     <  1:1311598/91‑1 (MQ=255)
                                                                       aTCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCa                     >  2:555263/1‑91 (MQ=255)
                                                                       aTCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCa                     >  1:364981/1‑91 (MQ=255)
                                                                        tCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCAc                    <  1:1218699/91‑1 (MQ=255)
                                                                         cAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACt                   >  1:713482/1‑91 (MQ=255)
                                                                          aaaCCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACta                  <  1:1868450/91‑1 (MQ=255)
                                                                            aCCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACtata                >  1:3592354/1‑91 (MQ=255)
                                                                              cAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACa              <  1:1090750/91‑1 (MQ=255)
                                                                              cAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACa              <  2:2278019/91‑1 (MQ=255)
                                                                                aGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAat            >  2:1805151/1‑91 (MQ=255)
                                                                                aGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAat            >  1:3744352/1‑91 (MQ=255)
                                                                                  cTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAatca          >  1:1983783/1‑91 (MQ=255)
                                                                                  cTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAatca          <  1:2715476/91‑1 (MQ=255)
                                                                                    gATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAatcatc        >  1:1331127/1‑91 (MQ=255)
                                                                                    gATATTAAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAatcatc        >  2:3032076/1‑91 (MQ=255)
                                                                                     atatTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAAtcatca       <  1:1991978/91‑1 (MQ=255)
                                                                                      tatTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAATCATCaa      >  1:3476960/1‑91 (MQ=255)
                                                                                      tatTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAATCATCaa      <  1:3466597/91‑1 (MQ=255)
                                                                                      tatTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAATCATCaa      <  2:2216634/91‑1 (MQ=255)
                                                                                      tatTAAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAATCATCaa      >  1:934190/1‑91 (MQ=255)
                                                                                        tttAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATGCACTATACAATCATCAAAc    >  2:345768/1‑91 (MQ=255)
                                                                                        tttAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAATCATCAAAc    >  2:1961412/1‑91 (MQ=255)
                                                                                        tttAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAATCATCAAAc    >  1:3095503/1‑91 (MQ=255)
                                                                                          tAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAATCATCAAACaa  >  2:769621/1‑91 (MQ=255)
                                                                                          |                                                                                          
ATTTTTCCAAAATTCGCTGCTTTTCCAGTTGATCAATTGCCTTTAACAATTTTTCATTTTCGTTTTTACCAATCAAACCAAGCTGATATTTAATATCAGCATAACTCAACTCATTTGTCACTTGATTTAAGGCATAGTCTGGTAAGCGATGTGCTTCATCCACTATACAATCATCAAACAA  >  USA300TCH1516_ALE/1526442‑1526622

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