Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,218,719 Δ1 bp 9.3% intergenic (+210/+53) 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,7190C. 9.3% 659.3 / 54.8 194intergenic (+210/+53)USA300HOU_RS06065/USA300HOU_RS06070hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (79/97);  new base . (2/16);  total (81/113)
Fisher's exact test for biased strand distribution p-value = 5.27e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

Only 100 of 202 total aligned reads displayed.

AAGATAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAG  >  USA300TCH1516_ALE/1218629‑1218809
                                                                                          |                                                                                           
aaGATAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCAc                                                                                             <  1:2759825/91‑1 (MQ=255)
  gATAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAGCGCATTAAATAATAAATAAACACCTCACCa                                                                                           >  2:2811165/1‑91 (MQ=255)
  gATAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCa                                                                                           >  1:1091055/1‑91 (MQ=255)
   aTAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCAt                                                                                          <  2:3860974/91‑1 (MQ=255)
   aTAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCAt                                                                                          <  1:1996697/91‑1 (MQ=255)
      aTTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTCAATAATAAATAAACACCTCACCATAAg                                                                                       >  2:262713/1‑91 (MQ=255)
      aTTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAg                                                                                       >  1:1874390/1‑91 (MQ=255)
        ttGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGaa                                                                                     <  1:336840/91‑1 (MQ=255)
             aCTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAgtag                                                                                <  2:3797258/91‑4 (MQ=255)
              cTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGact                                                                               >  2:2974548/1‑91 (MQ=255)
              cTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGact                                                                               >  1:297554/1‑91 (MQ=255)
               ttAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGacta                                                                              <  2:1448929/91‑1 (MQ=255)
                 aaaTGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGactact                                                                            <  1:897259/91‑1 (MQ=255)
                 aaaTGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGactact                                                                            <  2:1567803/91‑1 (MQ=255)
                    tgatgaTACGTTATTTACAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttc                                                                         <  1:3387426/91‑1 (MQ=255)
                     gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttct                                                                        <  1:3944412/91‑1 (MQ=255)
                     gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttct                                                                        <  1:1200044/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                                                                      <  2:1249390/91‑2 (MQ=21)
                       tgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttcttat                                                                     <  2:2036942/91‑3 (MQ=21)
                          tACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTtgt                                                                  <  2:1649225/91‑1 (MQ=255)
                            cGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGa                                                                <  1:2374706/91‑1 (MQ=255)
                            cGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTAtggtga                                                               <  2:3889438/91‑5 (MQ=18)
                            cGTTATTCAACAACACGATGCTTCTTTCCACAGGACGCATTAAAGAATAAATAAACTCCTCACCATAAGAAAGAAATACTTTCTTTTGTGa                                                                <  2:3078643/91‑1 (MQ=255)
                               tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgt                                                             <  2:8990/91‑1 (MQ=255)
                               tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgt                                                             >  2:2658682/1‑91 (MQ=255)
                               tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGgt                                                             <  1:330213/91‑1 (MQ=17)
                               tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGgt                                                             <  2:2922554/91‑1 (MQ=17)
                               tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGgt                                                             <  1:3837020/91‑1 (MQ=17)
                                aTTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtg                                                            <  2:2006080/91‑1 (MQ=255)
                                 tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                           <  1:3225306/91‑1 (MQ=255)
                                 tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                           <  1:2711618/91‑1 (MQ=255)
                                 tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                           <  2:3849328/91‑1 (MQ=255)
                                  ttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGtt                                                          <  1:2209285/91‑1 (MQ=255)
                                   tAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGttt                                                         <  1:77144/91‑1 (MQ=255)
                                   tAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGttt                                                         <  1:258738/91‑1 (MQ=255)
                                      aaaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAtt                                                      >  2:3409846/1‑91 (MQ=255)
                                      aaaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAtt                                                      <  2:2671691/91‑1 (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:2582816/1‑91 (MQ=255)
                                         cacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGt                                                   <  2:1184768/91‑1 (MQ=255)
                                          acGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTc                                                  >  1:1896132/1‑91 (MQ=255)
                                           cGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCg                                                 >  1:642275/1‑91 (MQ=255)
                                             aTACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGAc                                               <  2:2909669/91‑1 (MQ=255)
                                             aTACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGAc                                               >  1:3612745/1‑91 (MQ=255)
                                              tACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACa                                              <  2:3172885/91‑1 (MQ=255)
                                                  tCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGAGAGGTGTTTATTTGTCGACAACtt                                          >  2:1650773/1‑91 (MQ=255)
                                                   cGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTg                                         >  2:260586/1‑91 (MQ=255)
                                                   cGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTg                                         <  2:3580994/91‑1 (MQ=255)
                                                    gTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGa                                        >  1:577277/1‑91 (MQ=255)
                                                    gTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGa                                        >  1:2048889/1‑91 (MQ=255)
                                                    gTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGa                                        >  1:3088542/1‑91 (MQ=255)
                                                      ttCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGAtt                                      >  1:1010180/1‑91 (MQ=255)
                                                       tCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTa                                     <  1:1831263/91‑1 (MQ=255)
                                                         aaTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATg                                   >  2:2278476/1‑91 (MQ=255)
                                                          aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                  >  2:1836424/1‑91 (MQ=255)
                                                          aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                  >  1:3010014/1‑91 (MQ=255)
                                                            gAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGAtt                                <  2:1682575/91‑1 (MQ=255)
                                                             aaCGCATTAAATAATAAATAAGCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAAGTGTTTATTTGACGACAACTTGATTATGATTg                               <  2:483669/91‑1 (MQ=255)
                                                              aCGCATTAAATAATAAATAAACACCTCGCCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                              <  2:2729102/91‑1 (MQ=255)
                                                              aCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                              >  1:2514477/1‑91 (MQ=255)
                                                              aCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                              >  1:115873/1‑91 (MQ=255)
                                                               cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                             <  1:3325662/91‑1 (MQ=255)
                                                               cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                             <  2:60259/91‑1 (MQ=255)
                                                               cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                             <  1:3916519/91‑1 (MQ=255)
                                                               cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                             <  2:1027553/91‑1 (MQ=255)
                                                               cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                             <  2:2128678/91‑1 (MQ=255)
                                                               cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACGTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                             <  1:2280618/91‑1 (MQ=255)
                                                                  aTTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCa                          >  2:3827661/1‑91 (MQ=255)
                                                                  aTTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCa                          >  2:1837502/1‑91 (MQ=255)
                                                                      aataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtttt                      <  1:1641885/91‑1 (MQ=255)
                                                                       ataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttttt                     >  1:357487/1‑91 (MQ=255)
                                                                        taataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtttttt                    <  1:1379458/91‑1 (MQ=255)
                                                                         aataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTg                   <  2:2073925/91‑1 (MQ=255)
                                                                          ataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGa                  <  1:1612640/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)
                                                                               taaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTCTTTGAATacg             >  2:3620574/1‑89 (MQ=255)
                                                                                 aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt           >  1:46993/1‑91 (MQ=255)
                                                                                 aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt           <  2:1963755/91‑1 (MQ=255)
                                                                                 aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt           <  2:77982/91‑1 (MQ=255)
                                                                                 aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt           <  1:701587/91‑1 (MQ=255)
                                                                                    acCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCAt        <  1:2638266/91‑1 (MQ=255)
                                                                                     ccTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATa       >  2:434273/1‑91 (MQ=255)
                                                                                     ccTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATa       <  2:1144837/91‑1 (MQ=255)
                                                                                      cTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGTTTTTGTTCATat      >  2:3737702/1‑90 (MQ=255)
                                                                                      cTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATaa      >  1:1275346/1‑91 (MQ=255)
                                                                                        cACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATc    <  1:2456342/91‑1 (MQ=255)
                                                                                         aCCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCa   >  2:2489810/1‑91 (MQ=255)
                                                                                         aCCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCa   <  2:1531724/91‑1 (MQ=255)
                                                                                          ccATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg  >  1:3577161/1‑91 (MQ=255)
                                                                                          ccATAAGAAAGGACTACTTTCTATTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg  >  1:3709309/1‑91 (MQ=255)
                                                                                          |                                                                                           
AAGATAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAG  >  USA300TCH1516_ALE/1218629‑1218809

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: