Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,123,779 T→A 5.5% intergenic (‑39/‑94) USA300HOU_RS10825 ← / → lexA_2 hypothetical protein/LexA repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,123,7790TA 5.5% 916.6 / 26.1 269intergenic (‑39/‑94)USA300HOU_RS10825/lexA_2hypothetical protein/LexA repressor
Reads supporting (aligned to +/- strand):  ref base T (136/117);  new base A (8/8);  total (144/125)
Fisher's exact test for biased strand distribution p-value = 8.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 284 total aligned reads displayed.

TCTTTTTGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAAGGAGGTAAGT  >  USA300TCH1516_ALE/2123691‑2123869
                                                                                        |                                                                                          
tctTTTTGAGTCAATCCGTAATCCGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAGCCTAATATAATACGAAAttt                                                                                          <  2:1437516/91‑1 (MQ=255)
 ctTTTTGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAAtttt                                                                                         <  2:1137249/91‑1 (MQ=255)
 ctTTGTGAGTCTATGCGTGATCAGGCCTCGGCTCTTTTAATGTGTTCATGTTGTCACCGCCTTTCGTAAACCTAGTATAATACGAAAtttt                                                                                         <  2:1364966/91‑1 (MQ=255)
   ttttGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCg                                                                                       <  1:1421374/91‑1 (MQ=255)
    tttGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCATTTCGTAAACCTAATATAATACGAAATTTTCGt                                                                                      >  1:3237490/1‑91 (MQ=255)
      tGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTAt                                                                                    <  1:3503559/91‑1 (MQ=255)
       gAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTAGTAAACCTAATATAATACGAAATTTTCGTAtt                                                                                   >  2:2304129/1‑91 (MQ=255)
        aGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTg                                                                                  >  2:1410937/1‑91 (MQ=255)
        aGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTg                                                                                  <  2:588630/91‑1 (MQ=255)
         gTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGt                                                                                 <  1:175925/91‑1 (MQ=255)
         gTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGt                                                                                 >  1:2811296/1‑91 (MQ=255)
         gTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGt                                                                                 >  2:3142736/1‑91 (MQ=255)
          tCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTc                                                                                >  2:707346/1‑91 (MQ=255)
          tCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTc                                                                                >  1:2337947/1‑91 (MQ=255)
           cAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCa                                                                               >  2:1226452/1‑91 (MQ=255)
           cAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAAATTTCGTAttatat                                                                               <  2:1767732/91‑5 (MQ=255)
             aTCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATATTACGAAATTTTCGTATTGTCAAc                                                                             >  2:2640947/1‑91 (MQ=255)
             aTCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAAc                                                                             >  2:3716177/1‑91 (MQ=255)
             aTCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAAc                                                                             <  2:226536/91‑1 (MQ=255)
                cGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACAtt                                                                          >  2:1195864/1‑91 (MQ=255)
                 gTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTa                                                                         <  2:2156184/91‑1 (MQ=255)
                  tAATCATTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATCTTCGTACTATCAACATTaa                                                                        >  2:33609/1‑91 (MQ=255)
                  tAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTaa                                                                        >  2:2707476/1‑91 (MQ=255)
                  tAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTaa                                                                        >  2:188773/1‑91 (MQ=255)
                  tAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTaa                                                                        >  2:3586169/1‑91 (MQ=255)
                   aaTCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTaaa                                                                       <  1:609473/91‑1 (MQ=255)
                     tCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATa                                                                     >  2:2042316/1‑91 (MQ=255)
                     tCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATa                                                                     >  2:2042935/1‑91 (MQ=255)
                      ccgTCCTCAACTCTTTTACTGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATAc                                                                    >  2:3272303/3‑91 (MQ=255)
                       aGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACg                                                                   >  2:1599309/1‑91 (MQ=255)
                          ccTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGttt                                                                >  1:515775/1‑91 (MQ=255)
                              aaCTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGttttttt                                                            <  1:3086045/91‑1 (MQ=255)
                              aaCTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGttttttt                                                            >  1:1644606/1‑91 (MQ=255)
                               aCTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTc                                                           >  1:505098/1‑91 (MQ=255)
                                ctctTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCg                                                          >  1:1024940/1‑91 (MQ=255)
                                atctTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCg                                                          <  2:3928672/90‑1 (MQ=255)
                                  ctTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTa                                                        <  1:3127700/91‑1 (MQ=255)
                                  ctTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTa                                                        >  1:1164069/1‑91 (MQ=255)
                                   ttttAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTaa                                                       <  1:2700282/91‑1 (MQ=255)
                                   ttttAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTACgaaa                                                       <  1:2794869/91‑4 (MQ=255)
                                    tttAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTaag                                                      <  1:940545/91‑2 (MQ=255)
                                    tttAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTaaa                                                      >  2:2516910/1‑91 (MQ=255)
                                     ttAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTaata                                                     >  2:1593839/1‑89 (MQ=255)
                                     ttAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTaaaa                                                     <  2:812286/91‑1 (MQ=255)
                                     ttAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTaaaa                                                     >  2:433323/1‑91 (MQ=255)
                                     ttAATGATTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTcaaa                                                     >  1:3068466/1‑91 (MQ=255)
                                       aaTGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTaaaaac                                                   <  1:1215371/91‑2 (MQ=255)
                                       aaTGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTaaaaaa                                                   >  2:2968018/1‑91 (MQ=255)
                                       aaTGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTaaaaaa                                                   <  1:3821427/91‑1 (MQ=255)
                                       aaTGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTACGAaaaaaa                                                   <  1:1197843/91‑1 (MQ=255)
                                        aTGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAAc                                                  <  2:361304/91‑1 (MQ=255)
                                         tGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACt                                                 <  1:1102659/91‑1 (MQ=255)
                                          gTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACtt                                                >  2:3436318/1‑91 (MQ=255)
                                          gTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTAACAAAAAACtt                                                <  2:2602747/91‑1 (MQ=255)
                                            tttCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTa                                              >  1:1813809/1‑91 (MQ=255)
                                            tttCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTa                                              >  1:2907944/1‑91 (MQ=255)
                                            tttCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTa                                              >  1:331422/1‑91 (MQ=255)
                                             ttCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTAc                                             >  1:447729/1‑91 (MQ=255)
                                             ttCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTAc                                             >  1:1590772/1‑91 (MQ=255)
                                             ttCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTAc                                             >  1:3653456/1‑91 (MQ=255)
                                              tCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACt                                            <  1:2652161/91‑1 (MQ=255)
                                              tCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACt                                            >  1:1461939/1‑91 (MQ=255)
                                               cATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCTTAAAAAACTTTACta                                           >  2:911864/1‑91 (MQ=255)
                                                 tgttgtCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatg                                         >  1:533298/1‑91 (MQ=255)
                                                   ttgtCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgat                                       <  2:862226/91‑1 (MQ=255)
                                                   ttgtCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgat                                       >  1:1862488/1‑91 (MQ=255)
                                                     gtCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgatat                                     <  1:2569126/91‑1 (MQ=255)
                                                     gtCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgatat                                     >  1:31766/1‑91 (MQ=255)
                                                      tCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAATAACTTTACtatgatatg                                    >  2:284302/1‑91 (MQ=255)
                                                      tCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgatatg                                    <  2:3802636/91‑1 (MQ=255)
                                                      tCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgatatg                                    <  2:1025714/91‑1 (MQ=255)
                                                      tCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgatatg                                    <  1:3882076/91‑1 (MQ=255)
                                                        aCCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGaa                                  <  1:3852474/91‑1 (MQ=255)
                                                         ccGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGaaa                                 <  1:954844/91‑1 (MQ=255)
                                                         ccGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGaaa                                 >  2:324797/1‑91 (MQ=255)
                                                            ccTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAAtt                              <  2:1076529/91‑1 (MQ=255)
                                                            ccTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAAtt                              <  1:823337/91‑1 (MQ=255)
                                                            ccTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAAtt                              >  1:373847/1‑91 (MQ=255)
                                                                tCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGt                          <  1:1250883/91‑1 (MQ=255)
                                                                 cGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGta                         >  2:2758619/1‑91 (MQ=255)
                                                                 cGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGta                         <  1:3624788/91‑1 (MQ=255)
                                                                 cGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGta                         >  2:1303737/1‑91 (MQ=255)
                                                                 cGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGta                         <  2:1565651/91‑1 (MQ=255)
                                                                  gTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTCCTATGATATGAAAATTTCGtat                        >  2:2399464/1‑91 (MQ=255)
                                                                     aaCCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTataat                     <  1:2416836/91‑1 (MQ=255)
                                                                     aaCCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTataat                     <  1:2417343/91‑1 (MQ=255)
                                                                        cTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAg                  >  1:639808/1‑91 (MQ=255)
                                                                         taatataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGa                 >  2:2164258/1‑91 (MQ=255)
                                                                           atataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATATGaaa               <  2:1985292/91‑1 (MQ=255)
                                                                           atataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGaaa               <  1:1005554/91‑1 (MQ=255)
                                                                           atataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGaaa               >  1:3601122/1‑91 (MQ=255)
                                                                             ataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGaaaaa             <  2:44128/91‑1 (MQ=255)
                                                                               aataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTTGTTTTTAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAag           >  2:1942352/1‑91 (MQ=255)
                                                                                ataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAagg          <  1:2161996/91‑1 (MQ=255)
                                                                                   cGAAAATTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAaggagg       >  1:3266999/1‑91 (MQ=255)
                                                                                   cGAAAATTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAaggagg       >  2:3006144/1‑91 (MQ=255)
                                                                                       aTTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAAGGAGgtaag   <  1:3903712/91‑1 (MQ=255)
                                                                                       aTTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCAAAAAAAAAAAAAAAAGGAGgtaag   <  2:1313544/91‑1 (MQ=255)
                                                                                        ttttCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAAGGAGgtaagt  <  1:2978933/91‑1 (MQ=255)
                                                                                        ttttCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAAGGAGgtaagt  >  1:402062/1‑91 (MQ=255)
                                                                                        |                                                                                          
TCTTTTTGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAAGGAGGTAAGT  >  USA300TCH1516_ALE/2123691‑2123869

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: