Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,870,863 C→T 7.0% intergenic (‑350/+200) srkA ← / ← dat Stress response kinase A/D‑alanine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,870,8630CT 7.0% 768.4 / 30.6 242intergenic (‑350/+200)srkA/datStress response kinase A/D‑alanine aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (117/107);  new base T (7/10);  total (124/118)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 265 total aligned reads displayed.

ACTCTAGTTTTCAAAATAATTGCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGCACTTTTTACTAAT  >  USA300TCH1516_ALE/1870773‑1870953
                                                                                          |                                                                                          
aCTCTAGTTTTCAAAATAATTGCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGc                                                                                            <  2:219953/91‑1 (MQ=255)
  tctAGTTTTCAAAATAATTGCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGcac                                                                                          <  2:3868154/91‑1 (MQ=255)
  tctAGTTTTCAAAATAATTGCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGcac                                                                                          >  1:2958200/1‑91 (MQ=255)
   ctAGTTTTCAAAATAATTGCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGcaca                                                                                         <  1:94196/91‑1 (MQ=255)
   ctAGTTTTCAAAATAATTGCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGcaca                                                                                         <  2:2122618/91‑1 (MQ=255)
   ctAGTTTTCAAAATAATTGCTATCCCATTTTCATGGTTAGAATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGcaca                                                                                         >  1:1544368/1‑91 (MQ=255)
   ctAGTTTTCAAAAGAATTGCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGcaca                                                                                         <  2:1729504/91‑1 (MQ=255)
     aGTTTTCAAAATAATTGCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGTGcatc                                                                                       <  1:2591447/91‑1 (MQ=255)
       ttttCAAAATAATTGCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGTGcatcat                                                                                     <  2:2559065/91‑1 (MQ=255)
       ttttCAAAATAATTGCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCAcatcat                                                                                     >  1:532538/1‑91 (MQ=255)
        tttCAAAATAATTGCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGTGCATCata                                                                                    <  1:2420086/91‑1 (MQ=255)
        tttCAAAATAATTGCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCata                                                                                    <  1:1145537/91‑1 (MQ=255)
         ttCAAAATAATTGCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCatat                                                                                   >  2:1566307/1‑91 (MQ=255)
            aaaaTAATTGCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGTGCATCATAtttt                                                                                <  1:2123825/91‑2 (MQ=255)
               ataatTGCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTaa                                                                             >  2:3406959/1‑91 (MQ=255)
                  atTGCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTc                                                                          >  2:2854347/1‑91 (MQ=255)
                   ttGCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCa                                                                         <  1:976388/91‑1 (MQ=255)
                   ttGCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCa                                                                         >  1:390035/1‑91 (MQ=255)
                     gCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATa                                                                       <  2:1671933/91‑1 (MQ=255)
                     gCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATa                                                                       >  2:3447356/1‑91 (MQ=255)
                     gCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATa                                                                       <  1:736708/91‑1 (MQ=255)
                     gCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATa                                                                       >  1:904780/1‑91 (MQ=255)
                      cTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAg                                                                      <  1:1448855/91‑1 (MQ=255)
                            cATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAAtt                                                                <  1:555727/91‑1 (MQ=255)
                            cATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAAtt                                                                <  1:3091428/91‑1 (MQ=255)
                             aTTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAAttt                                                               >  2:2102132/1‑91 (MQ=255)
                             aTTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAAttt                                                               <  1:576536/91‑1 (MQ=255)
                             aTTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAAttt                                                               <  2:2765690/91‑1 (MQ=255)
                              ttttCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAAtttt                                                              <  1:3199011/91‑1 (MQ=255)
                              ttttCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAAtttt                                                              <  1:727077/91‑1 (MQ=255)
                              ttttCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAAtttt                                                              <  1:1592170/91‑1 (MQ=255)
                              ttttCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAAtttt                                                              <  2:2440835/91‑1 (MQ=255)
                               tttCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTta                                                             >  2:1232977/1‑91 (MQ=255)
                                ttCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTtat                                                            <  1:1440742/91‑1 (MQ=255)
                                ttCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTtat                                                            <  2:2486280/91‑1 (MQ=255)
                                ttCATGGTTAGCATAAATTTATGAACTGTAACATATACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTtat                                                            <  2:823405/91‑1 (MQ=255)
                                 tCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTtata                                                           <  1:1245005/91‑1 (MQ=255)
                                 tCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTtata                                                           >  2:1462158/1‑91 (MQ=255)
                                 tCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTtata                                                           >  2:3032294/1‑91 (MQ=255)
                                  cATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATaa                                                          <  2:3671218/91‑1 (MQ=255)
                                    tGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAAtt                                                        >  1:3504177/1‑91 (MQ=255)
                                    tGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAAtt                                                        >  1:1960985/1‑91 (MQ=255)
                                        tAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTa                                                    <  2:2284221/91‑1 (MQ=255)
                                        tAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTa                                                    <  2:3079890/91‑1 (MQ=255)
                                          gCATAAATTTATGAACTGTAACATTTTCGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATc                                                  >  2:3217987/1‑91 (MQ=255)
                                          gCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGTTGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATc                                                  <  1:3495485/91‑1 (MQ=255)
                                          gCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATc                                                  <  2:2718652/91‑1 (MQ=255)
                                          gCAGAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATc                                                  <  1:267476/91‑1 (MQ=255)
                                           cATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCa                                                 <  1:2485817/91‑1 (MQ=255)
                                             tAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCAtt                                               <  1:2887675/91‑1 (MQ=255)
                                              aaaTTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATta                                              <  2:2019612/91‑1 (MQ=255)
                                               aaTTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATtat                                             >  1:3387120/1‑91 (MQ=255)
                                                aTTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATtata                                            >  1:1977237/1‑91 (MQ=255)
                                                   tATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATAttt                                         <  2:3521328/91‑1 (MQ=255)
                                                    aTGAACTGTAACATTTACGTACTTAGTAAAATATGATTCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTc                                        >  2:2525131/1‑91 (MQ=255)
                                                    aTGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTc                                        <  2:2549934/91‑1 (MQ=255)
                                                     tGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCa                                       <  2:580905/91‑1 (MQ=255)
                                                     tGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCa                                       <  2:3775020/91‑1 (MQ=255)
                                                      gAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCaa                                      <  1:1313690/91‑1 (MQ=255)
                                                       aaCTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAAc                                     >  1:2045831/1‑91 (MQ=255)
                                                       aaCTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAAc                                     >  2:1162969/1‑91 (MQ=255)
                                                        aCTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACt                                    <  1:3064256/91‑1 (MQ=255)
                                                           gTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGaa                                 >  2:3341079/1‑91 (MQ=255)
                                                           gTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGaa                                 >  2:2042259/1‑91 (MQ=255)
                                                            tAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGaaa                                <  1:1813215/91‑1 (MQ=255)
                                                              aCATTTACGTACTTAGTAAAATATGATGCACATCATACTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAAt                              >  1:3258049/1‑91 (MQ=255)
                                                                 tttACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATgag                           <  1:3541117/91‑1 (MQ=255)
                                                                 tttACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATgag                           >  2:643309/1‑91 (MQ=255)
                                                                 tttACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACCGAAAATgag                           >  1:3718508/1‑91 (MQ=255)
                                                                  ttACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATgaga                          >  1:3807365/1‑91 (MQ=255)
                                                                   tACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGaa                         <  1:2778566/91‑1 (MQ=255)
                                                                   tACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGaa                         <  1:193851/91‑1 (MQ=255)
                                                                     cGTACTTAGTAAAATATGATGCACATTATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAAc                       >  2:2279011/1‑91 (MQ=255)
                                                                     cGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATTGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAAc                       >  2:2326267/1‑91 (MQ=255)
                                                                         cTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACaaaa                   >  2:755473/1‑91 (MQ=255)
                                                                          ttAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAAt                  <  1:2504238/91‑1 (MQ=255)
                                                                          ttAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAAt                  <  2:2894088/91‑1 (MQ=255)
                                                                          ttAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAAt                  <  1:3346206/91‑1 (MQ=255)
                                                                          ctAGTAAAATATGATGCACATCATATTTGTAACTCAAAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAAt                  <  2:3845656/90‑1 (MQ=255)
                                                                             gTAAAATATGATGCACATCATATTTGTAACTCATAGTAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGc               <  2:3504177/91‑1 (MQ=255)
                                                                             gTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGc               >  1:660292/1‑91 (MQ=255)
                                                                             gTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGc               >  1:3773479/1‑91 (MQ=255)
                                                                             gTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGc               >  2:3584840/1‑91 (MQ=255)
                                                                             gTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGc               <  1:1534637/91‑1 (MQ=255)
                                                                              tAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGCa              <  2:3630181/91‑1 (MQ=255)
                                                                              tAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGCa              >  1:3630095/1‑91 (MQ=255)
                                                                               caaaTATGATGTGCATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGCAc             >  1:367667/2‑91 (MQ=255)
                                                                               aaaaTATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGCAc             >  1:238817/1‑91 (MQ=255)
                                                                               aaaaTATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGCAc             >  1:3461928/1‑91 (MQ=255)
                                                                                 aaTATGATGTGCATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGCACtt           >  2:162395/1‑91 (MQ=255)
                                                                                 aaTATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGCACtt           >  2:1007927/1‑91 (MQ=255)
                                                                                  atatGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGCACttt          >  1:2625601/1‑91 (MQ=255)
                                                                                   tatGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGCACtttt         <  2:2606272/91‑1 (MQ=255)
                                                                                    atgatgTGCATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGCACttttt        >  1:302077/1‑91 (MQ=255)
                                                                                     tgatgTGCATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGCACTTTTTa       >  1:1471702/1‑91 (MQ=255)
                                                                                     tgatgCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGCACTTTTTa       >  1:232761/1‑91 (MQ=255)
                                                                                       atgCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGCACTTTTTACt     >  2:2091250/1‑91 (MQ=255)
                                                                                        tgCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGCACTTTTTACTa    >  1:2001515/1‑91 (MQ=255)
                                                                                         gCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGc               <  1:354904/79‑1 (MQ=255)
                                                                                          cacaTCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGCACTTTTTACTAat  <  2:2641711/91‑1 (MQ=255)
                                                                                          |                                                                                          
ACTCTAGTTTTCAAAATAATTGCTATCCCATTTTCATGGTTAGCATAAATTTATGAACTGTAACATTTACGTACTTAGTAAAATATGATGCACATCATATTTGTAACTCATAGAAAATTTTATAATTTTTATCATTATATTTCAACTGAAAATGAGAAACAAAATGGCACTTTTTACTAAT  >  USA300TCH1516_ALE/1870773‑1870953

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