Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,559,831 G→A 11.5% intergenic (+25/+98) bcr_2 → / ← USA300HOU_RS13190 Bicyclomycin resistance protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,559,8310GA11.5% 598.9 / 42.0 204intergenic (+25/+98)bcr_2/USA300HOU_RS13190Bicyclomycin resistance protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (104/75);  new base A (3/21);  total (107/97)
Fisher's exact test for biased strand distribution p-value = 1.98e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 207 total aligned reads displayed.

GCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATT  >  USA300TCH1516_ALE/2559741‑2559919
                                                                                          |                                                                                        
gCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAg                                                                                          >  1:3467032/1‑91 (MQ=255)
gCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAg                                                                                          <  1:2805189/91‑1 (MQ=255)
 ccATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCTCATAACACTTCAACATAATTAGAACCCTAGc                                                                                         >  1:3181955/1‑91 (MQ=255)
  cATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCa                                                                                        >  2:1161647/1‑91 (MQ=255)
   aTTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCaa                                                                                       >  1:1428174/1‑91 (MQ=255)
    ttCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCaaa                                                                                      >  2:2113491/1‑91 (MQ=255)
    ttCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCaaa                                                                                      >  1:2684274/1‑91 (MQ=255)
      cTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGTAAAGa                                                                                    <  1:595656/91‑1 (MQ=255)
      cTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGa                                                                                    >  2:1196259/1‑91 (MQ=255)
       tATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGat                                                                                   >  1:5179/1‑91 (MQ=255)
        aTTAGTCAGTCTACAAATCATTTACTTTAAAGTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGata                                                                                  <  2:1280764/91‑1 (MQ=255)
        aTTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAACGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGata                                                                                  <  2:3481706/91‑1 (MQ=255)
        aTTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGata                                                                                  >  1:614118/1‑91 (MQ=255)
        aTTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGata                                                                                  >  1:2426711/1‑91 (MQ=255)
         ttAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGatat                                                                                 <  2:689872/91‑1 (MQ=255)
          tAGTCAGTCTACCAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatc                                                                                >  1:3023122/1‑91 (MQ=255)
           agtcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatct                                                                               <  1:3284926/91‑1 (MQ=255)
            gtcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatctc                                                                              >  2:3193749/1‑90 (MQ=255)
            gtcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatcta                                                                              >  1:1943076/1‑91 (MQ=255)
             tcagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatctat                                                                             <  2:768200/91‑1 (MQ=255)
             tcagtcTACAAATCACTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACGTCAACATAATTAGAACCCTCGCAACGAtatctat                                                                             <  2:1426144/91‑1 (MQ=255)
              cagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatctatc                                                                            >  1:2341925/1‑91 (MQ=255)
               agtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATagatatct                                                                           <  2:2649706/91‑1 (MQ=38)
               agtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATatctatct                                                                           <  2:2608449/91‑1 (MQ=255)
               agtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATatctatct                                                                           >  2:1724276/1‑91 (MQ=255)
                gtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCtt                                                                          >  1:3239042/1‑91 (MQ=255)
                 tcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCttt                                                                         <  1:2733959/91‑1 (MQ=38)
                   tACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgc                                                                       <  2:3125098/91‑2 (MQ=38)
                   tACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgc                                                                       <  2:1908064/91‑2 (MQ=37)
                   tACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgc                                                                       <  1:3870316/91‑2 (MQ=37)
                   tACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgc                                                                       <  2:3147748/91‑2 (MQ=37)
                   tACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGt                                                                       <  1:365674/91‑1 (MQ=255)
                    aCAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgct                                                                      <  2:2808278/91‑3 (MQ=37)
                     cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgcta                                                                     <  1:1499121/91‑4 (MQ=37)
                     cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATCTCTATCTTTGTCa                                                                     >  2:1490481/1‑91 (MQ=255)
                     cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCa                                                                     >  1:3790848/1‑91 (MQ=255)
                     cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCa                                                                     >  2:1098600/1‑91 (MQ=255)
                     cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCa                                                                     >  2:2180429/1‑91 (MQ=255)
                      aaaTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAAGTAGAACCCTAGCAAAGATATCTATCTTTGTCAg                                                                    <  2:3577280/91‑1 (MQ=255)
                        aTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAggg                                                                  >  1:2786493/1‑91 (MQ=255)
                         tCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGt                                                                 <  2:41694/91‑1 (MQ=255)
                          cATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTGCTAGGGtt                                                                <  2:608065/91‑1 (MQ=25)
                              tACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttcttc                                                            >  1:1207429/1‑91 (MQ=255)
                               aCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCGACATAACTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTtcttct                                                           <  2:2164326/91‑1 (MQ=255)
                               aCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTGCTAGGGttctaat                                                           <  1:1592485/91‑4 (MQ=18)
                               aCTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttcttgt                                                           >  2:535930/1‑89 (MQ=255)
                                 tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACCTAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGCTCTTCttt                                                         >  2:1595390/1‑91 (MQ=255)
                                 tttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCttt                                                         <  1:3045180/91‑1 (MQ=255)
                                      aaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGaa                                                    >  2:3377937/1‑91 (MQ=255)
                                      aaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGaa                                                    <  1:2295551/91‑1 (MQ=255)
                                        tGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAAtt                                                  >  1:424568/1‑91 (MQ=255)
                                         gATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTa                                                 <  1:1626369/91‑1 (MQ=255)
                                          aTTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTAt                                                <  1:1647101/91‑1 (MQ=255)
                                          aTTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTAt                                                >  2:1850266/1‑91 (MQ=255)
                                              aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                            >  1:981798/1‑91 (MQ=255)
                                                aaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATc                                          >  1:255843/1‑91 (MQ=255)
                                                 aGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCg                                         <  2:3713624/91‑1 (MQ=255)
                                                 aGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCg                                         <  2:1952496/91‑1 (MQ=255)
                                                     aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                     >  2:3164933/1‑91 (MQ=255)
                                                     aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                     >  1:2150783/1‑91 (MQ=255)
                                                     aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                     >  1:2152406/1‑91 (MQ=255)
                                                      cATGTCGCATAACACTTCAACATAACTAGAACCCTAGCACAGATATCTATCTTTGTCAGGGACCTTCCTTATGAATTATGAGATCGaaccc                                    >  2:1489728/1‑88 (MQ=255)
                                                       aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                   >  1:1006819/1‑91 (MQ=255)
                                                       aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                   >  2:3690206/1‑91 (MQ=255)
                                                        tGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTc                                  <  1:2991315/91‑1 (MQ=255)
                                                            gCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGCTCGAATCTTCAACt                              >  1:371988/1‑91 (MQ=255)
                                                            gCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACt                              <  2:1638865/91‑1 (MQ=255)
                                                            gCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACt                              <  2:3815624/91‑1 (MQ=255)
                                                            gCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACt                              >  1:1670994/1‑91 (MQ=255)
                                                                aaCACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaaaa                          >  2:2127180/1‑91 (MQ=255)
                                                                 acacTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAAt                         <  1:858700/91‑1 (MQ=255)
                                                                  cacTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAAtt                        <  2:275805/91‑1 (MQ=255)
                                                                  cacTTCAACAGAAGTGGAACCCGAGCAAAGATAGCTATCTTGGTCAGGGTTCTTCTGTATGAATTATGAGATCGAATCTTCAACTAAAAtt                        <  2:1556645/91‑1 (MQ=255)
                                                                   acTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTa                       >  2:2970623/1‑91 (MQ=255)
                                                                    cTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTAc                      >  1:1321265/1‑91 (MQ=255)
                                                                      tCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGc                    <  2:3187123/91‑1 (MQ=255)
                                                                      tCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGc                    <  1:821915/91‑1 (MQ=255)
                                                                        aaCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCt                  <  2:1348184/91‑1 (MQ=255)
                                                                         aCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCtt                 >  2:882000/1‑91 (MQ=255)
                                                                         aCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCtt                 <  1:1098600/91‑1 (MQ=255)
                                                                           aTAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAGATTACGCCTTCa               >  2:3290953/1‑91 (MQ=255)
                                                                           aTAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAACATTACGCCTTCa               >  2:3290558/1‑91 (MQ=255)
                                                                             aaTTAGAACCCTAGCACAGATATCTCTCTTCGTCAGGGTTCTTCTTTATGAATCATGAGACCGAATCTCCAACTAAACTTACGCCTTCATa             >  2:2595401/1‑91 (MQ=255)
                                                                             aaTTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATa             >  1:1122530/1‑91 (MQ=255)
                                                                              aTTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAg            >  1:2629824/1‑91 (MQ=255)
                                                                              aTTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGATCTTCTTTATGAACTATGAGATCGAATGTTCAACTAAAATTACGCCTCCATAg            >  1:3673613/1‑91 (MQ=255)
                                                                               gaaGAACCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGc           >  1:2111366/3‑91 (MQ=18)
                                                                                tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa          <  2:3006748/91‑1 (MQ=255)
                                                                                tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa          >  1:1234304/1‑91 (MQ=255)
                                                                                 aGAACCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCaa         >  2:83731/1‑91 (MQ=21)
                                                                                  gAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg        >  1:3533859/1‑91 (MQ=255)
                                                                                  gAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg        <  1:3378388/91‑1 (MQ=255)
                                                                                  gAACCCTAGCAAAGATATCTATCTTGGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg        <  2:22124/91‑1 (MQ=255)
                                                                                   aaCCCTAGCAAAGATATCTATCTTTGTCGGGGTTCTTCTTTATGAATTATGAAATCCAATCTTCAACTAAAATTACGACTTCATATCAAgt       >  2:2306418/1‑90 (MQ=255)
                                                                                   aaCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAgg       <  1:1379123/91‑1 (MQ=255)
                                                                                    aCCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa      >  2:478249/1‑91 (MQ=25)
                                                                                    aCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa      <  2:340117/91‑1 (MQ=255)
                                                                                     cccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGAc     <  1:2294864/91‑1 (MQ=255)
                                                                                     cccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGAc     <  1:1180638/91‑1 (MQ=255)
                                                                                        tAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAtt  <  1:1020172/91‑1 (MQ=255)
                                                                                          |                                                                                        
GCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATT  >  USA300TCH1516_ALE/2559741‑2559919

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: