Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 70,101 A→T 6.1% I228N (ATT→AAT)  argF_1 ← Ornithine carbamoyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE70,1010AT 6.1% 661.6 / 23.2 195I228N (ATT→AAT) argF_1Ornithine carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base A (89/94);  new base T (5/7);  total (94/101)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 206 total aligned reads displayed.

TCTTTAGTTACACGGTATGGTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTATATCATTACACCGATTTA  >  USA300TCH1516_ALE/70011‑70191
                                                                                           |                                                                                          
 cTTTAGTTACACGGTATGGTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAAt                                                                                           >  2:1494888/1‑91 (MQ=255)
  tttAGTTACACGGTATGGTTCTAATAATTTGATACGTTTTACCCAAACTTCATCAGGTTCACCCATCGATACCCAAACATCTGTGTAAAtt                                                                                          >  2:3600295/1‑91 (MQ=255)
    tAGTTACACGGTATGGTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTAc                                                                                        >  2:586371/1‑91 (MQ=255)
    tAGTTACACGGTATGGTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTAc                                                                                        <  2:280181/91‑1 (MQ=255)
       ttACACGGTATGGTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGATGTAATTTACac                                                                                      <  1:2565850/91‑2 (MQ=38)
       ttACACGGTATGGTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATc                                                                                     <  2:1297397/91‑1 (MQ=255)
        tACACGGTATGGTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGATGTAATTTACaca                                                                                     <  2:3115218/91‑3 (MQ=38)
        tACACGGTATGGTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGATGTAATTTACaca                                                                                     <  1:1708528/91‑3 (MQ=38)
         acacGGTATGGTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAg                                                                                   >  1:364895/1‑91 (MQ=255)
         acacGGTATGGTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAg                                                                                   <  1:3859794/91‑1 (MQ=255)
             ggTATGGTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATcc                                                                               >  1:1432270/1‑91 (MQ=255)
               tATGGTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCtt                                                                             >  2:633336/1‑91 (MQ=255)
                aTGGTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCttt                                                                            <  1:476207/91‑1 (MQ=255)
                aTGGTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCttt                                                                            <  2:3256028/91‑1 (MQ=255)
                aTGGTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCttt                                                                            >  1:1375030/1‑91 (MQ=255)
                  ggTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGa                                                                          <  2:601415/91‑1 (MQ=255)
                  ggTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGa                                                                          <  2:2004311/91‑1 (MQ=255)
                   gTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGac                                                                         >  1:139064/1‑91 (MQ=255)
                   gTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGac                                                                         <  1:236163/91‑1 (MQ=255)
                        aataatTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCtt                                                                    >  1:2112002/1‑91 (MQ=255)
                        aataatTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCtt                                                                    <  2:2084113/91‑1 (MQ=255)
                         ataatTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTtc                                                                   <  2:283295/91‑1 (MQ=255)
                            atTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTtcatc                                                                <  2:2624686/91‑1 (MQ=255)
                            atTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTtcatc                                                                >  2:489367/1‑91 (MQ=255)
                            atTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTtcatc                                                                >  2:430667/1‑91 (MQ=255)
                             tttGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTtcatca                                                               >  1:3174456/1‑91 (MQ=255)
                               tGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAat                                                             >  1:2848777/1‑91 (MQ=255)
                                    cGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATatcat                                                        <  1:160178/91‑1 (MQ=255)
                                     gTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATatcatc                                                       <  1:3737509/91‑1 (MQ=255)
                                     gTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATatcatc                                                       <  1:1079242/91‑1 (MQ=255)
                                       ttttCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTTACACCTTCATCAATATCATCAg                                                     >  1:2010397/1‑91 (MQ=255)
                                        tttCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGt                                                    <  2:1011545/91‑1 (MQ=255)
                                        tttCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGt                                                    >  2:293197/1‑91 (MQ=255)
                                        tttCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGt                                                    <  2:364981/91‑1 (MQ=255)
                                          tCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTta                                                  >  2:2697489/1‑91 (MQ=255)
                                          tCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTta                                                  >  1:186595/1‑91 (MQ=255)
                                          tCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTta                                                  <  2:2812434/91‑1 (MQ=255)
                                           cccAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTtat                                                 <  2:2019198/91‑1 (MQ=255)
                                            ccAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTtata                                                <  1:911668/91‑1 (MQ=255)
                                             cAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATaa                                               <  2:3424127/91‑1 (MQ=255)
                                             cAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATaa                                               >  1:1723679/1‑91 (MQ=255)
                                              aaaCTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAg                                              >  1:3759170/1‑91 (MQ=255)
                                              aaaCTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAg                                              >  2:3025857/1‑91 (MQ=255)
                                               aaCTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAgg                                             <  1:3583515/91‑1 (MQ=255)
                                               aaCTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAgg                                             <  2:1692489/91‑1 (MQ=255)
                                                     atcaGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTca                                       >  1:1051217/1‑91 (MQ=255)
                                                     atcaGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTca                                       >  2:1881778/1‑91 (MQ=255)
                                                     atcaGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTca                                       >  1:1238253/1‑91 (MQ=255)
                                                     atcaGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTca                                       >  2:2498929/1‑91 (MQ=255)
                                                     atcaGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTca                                       >  2:1596567/1‑91 (MQ=255)
                                                     atcaGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTca                                       >  1:1615716/1‑91 (MQ=255)
                                                     atcaGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTca                                       >  1:2172240/1‑91 (MQ=255)
                                                      tcaGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTcac                                      >  1:2152910/1‑91 (MQ=255)
                                                      tcaGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTcac                                      >  2:2999797/1‑91 (MQ=255)
                                                       caGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTcacc                                     <  2:283818/91‑1 (MQ=255)
                                                       caGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTcacc                                     <  1:2446401/91‑1 (MQ=255)
                                                         ggTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTcaccac                                   <  1:536427/91‑1 (MQ=255)
                                                         ggTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTcaccac                                   <  2:865392/91‑1 (MQ=255)
                                                          gTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCaccacc                                  <  2:1379620/91‑1 (MQ=255)
                                                             cACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGtt                               <  2:2292345/91‑1 (MQ=255)
                                                              aCCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGttt                              >  2:2244124/1‑91 (MQ=255)
                                                              aCCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGttt                              <  2:920736/91‑1 (MQ=255)
                                                              aCCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGttt                              <  1:2640293/91‑1 (MQ=255)
                                                              aCCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGttt                              >  1:3862281/1‑91 (MQ=255)
                                                               cccATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGtttt                             <  1:1015482/91‑1 (MQ=255)
                                                                ccATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGttttt                            <  1:3362034/91‑1 (MQ=255)
                                                                  aTAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGt                          <  2:2155094/91‑1 (MQ=255)
                                                                  aTAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGt                          <  2:2152498/91‑1 (MQ=255)
                                                                  aTAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGt                          <  2:2870854/91‑1 (MQ=255)
                                                                   tAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTc                         >  2:6786/1‑91 (MQ=255)
                                                                   tAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTc                         <  1:1216566/91‑1 (MQ=255)
                                                                   tAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTc                         >  2:1127387/1‑91 (MQ=255)
                                                                   tAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTc                         >  1:2336141/1‑91 (MQ=255)
                                                                      aTACCGAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGc                      <  2:3575357/91‑1 (MQ=255)
                                                                       tACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCt                     >  2:2270886/1‑91 (MQ=255)
                                                                       tACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCt                     <  1:1494888/91‑1 (MQ=255)
                                                                       tACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCt                     <  1:2810649/91‑1 (MQ=255)
                                                                         cccAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCtat                   <  2:1573457/91‑1 (MQ=255)
                                                                         cccAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCtat                   >  2:3357807/1‑91 (MQ=255)
                                                                         cccAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCtat                   <  1:3679980/91‑1 (MQ=255)
                                                                           cAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTatat                 <  2:1783832/91‑1 (MQ=255)
                                                                           cAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTatat                 <  1:296427/91‑1 (MQ=255)
                                                                             aaCATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTATATCa               <  2:2202913/91‑1 (MQ=255)
                                                                              aCATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTATATCAt              <  1:1142582/91‑1 (MQ=255)
                                                                               cATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTATATCAtt             <  2:1160791/91‑1 (MQ=255)
                                                                                 tCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTGTATCATTac           >  1:3483603/1‑91 (MQ=255)
                                                                                  cTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTATATCATTaca          <  2:1451434/91‑1 (MQ=255)
                                                                                   tgtgTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTATATCATTacac         >  1:3817511/1‑91 (MQ=255)
                                                                                   tgtgTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTATATCATTacac         <  1:2498929/91‑1 (MQ=255)
                                                                                    gatgtAATTTACACAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTATATCATTACAcc        >  2:1243046/3‑91 (MQ=255)
                                                                                    gatgtAATTTACACAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTATATCATTACAcc        >  2:2901590/3‑91 (MQ=255)
                                                                                    gatgtAATTTACACAGATCCTTTCACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTATATCATTACAcc        >  2:3507759/3‑91 (MQ=255)
                                                                                    gtgtAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTATATCATTACAcc        >  2:776657/1‑91 (MQ=255)
                                                                                    gtgtAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTATATCATTACAcc        <  1:3436406/91‑1 (MQ=255)
                                                                                     atgtAATTTACACAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTATATCATTACACCg       >  1:3693689/2‑91 (MQ=255)
                                                                                      tgtAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTATATCATTACACCg       >  1:1574100/1‑91 (MQ=255)
                                                                                      tgtAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTATATCATTACACCg       <  1:2662787/91‑1 (MQ=255)
                                                                                         aaaTTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTATATCATTACACCGAtt    <  1:1231302/91‑1 (MQ=255)
                                                                                         aaaTTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTATATCATTACACCGAtt    <  2:3381843/91‑1 (MQ=255)
                                                                                           aTTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTATATCATTACACCGATTta  <  1:619657/91‑1 (MQ=255)
                                                                                           |                                                                                          
TCTTTAGTTACACGGTATGGTTCTAATAATTTGATACGTTTTTCCCAAACTTCATCAGGTTCACCCATAGATACCCAAACATCTGTGTAAATTACATCAGATCCTTTGACACCTTCATCAATATCATCAGTTATAAGGATTTCACCACCGTTTTTGTCAGCTATATCATTACACCGATTTA  >  USA300TCH1516_ALE/70011‑70191

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