Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,218,722 T→A 9.3% intergenic (+213/+50) 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,7220TA 9.3% 640.4 / 38.9 195intergenic (+213/+50)USA300HOU_RS06065/USA300HOU_RS06070hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (81/95);  new base A (2/16);  total (83/112)
Fisher's exact test for biased strand distribution p-value = 4.83e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 201 total aligned reads displayed.

ATAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTG  >  USA300TCH1516_ALE/1218632‑1218812
                                                                                          |                                                                                           
aTAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCAt                                                                                             <  2:2650747/91‑1 (MQ=255)
aTAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCAt                                                                                             <  2:3860974/91‑1 (MQ=255)
aTAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCAt                                                                                             <  1:1996697/91‑1 (MQ=255)
   aTTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAg                                                                                          >  1:1874390/1‑91 (MQ=255)
     ttGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGaa                                                                                        <  1:336840/91‑1 (MQ=255)
          aCTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAgtag                                                                                   <  2:3797258/91‑4 (MQ=255)
           cTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAgtagt                                                                                  <  1:11246/91‑5 (MQ=255)
           cTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGact                                                                                  >  1:297554/1‑91 (MQ=255)
           cTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGact                                                                                  >  2:2974548/1‑91 (MQ=255)
           cTTAAATGATAATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAgtagt                                                                                  <  2:1886187/91‑5 (MQ=255)
                 tgatgaTACGTTATTTACAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttc                                                                            <  1:3387426/91‑1 (MQ=255)
                  gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttct                                                                           <  1:3944412/91‑1 (MQ=255)
                  gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctt                                                                          <  2:3108491/91‑1 (MQ=25)
                  gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctt                                                                          <  2:2913815/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)
                      atACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttcttatgg                                                                      <  2:1520314/91‑5 (MQ=21)
                         cGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGa                                                                   <  1:2374706/91‑1 (MQ=255)
                         cGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTAtggtga                                                                  <  2:3889438/91‑5 (MQ=18)
                         cGTTATTTAAAAACACGATACTTCGTGTCAATAAACGCATTAAATATTATATAAACTCCTCACCATAAGGTAGGACTACTTTCTTTTGTGa                                                                   <  2:1239441/91‑1 (MQ=255)
                         cGTTATTCAACAACACGATGCTTCTTTCCACAGGACGCATTAAAGAATAAATAAACTCCTCACCATAAGAAAGAAATACTTTCTTTTGTGa                                                                   <  2:3078643/91‑1 (MQ=255)
                            tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgt                                                                <  2:8990/91‑1 (MQ=255)
                            tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGgt                                                                <  1:330213/91‑1 (MQ=17)
                            tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGgt                                                                <  1:3837020/91‑1 (MQ=17)
                            tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGgt                                                                <  2:2922554/91‑1 (MQ=17)
                              tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                              <  1:2711618/91‑1 (MQ=255)
                              tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                              <  1:3225306/91‑1 (MQ=255)
                              tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                              <  2:3849328/91‑1 (MQ=255)
                               ttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGtt                                                             <  1:2209285/91‑1 (MQ=255)
                                tAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGttt                                                            <  1:2772467/91‑1 (MQ=255)
                                tAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGttt                                                            <  1:258738/91‑1 (MQ=255)
                                   aaaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAtt                                                         >  2:3409846/1‑91 (MQ=255)
                                    aaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAttt                                                        >  1:3925925/1‑91 (MQ=255)
                                     acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGGGGTGTTTATTTg                                                       <  1:3054373/91‑1 (MQ=255)
                                     acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTg                                                       <  2:3654933/91‑1 (MQ=255)
                                     acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTAttta                                                       <  1:3448576/91‑2 (MQ=11)
                                      cacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGt                                                      <  2:1558508/91‑1 (MQ=255)
                                      cacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGt                                                      >  1:11082/1‑91 (MQ=255)
                                       acGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTc                                                     >  1:1896132/1‑91 (MQ=255)
                                        cGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCg                                                    >  1:642275/1‑91 (MQ=255)
                                        cGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCg                                                    <  1:902956/91‑1 (MQ=255)
                                         gATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGa                                                   >  2:2592315/1‑91 (MQ=255)
                                          aTACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGAc                                                  >  1:3612745/1‑91 (MQ=255)
                                          aTACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGAc                                                  <  2:2909669/91‑1 (MQ=255)
                                           tACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGCAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACa                                                 >  2:266419/1‑91 (MQ=255)
                                           tACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACa                                                 >  2:2798041/1‑91 (MQ=255)
                                           tACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACa                                                 <  2:3172885/91‑1 (MQ=255)
                                              ttCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACt                                              <  1:716904/91‑1 (MQ=255)
                                               tCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACtt                                             <  2:2079422/91‑1 (MQ=255)
                                               tCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACtt                                             >  2:1650738/1‑91 (MQ=255)
                                                cGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTg                                            <  2:3580994/91‑1 (MQ=255)
                                                 gTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGa                                           >  1:2048889/1‑91 (MQ=255)
                                                 gTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGa                                           >  1:577277/1‑91 (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)
                                                     cAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTAt                                       >  2:1783444/1‑91 (MQ=255)
                                                       aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                     >  1:3010014/1‑91 (MQ=255)
                                                       aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                     <  2:1674367/91‑1 (MQ=255)
                                                       aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                     >  2:1836425/1‑91 (MQ=255)
                                                          aaCGCATTAAATAATAAATAAGCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAAGTGTTTATTTGACGACAACTTGATTATGATTg                                  <  2:483669/91‑1 (MQ=255)
                                                           aCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                                 <  2:2728711/91‑1 (MQ=255)
                                                           aCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                                 >  1:2514477/1‑91 (MQ=255)
                                                            cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                <  2:60259/91‑1 (MQ=255)
                                                            cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                <  1:3916519/91‑1 (MQ=255)
                                                            cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                <  2:2004021/91‑1 (MQ=255)
                                                            cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                <  2:1027553/91‑1 (MQ=255)
                                                            cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                <  2:1177524/91‑1 (MQ=255)
                                                            cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                <  1:3325662/91‑1 (MQ=255)
                                                            cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACGTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                <  1:2280618/91‑1 (MQ=255)
                                                              cATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTc                              <  2:2033507/91‑1 (MQ=255)
                                                               aTTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCa                             >  2:3827661/1‑91 (MQ=255)
                                                                  aaataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttt                          <  1:1197190/91‑1 (MQ=255)
                                                                  aaataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttt                          <  2:1365437/91‑1 (MQ=255)
                                                                   aataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtttt                         <  1:1641885/91‑1 (MQ=255)
                                                                    ataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttttt                        >  1:357487/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)
                                                                          aataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATa                  >  1:472804/1‑91 (MQ=255)
                                                                           ataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATaa                 >  2:245107/1‑91 (MQ=255)
                                                                           ataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATaa                 <  1:2930624/91‑1 (MQ=255)
                                                                            taaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTCTTTGAATacg                >  2:3620574/1‑89 (MQ=255)
                                                                              aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt              >  1:46993/1‑91 (MQ=255)
                                                                              aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt              <  2:77982/91‑1 (MQ=255)
                                                                                cacCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCa            >  2:1581888/1‑91 (MQ=17)
                                                                                 acCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCAt           <  1:2638266/91‑1 (MQ=255)
                                                                                  ccTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATa          >  2:434273/1‑91 (MQ=255)
                                                                                   cTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGTTTTTGTTCATat         >  2:3737702/1‑90 (MQ=255)
                                                                                   cTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATaa         >  1:1275346/1‑91 (MQ=255)
                                                                                     cACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATc       <  1:2456342/91‑1 (MQ=255)
                                                                                       ccATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg     >  1:3577161/1‑91 (MQ=255)
                                                                                       ccATAAGAAAGGACTACTTTCTATTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg     >  1:3709309/1‑91 (MQ=255)
                                                                                        cATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAgg    <  2:2398691/91‑1 (MQ=255)
                                                                                          tAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTTATTTTTTGAATAAGTTCATAATCAGgtg  >  2:1040595/1‑91 (MQ=255)
                                                                                          tAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtg  >  1:1206183/1‑91 (MQ=255)
                                                                                          |                                                                                           
ATAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTG  >  USA300TCH1516_ALE/1218632‑1218812

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: