Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,934,780 A→T 5.9% N156K (AAT→AAA)  USA300HOU_RS09605 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,934,7800AT 5.9% 1206.3 / 42.6 370N156K (AAT→AAA) USA300HOU_RS09605hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (177/168);  new base T (12/11);  total (189/181)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 391 total aligned reads displayed.

ACTTAAGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATGAAATA  >  USA300TCH1516_ALE/1934690‑1934870
                                                                                          |                                                                                          
aCTTAAGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGaaa                                                                                            >  1:3018738/1‑91 (MQ=255)
 cTTAAGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAAt                                                                                           >  1:99740/1‑91 (MQ=255)
 cTTAAGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAAt                                                                                           >  1:111828/1‑91 (MQ=255)
 cTTAAGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAAt                                                                                           >  2:2927670/1‑91 (MQ=255)
 cTTAAGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAAt                                                                                           <  1:2265494/91‑1 (MQ=255)
    aaGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCa                                                                                        >  1:925427/1‑91 (MQ=255)
    aaGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCa                                                                                        >  1:1813347/1‑91 (MQ=255)
     aGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCat                                                                                       >  2:3341203/1‑90 (MQ=255)
        aTTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCaaata                                                                                    <  2:1295992/91‑1 (MQ=255)
          tttttCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCgaataaa                                                                                  <  1:3076559/91‑1 (MQ=255)
           ttttcttttcTTCTAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAATTTAGAGAATATTTTTTTGAAATTCaaataaat                                                                                 <  2:2515137/91‑1 (MQ=255)
           ttttcttttcTTCAAACTCTTCTTCCAGTTCCCAATAATTATTGACTGATTGTTCTAACTTAGAGAATATTTTTTTGAAATTCaaataaat                                                                                 >  2:3653608/1‑91 (MQ=255)
           ttttcttttcTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGCATATTTTTTTGAAATTCaaataaat                                                                                 >  1:3688342/1‑91 (MQ=255)
           ttttcttttcTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCaaataaat                                                                                 >  1:1731311/1‑91 (MQ=255)
             ttcttttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAatat                                                                               >  1:1284673/1‑91 (MQ=255)
               cttttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTc                                                                             >  1:2457461/1‑91 (MQ=255)
                 tttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGTGAATATTTTTTTGAAATTCAAATATATATTCtt                                                                           <  2:142481/91‑1 (MQ=255)
                 tttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCtt                                                                           >  2:1386623/1‑91 (MQ=255)
                   tcttcAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTg                                                                         <  1:1076240/91‑1 (MQ=255)
                   tcttcAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTATTTGAATTTCAAAAAAATATTCtcta                                                                         <  1:3360392/91‑4 (MQ=25)
                    cttcAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGc                                                                        >  2:2799964/1‑91 (MQ=255)
                    cttcAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGc                                                                        >  2:3230556/1‑91 (MQ=255)
                    cttcAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATAGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGc                                                                        >  1:1570281/1‑91 (MQ=255)
                     ttcAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCg                                                                       <  1:1638816/91‑1 (MQ=255)
                        aaaCTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTAt                                                                    <  2:3456136/91‑1 (MQ=255)
                         aaCTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTAtt                                                                   >  2:1223481/1‑91 (MQ=255)
                         aaCTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTAtt                                                                   >  2:283589/1‑91 (MQ=255)
                         aaCTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTAtt                                                                   >  2:559719/1‑91 (MQ=255)
                         aaCTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTAtt                                                                   <  2:3414079/91‑1 (MQ=255)
                          aCTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTTAAATTCAAATAAATATTCTTTGCGTAttt                                                                  >  2:1668265/1‑91 (MQ=255)
                           ctctTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTAtttt                                                                 >  1:2057196/1‑91 (MQ=255)
                           ctctTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTAtttt                                                                 >  1:1392620/1‑91 (MQ=255)
                           ctctTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGAGTAtttt                                                                 >  1:1959717/1‑91 (MQ=255)
                             cttcttCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGa                                                               <  1:2942570/91‑1 (MQ=255)
                                cttcCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGAttt                                                            <  1:2203814/91‑1 (MQ=255)
                                cttcCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGAttt                                                            >  2:340148/1‑91 (MQ=255)
                                 ttcCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTc                                                           >  2:1733950/1‑91 (MQ=255)
                                  tcCAGTTCCAAATAATTATTGACTGATTGTTCTCAATAAGAGAATACTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCa                                                          >  1:3151201/1‑91 (MQ=255)
                                     aGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCt                                                       <  2:2871241/91‑1 (MQ=255)
                                     aGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCt                                                       <  2:435939/91‑1 (MQ=255)
                                     aGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCt                                                       <  1:3869848/91‑1 (MQ=255)
                                      gTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCtt                                                      >  2:94293/1‑91 (MQ=255)
                                      gTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCtt                                                      <  1:2869789/91‑1 (MQ=255)
                                      gTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCtt                                                      >  1:3764781/1‑91 (MQ=255)
                                       ttCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTa                                                     <  2:3140157/91‑1 (MQ=255)
                                          cAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCt                                                  <  2:679461/91‑1 (MQ=255)
                                          cAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCt                                                  <  2:1957890/91‑1 (MQ=255)
                                          cAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCt                                                  >  2:2219340/1‑91 (MQ=255)
                                          cAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCt                                                  >  1:3373441/1‑91 (MQ=255)
                                          cAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCt                                                  >  1:2768396/1‑91 (MQ=255)
                                            aataatTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTaa                                                >  1:3261244/1‑91 (MQ=255)
                                            aataatTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTaa                                                >  1:2592645/1‑91 (MQ=255)
                                             ataatTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTaaa                                               >  2:504367/1‑91 (MQ=255)
                                                attattGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAata                                            >  2:2029578/1‑91 (MQ=255)
                                                attattGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAata                                            >  1:2509502/1‑91 (MQ=255)
                                                 ttattGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAatat                                           <  1:1474013/91‑1 (MQ=255)
                                                  tattGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATc                                          >  1:3804571/1‑91 (MQ=255)
                                                  tattGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATc                                          <  1:1779374/91‑1 (MQ=255)
                                                  tattGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATc                                          <  2:1881251/91‑1 (MQ=255)
                                                    ttGACTGATTGTTCTAAATTAGAGAATATTTTTTTTATTTTTTTTTAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCAt                                        >  2:2636536/1‑91 (MQ=255)
                                                    ttGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCAt                                        <  2:3327467/91‑1 (MQ=255)
                                                    ttGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCAt                                        >  2:3622908/1‑91 (MQ=255)
                                                    ttGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCAt                                        >  1:1023225/1‑91 (MQ=255)
                                                     tGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCAtt                                       <  1:2093629/91‑1 (MQ=255)
                                                     tGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCAtt                                       >  1:2352116/1‑91 (MQ=255)
                                                        cTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTc                                    >  1:595149/1‑91 (MQ=255)
                                                        cTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTc                                    <  1:617669/91‑1 (MQ=255)
                                                             tGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTTGatat                               >  2:3857046/1‑91 (MQ=255)
                                                             tGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGatat                               <  1:971110/91‑1 (MQ=255)
                                                               ttCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATaa                             >  2:878239/1‑91 (MQ=255)
                                                                tCTAAATTAGAGAATATTTTTTTGATATTCAAATAAATATTCTTTGCCTATTTTGATTTCATCTTATCTAAAATATTATTTTCGATATAAc                            >  2:259873/1‑91 (MQ=255)
                                                                tCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAAc                            >  1:653906/1‑91 (MQ=255)
                                                                 cTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAAcc                           >  2:2092989/1‑91 (MQ=255)
                                                                  tAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCt                          >  2:629298/1‑91 (MQ=255)
                                                                   aaaTTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCtt                         >  1:790831/1‑91 (MQ=255)
                                                                     aTTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCtttt                       <  2:77068/91‑1 (MQ=255)
                                                                       tAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTaa                     <  2:2276491/91‑1 (MQ=255)
                                                                        agagAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAACATATCATTTTCGATATAACCTTTTAAt                    >  2:2395437/1‑91 (MQ=255)
                                                                        agagAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATATCCTTTTAAt                    >  2:2943852/1‑91 (MQ=255)
                                                                        agagAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAAt                    <  1:479412/91‑1 (MQ=255)
                                                                        agagAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAAt                    >  1:3507656/1‑91 (MQ=255)
                                                                        agagAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATAATTTTCGATATAACCTTTTAAt                    >  2:3182324/1‑91 (MQ=255)
                                                                         gagaATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAAtt                   <  1:241148/91‑1 (MQ=255)
                                                                           gaATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCt                 >  1:556547/1‑91 (MQ=255)
                                                                           gaATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCt                 >  2:1333196/1‑91 (MQ=255)
                                                                            aaTATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCtt                <  2:3757153/91‑1 (MQ=255)
                                                                            aaTATTTATTTGAATTTCAAAAAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCtt                >  2:1208149/1‑91 (MQ=38)
                                                                             atatTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCttt               <  2:2408643/91‑1 (MQ=255)
                                                                              tatTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCtttt              >  1:2018003/1‑91 (MQ=255)
                                                                              tatTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCtttt              >  2:3795170/1‑91 (MQ=255)
                                                                               atTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTa             >  2:2148202/1‑91 (MQ=255)
                                                                                  tttttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGat          <  1:1704293/91‑1 (MQ=255)
                                                                                  tttttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGat          <  1:1221510/91‑1 (MQ=255)
                                                                                   ttttGAAATTCAAATAAATATTCTTTGCGTGTTTTGATTTCATCTTATCTAAAAGATCATTTTCGATATAACCTTTTAAGTGTTATAGata         <  2:1162352/91‑1 (MQ=255)
                                                                                    tttGAATTTCAAAAAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGatat        >  1:2377432/1‑91 (MQ=255)
                                                                                    tttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGatat        <  2:1595491/91‑1 (MQ=255)
                                                                                      tGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATGa      <  2:1038832/91‑1 (MQ=255)
                                                                                         aaTTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATGAAat   >  1:1343522/1‑91 (MQ=255)
                                                                                         aaTTCAAATAAATATTCTTTGCGTATTTTGATCTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATGAAat   >  2:3887448/1‑91 (MQ=255)
                                                                                          aTTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATGAAata  >  2:719088/1‑91 (MQ=255)
                                                                                          |                                                                                          
ACTTAAGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATGAAATA  >  USA300TCH1516_ALE/1934690‑1934870

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