Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,218,741:1 +A 5.6% intergenic (+232/+31) USA300HOU_RS06065 → / ← USA300HOU_RS06070 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,218,7411.A 5.6% 832.6 / 39.1 197intergenic (+232/+31)USA300HOU_RS06065/USA300HOU_RS06070hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (97/89);  new base A (7/4);  total (104/93)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 235 total aligned reads displayed.

GATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGTTTTTT  >  USA300TCH1516_ALE/1218650‑1218830
                                                                                            |                                                                                         
gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctt                                                                                            <  2:2913894/91‑1 (MQ=25)
gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctt                                                                                            <  2:3108491/91‑1 (MQ=25)
 atgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctta                                                                                           <  2:1249390/91‑2 (MQ=21)
   gatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCtttt                                                                                         <  1:3486335/91‑1 (MQ=255)
   gatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttcttatg                                                                                         <  1:2222774/91‑4 (MQ=21)
   gatACGTAATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttcttatg                                                                                         <  2:3515513/91‑4 (MQ=18)
    atACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttcttatgg                                                                                        <  2:1520314/91‑5 (MQ=21)
       cGTTATTCAACAACACGATGCTTCTTTCCACAGGACGCATTAAAGAATAAATAAACTCCTCACCATAAGAAAGAAATACTTTCTTTTGTGa                                                                                     <  2:3078643/91‑1 (MQ=255)
         ttATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAgg                                                                                   >  2:745898/1‑91 (MQ=255)
          tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgt                                                                                  >  2:2658682/1‑91 (MQ=255)
          tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGgt                                                                                  <  1:3837020/91‑1 (MQ=17)
           aTTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtg                                                                                 <  2:2006080/91‑1 (MQ=255)
           aTTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtg                                                                                 >  1:3919822/1‑91 (MQ=255)
            tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                                                <  1:2711618/91‑1 (MQ=255)
            tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                                                <  2:3849328/91‑1 (MQ=255)
              tAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGttt                                                                              <  1:258738/91‑1 (MQ=255)
              tAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGttt                                                                              <  1:77144/91‑1 (MQ=255)
              tAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGttt                                                                              <  1:2772467/91‑1 (MQ=255)
                   acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGGGGTGTTTATTTg                                                                         <  1:3054373/91‑1 (MQ=255)
                   acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTg                                                                         <  2:3654882/91‑1 (MQ=255)
                   acacGAGACTTCGTTTCAATGGGCGCATTAAATAATGAGTAAACGTCTCACCAGAAGAAAGGGCTACGTTCTTTTGTGAGGTGTTTATTTg                                                                         <  2:3656257/91‑1 (MQ=255)
                    cacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGt                                                                        >  2:2582816/1‑91 (MQ=255)
                    cacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGt                                                                        <  2:1558508/91‑1 (MQ=255)
                     acGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTc                                                                       >  1:1895818/1‑91 (MQ=255)
                       gATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGa                                                                     >  1:2341370/1‑91 (MQ=255)
                         tACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACa                                                                   <  2:3172885/91‑1 (MQ=255)
                          aCTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACaa                                                                  >  1:1698586/1‑91 (MQ=255)
                            ttCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACt                                                                <  1:716904/91‑1 (MQ=255)
                             tCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACtt                                                               >  2:1650738/1‑91 (MQ=255)
                              cGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTg                                                              >  2:260586/1‑91 (MQ=255)
                                 ttCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGAtt                                                           >  1:1010180/1‑91 (MQ=255)
                                   cAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTAt                                                         >  2:1783444/1‑91 (MQ=255)
                                    aaTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATg                                                        >  1:3939846/1‑91 (MQ=255)
                                     aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                                       >  2:1836425/1‑91 (MQ=255)
                                     aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                                       <  2:1674367/91‑1 (MQ=255)
                                     aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                                       >  1:3010014/1‑91 (MQ=255)
                                      tGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGAt                                                      >  2:821556/1‑91 (MQ=255)
                                         aCGCATTAAATAATAAATAAACACCTCGCCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                                                   <  2:2729102/91‑1 (MQ=255)
                                         aCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                                                   >  1:2514477/1‑91 (MQ=255)
                                          cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                                  <  2:2128678/91‑1 (MQ=255)
                                          cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                                  <  1:3099512/91‑1 (MQ=255)
                                          cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                                  <  1:1570050/91‑1 (MQ=255)
                                          cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                                  <  2:1177524/91‑1 (MQ=255)
                                          cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                                  <  2:1027553/91‑1 (MQ=255)
                                          cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                                  <  2:60259/91‑1 (MQ=255)
                                            cATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTc                                                <  2:2033507/91‑1 (MQ=255)
                                             aTTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCa                                               >  2:3827661/1‑91 (MQ=255)
                                                  ataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttttt                                          >  1:1933328/1‑91 (MQ=255)
                                                  ataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttttt                                          >  1:357499/1‑91 (MQ=255)
                                                    aataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTg                                        <  2:3612745/91‑1 (MQ=255)
                                                    aataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTg                                        <  2:2073925/91‑1 (MQ=255)
                                                      taaataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGaa                                      <  2:3382182/91‑1 (MQ=255)
                                                        aataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATa                                    >  1:472804/1‑91 (MQ=255)
                                                         ataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATaa                                   <  1:2930624/91‑1 (MQ=255)
                                                         ataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATaa                                   >  2:245107/1‑91 (MQ=255)
                                                          taaaCACCTTACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGATTAAg                                  >  1:1121913/1‑91 (MQ=255)
                                                            aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt                                >  1:46993/1‑91 (MQ=255)
                                                            aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt                                >  2:487213/1‑91 (MQ=255)
                                                            aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt                                <  2:1963755/91‑1 (MQ=255)
                                                               acCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCAt                             >  2:2913088/1‑91 (MQ=255)
                                                               acCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCAt                             <  1:2638266/91‑1 (MQ=255)
                                                                 cTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATaa                           >  1:1275346/1‑91 (MQ=255)
                                                                  tCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAAt                          <  1:1061682/91‑1 (MQ=255)
                                                                  tCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAAt                          <  1:1632972/91‑1 (MQ=255)
                                                                   cACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATc                         <  1:2456342/91‑1 (MQ=255)
                                                                    aCCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCa                        >  1:3302163/1‑91 (MQ=255)
                                                                    aCCATAAGAAAGGACTACTTTCTTTTGTGAGGCGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTtcatcatca                        >  2:3875809/1‑91 (MQ=255)
                                                                     ccATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg                       >  1:3578972/1‑91 (MQ=255)
                                                                     ccATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg                       >  1:2022192/1‑91 (MQ=255)
                                                                     ccATAAGAAAGGACTACTTTCTATTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg                       >  1:3709309/1‑91 (MQ=255)
                                                                      caaaaGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGCCAACTTGATTCTGATTGTTTCATTTTTTGAATAAGTTCATAATCAg                       >  2:1074534/4‑91 (MQ=18)
                                                                      caaaaGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg                       >  2:422806/4‑91 (MQ=21)
                                                                      cATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAgg                      <  2:2398691/91‑1 (MQ=255)
                                                                       aaaaGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAgg                      >  2:1056969/3‑91 (MQ=21)
                                                                        tAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtg                    >  1:1206183/1‑91 (MQ=255)
                                                                         aaGAAAGTGCTACTTTCTTTTGTGAGGTTTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtgt                   <  2:938597/91‑1 (MQ=255)
                                                                         aaGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtgt                   <  2:2906614/91‑1 (MQ=255)
                                                                         aaGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtgt                   >  1:987422/1‑91 (MQ=255)
                                                                         aaGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtgt                   <  1:3371915/91‑1 (MQ=255)
                                                                            aaaGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGc                >  2:3008547/1‑91 (MQ=255)
                                                                            aaaGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGc                >  2:2380536/1‑91 (MQ=255)
                                                                             aaGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCa               >  2:163644/1‑91 (MQ=255)
                                                                               ggACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATa             >  1:3088054/1‑91 (MQ=255)
                                                                                gACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATAc            >  2:2710342/1‑91 (MQ=255)
                                                                                 actactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGCTTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACa           >  2:888746/1‑91 (MQ=255)
                                                                                 actactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACa           <  1:966470/91‑1 (MQ=255)
                                                                                  gtcctttctttTGTGAGGTGTTTATTTGTCGACGACTTGATTATGAGTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACaa          >  1:3150237/4‑91 (MQ=255)
                                                                                  ctactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACaa          >  1:686629/1‑91 (MQ=255)
                                                                                  ctactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACaa          <  1:3122804/91‑1 (MQ=255)
                                                                                  ctactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTAGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACaa          >  1:3339803/1‑91 (MQ=255)
                                                                                   tactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACaaa         >  2:3703768/1‑91 (MQ=255)
                                                                                    actTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAg        >  2:458993/1‑91 (MQ=255)
                                                                                     ctttctttTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGt       <  1:3724483/91‑1 (MQ=255)
                                                                                     ctttctttTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGt       >  1:3167288/1‑91 (MQ=255)
                                                                                     ctttctttTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGt       <  1:521358/91‑1 (MQ=255)
                                                                                        tctttTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGtttt    <  1:12224/91‑1 (MQ=255)
                                                                                         tttttGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGttttt   <  2:3035972/90‑1 (MQ=255)
                                                                                         ctttTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGttttt   <  1:2184910/91‑1 (MQ=255)
                                                                                         ctttTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGttttt   >  2:2175139/1‑91 (MQ=255)
                                                                                         ctttTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGttttt   >  2:2680315/1‑91 (MQ=255)
                                                                                            |                                                                                         
GATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGTTTTTT  >  USA300TCH1516_ALE/1218650‑1218830

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: