Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,218,734 A→C 11.7% intergenic (+225/+38) 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,7340AC11.7% 604.7 / 68.2 202intergenic (+225/+38)USA300HOU_RS06065/USA300HOU_RS06070hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (89/89);  new base C (8/16);  total (97/105)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

Only 100 of 225 total aligned reads displayed.

TTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAG  >  USA300TCH1516_ALE/1218644‑1218824
                                                                                          |                                                                                            
  aaaTGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGactact                                                                                            <  2:452628/91‑1 (MQ=255)
      gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttct                                                                                        >  2:2932206/1‑91 (MQ=255)
      gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttct                                                                                        <  1:1200044/91‑1 (MQ=255)
      gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctt                                                                                       <  2:2913894/91‑1 (MQ=25)
       atgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttctt                                                                                       <  1:3111272/91‑1 (MQ=255)
       atgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctta                                                                                      <  1:798432/91‑2 (MQ=21)
        tgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAGCGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttcttat                                                                                     <  2:2038579/91‑3 (MQ=21)
         gatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCtttt                                                                                    <  1:124588/91‑1 (MQ=255)
          atACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttcttatgg                                                                                   <  2:1520314/91‑5 (MQ=21)
             cGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGa                                                                               <  1:2374706/91‑1 (MQ=255)
             cGTTATTTAAAAACACGATACTTCGTGTCAATAAACGCATTAAATATTATATAAACTCCTCACCATAAGGTAGGACTACTTTCTTTTGTGa                                                                               <  2:1239441/91‑1 (MQ=255)
               ttATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAgg                                                                             >  2:745898/1‑91 (MQ=255)
               ttATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAgg                                                                             <  1:2609939/91‑1 (MQ=255)
                tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGgt                                                                            <  1:3837020/91‑1 (MQ=17)
                tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGgt                                                                            <  1:330213/91‑1 (MQ=17)
                  tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                                          <  2:3849328/91‑1 (MQ=255)
                   ttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGtt                                                                         <  1:2209285/91‑1 (MQ=255)
                    tAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGttt                                                                        <  1:3349698/91‑1 (MQ=255)
                       aaaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAtt                                                                     <  2:2671691/91‑1 (MQ=255)
                        aaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAttt                                                                    >  1:3925925/1‑91 (MQ=255)
                         acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTg                                                                   <  2:3654882/91‑1 (MQ=255)
                         acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTg                                                                   >  1:2651647/1‑91 (MQ=255)
                         acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTAttta                                                                   <  1:3448576/91‑2 (MQ=11)
                         acacGAGACTTCGTTTCAATGGGCGCATTAAATAATGAGTAAACGTCTCACCAGAAGAAAGGGCTACGTTCTTTTGTGAGGTGTTTATTTg                                                                   <  2:3656257/91‑1 (MQ=255)
                          cacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGt                                                                  >  1:11082/1‑91 (MQ=255)
                          cacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGt                                                                  <  2:1558508/91‑1 (MQ=255)
                           acGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTc                                                                 >  1:1895818/1‑91 (MQ=255)
                            cGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCg                                                                >  1:642275/1‑91 (MQ=255)
                            cGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCg                                                                <  1:902956/91‑1 (MQ=255)
                            cGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCg                                                                <  1:433252/91‑1 (MQ=255)
                             gATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGa                                                               >  2:2592315/1‑91 (MQ=255)
                              aTACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGAc                                                              >  1:3612745/1‑91 (MQ=255)
                               tACTTCGTTTCAATGAACGCATTAAATAATAAATACACACCTCACCATACGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACa                                                             >  2:3165268/1‑91 (MQ=255)
                                  ttCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACt                                                          <  1:716904/91‑1 (MQ=255)
                                   tCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACtt                                                         >  2:1650738/1‑91 (MQ=255)
                                   tCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACtt                                                         <  2:2079422/91‑1 (MQ=255)
                                    cGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTg                                                        <  2:3580994/91‑1 (MQ=255)
                                     gTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGa                                                       >  1:3088542/1‑91 (MQ=255)
                                       ttCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGAtt                                                     >  1:1010180/1‑91 (MQ=255)
                                       ttCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGAtt                                                     <  2:2492768/91‑1 (MQ=255)
                                        tCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTa                                                    <  1:1831263/91‑1 (MQ=255)
                                          aaTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATg                                                  >  1:515841/1‑91 (MQ=255)
                                           aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                                 <  2:1674367/91‑1 (MQ=255)
                                           aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                                 >  2:1836425/1‑91 (MQ=255)
                                           aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                                 >  1:3010014/1‑91 (MQ=255)
                                            tGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGAt                                                >  2:821556/1‑91 (MQ=255)
                                             gAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGAtt                                               <  2:3106492/91‑1 (MQ=255)
                                              aaCGCATTAAATAATAAATAAGCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAAGTGTTTATTTGACGACAACTTGATTATGATTg                                              <  2:483669/91‑1 (MQ=255)
                                               aCGCATTAAATAATAAATAAACACCTCGCCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                                             <  2:2729102/91‑1 (MQ=255)
                                               aCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                                             >  1:115873/1‑91 (MQ=255)
                                                cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                            <  2:1027553/91‑1 (MQ=255)
                                                cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                            <  2:2004021/91‑1 (MQ=255)
                                                cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                            <  1:3325662/91‑1 (MQ=255)
                                                cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                            <  1:3916519/91‑1 (MQ=255)
                                                cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACGTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                            <  1:2280618/91‑1 (MQ=255)
                                                  cATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTc                                          <  1:3079392/91‑1 (MQ=255)
                                                   aTTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCa                                         >  2:3827661/1‑91 (MQ=255)
                                                      aaataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttt                                      <  2:1365437/91‑1 (MQ=255)
                                                       aataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtttt                                     <  1:1641885/91‑1 (MQ=255)
                                                        ataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttttt                                    >  1:357487/1‑91 (MQ=255)
                                                        ataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttttt                                    >  1:1933328/1‑91 (MQ=255)
                                                         taataaataaACTCCTCACCATGAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtttttt                                   <  2:2131532/91‑1 (MQ=255)
                                                         taataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtttttt                                   <  1:1379458/91‑1 (MQ=255)
                                                           ataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGa                                 <  1:1612640/91‑1 (MQ=255)
                                                            taaataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGaa                                <  2:2626018/91‑1 (MQ=255)
                                                            taaataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGaa                                <  2:3382182/91‑1 (MQ=255)
                                                               ataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATaa                             >  2:245107/1‑91 (MQ=255)
                                                               ataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATaa                             <  1:2930751/91‑1 (MQ=255)
                                                                taaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTCTTTGAATacg                            >  2:3620574/1‑89 (MQ=255)
                                                                     acCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCAt                       >  2:2913088/1‑91 (MQ=255)
                                                                      ccTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATa                      <  2:1144837/91‑1 (MQ=255)
                                                                         cACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATc                   <  1:2456342/91‑1 (MQ=255)
                                                                         cACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATc                   >  1:255728/1‑91 (MQ=255)
                                                                          aCCATAAGAAAGGACTACTTTCTTTTGTGAGGCGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTtcatcatca                  >  2:3875809/1‑91 (MQ=255)
                                                                           ccATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg                 >  1:3577161/1‑91 (MQ=255)
                                                                           ccATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg                 >  1:2022192/1‑91 (MQ=255)
                                                                           ccATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg                 <  1:470464/91‑1 (MQ=255)
                                                                           acaaaaGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTTGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCa                  >  2:2012257/5‑91 (MQ=18)
                                                                            caaaaGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg                 >  2:422806/4‑91 (MQ=21)
                                                                            cATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAgg                <  2:2808431/91‑1 (MQ=255)
                                                                            cATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAgg                >  1:2888643/1‑91 (MQ=255)
                                                                              tAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTTATTTTTTGAATAAGTTCATAATCAGgtg              >  2:1040595/1‑91 (MQ=255)
                                                                               aaGAAAGTGCTACTTTCTTTTGTGAGGTTTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtgt             <  2:938597/91‑1 (MQ=255)
                                                                               aaGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtg              >  2:3193226/1‑91 (MQ=21)
                                                                               aaGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtgt             <  1:3371915/91‑1 (MQ=255)
                                                                                aGAAAGGACTACTTTCTCTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATCGTTTCATTTCTTGAATAAGTTCATAATCAGGTGTa            >  2:1918019/1‑91 (MQ=255)
                                                                                aGAAAGGACTACGTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTa            >  1:3030287/1‑91 (MQ=255)
                                                                                 gAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAg           <  2:3733675/91‑1 (MQ=255)
                                                                                 gAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAg           >  2:107104/1‑91 (MQ=255)
                                                                                  aaaGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGc          >  2:3008547/1‑91 (MQ=255)
                                                                                  aaaGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGc          >  2:2380536/1‑91 (MQ=255)
                                                                                   aaGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCa         >  2:163644/1‑91 (MQ=255)
                                                                                   aaGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCa         >  2:1760552/1‑91 (MQ=255)
                                                                                     ggACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATa       >  2:51291/1‑91 (MQ=255)
                                                                                      gACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATAc      >  2:2710342/1‑91 (MQ=255)
                                                                                       actactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGCTTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACa     >  2:888746/1‑91 (MQ=255)
                                                                                       actactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACa     <  1:966501/91‑1 (MQ=255)
                                                                                       actactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACa     >  1:3721111/1‑91 (MQ=255)
                                                                                       actactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACa     >  1:2740194/1‑91 (MQ=255)
                                                                                        ctactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACaa    <  1:3122804/91‑1 (MQ=255)
                                                                                          |                                                                                            
TTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAG  >  USA300TCH1516_ALE/1218644‑1218824

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