Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,559,830 A→G 11.1% intergenic (+24/+99) 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,8300AG11.1% 578.1 / 46.1 207intergenic (+24/+99)bcr_2/USA300HOU_RS13190Bicyclomycin resistance protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (104/78);  new base G (3/20);  total (108/99)
Fisher's exact test for biased strand distribution p-value = 7.97e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 208 total aligned reads displayed.

AGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATT  >  USA300TCH1516_ALE/2559740‑2559919
                                                                                          |                                                                                         
 gCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAg                                                                                          <  1:2805189/91‑1 (MQ=255)
  ccATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGc                                                                                         <  2:2684257/91‑1 (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)
       cTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGa                                                                                    >  2:1196259/1‑91 (MQ=255)
         aTTAGTCAGTCTACAAATCATTTACTTTAAAGTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGata                                                                                  <  2:1280764/91‑1 (MQ=255)
         aTTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGata                                                                                  >  1:612839/1‑91 (MQ=255)
         aTTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGata                                                                                  >  1:2426711/1‑91 (MQ=255)
          ttAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGatat                                                                                 <  2:689872/91‑1 (MQ=255)
           tAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCAAAACACTTCAACATAATTAGAACCCTAGCAAAGataga                                                                                <  2:3379135/91‑3 (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:2947017/1‑91 (MQ=255)
               cagtcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGAtatctatc                                                                            >  1:2341925/1‑91 (MQ=255)
                  tcTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCttt                                                                         <  1:2733959/91‑1 (MQ=38)
                   cTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTGAACATAATTAGAACCCTGACAAAGATAGATATCTTTg                                                                        <  1:3326470/91‑1 (MQ=37)
                   cTACAAATCATTTACTGTAAAATGATTAAAAAGCAACATATCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTg                                                                        <  2:346135/91‑1 (MQ=37)
                    tACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgc                                                                       <  2:3147748/91‑2 (MQ=37)
                    tACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgc                                                                       <  2:1908064/91‑2 (MQ=37)
                    tACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgc                                                                       <  1:3870316/91‑2 (MQ=37)
                    tACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgc                                                                       <  2:3125098/91‑2 (MQ=38)
                    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                                                                     >  2:1098600/1‑91 (MQ=255)
                      cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCa                                                                     >  2:2180429/1‑91 (MQ=255)
                      cAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCa                                                                     >  1:3790848/1‑91 (MQ=255)
                      cAAAGCAGTGACTGTAAGATGATTAAAAAGCAACATGTCGCATAACGCGTCAACATAATTAGAACCCTGACAAAGATAGATATCTTTgcta                                                                     <  2:1555126/91‑4 (MQ=16)
                       aaaTCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAg                                                                    >  1:3493715/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)
                              ttACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttctt                                                             >  1:3001101/1‑91 (MQ=255)
                               tACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGttcttc                                                            <  2:1563988/91‑1 (MQ=255)
                               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)
                                   ttAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTa                                                        >  1:3075421/1‑91 (MQ=255)
                                       aaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGaa                                                    >  2:3377937/1‑91 (MQ=255)
                                       aaTGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGaa                                                    <  1:2295551/91‑1 (MQ=255)
                                          gATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTa                                                 <  1:1626369/91‑1 (MQ=255)
                                               aaaaGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATgaga                                            >  1:981755/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)
                                                    cAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGaa                                       >  1:3213469/1‑91 (MQ=255)
                                                      aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                     >  1:2150783/1‑91 (MQ=255)
                                                      aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                     >  2:3164933/1‑91 (MQ=255)
                                                      aCATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATc                                     >  1:2152406/1‑91 (MQ=255)
                                                       cATGTCGCATAACACTTCAACATAATTAGAACCCCAGCAAAGATATCTATCCCCGCCAGCGTCCTCCCTAACGAATTACGAGATCGCATCt                                    >  2:1634530/1‑91 (MQ=255)
                                                       cATGTCGCATAACACTTCAACATAACTAGAACCCTAGCACAGATATCTATCTTTGTCAGGGACCTTCCTTATGAATTATGAGATCGaaccc                                    >  2:1489728/1‑88 (MQ=255)
                                                        aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                   >  1:1006819/1‑91 (MQ=255)
                                                        aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                   >  1:402255/1‑91 (MQ=255)
                                                        aTGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCtt                                   >  2:3690206/1‑91 (MQ=255)
                                                           tCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCaa                                <  1:777359/91‑1 (MQ=255)
                                                             gCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGCTCGAATCTTCAACt                              >  1:371988/1‑91 (MQ=255)
                                                             gCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACt                              <  2:3815624/91‑1 (MQ=255)
                                                             gCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCAAATCTTCAACt                              <  2:1742884/91‑1 (MQ=255)
                                                              cATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTa                             <  1:468501/91‑1 (MQ=255)
                                                                 aaCACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTaaaa                          >  2:2127180/1‑91 (MQ=255)
                                                                    acTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTa                       >  2:2970623/1‑91 (MQ=255)
                                                                     cTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTAc                      >  1:1321265/1‑91 (MQ=255)
                                                                      ttCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACg                     >  2:2625660/1‑91 (MQ=255)
                                                                       tCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGc                    <  1:821945/91‑1 (MQ=255)
                                                                       tCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGc                    <  2:3187123/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)
                                                                          aCATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCtt                 >  1:1652391/1‑91 (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)
                                                                                ttAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGACTCTTCAACTAAAATTACGCCTTCATAGc           >  2:2606308/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)
                                                                                 tAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCa          <  1:3240269/91‑1 (MQ=255)
                                                                                  aGAACCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCaa         >  2:83731/1‑91 (MQ=21)
                                                                                   gAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg        <  1:3378397/91‑1 (MQ=255)
                                                                                   gAACCCTAGCAAAGATATCTATCTTGGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAg        <  2:22124/91‑1 (MQ=255)
                                                                                    aaCCCTAGCAAAGATATCTATCTTTGTCGGGGTTCTTCTTTATGAATTATGAAATCCAATCTTCAACTAAAATTACGACTTCATATCAAgt       >  2:2306418/1‑90 (MQ=255)
                                                                                    aaCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAgg       >  2:1710924/1‑91 (MQ=255)
                                                                                    aaCCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAgg       <  1:1379123/91‑1 (MQ=255)
                                                                                     aCCCTGACAAAGATAGATATCTTTGCTAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGa      >  2:478249/1‑91 (MQ=25)
                                                                                      cccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGAc     <  1:2294864/91‑1 (MQ=255)
                                                                                      cccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGAc     >  1:556250/1‑91 (MQ=255)
                                                                                      cccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGAc     <  1:1180638/91‑1 (MQ=255)
                                                                                       ccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACa    >  1:1772121/1‑91 (MQ=255)
                                                                                       ccTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTAATAGCAAGGACa    <  1:3024613/91‑1 (MQ=255)
                                                                                         tAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACAtt  <  1:1020172/91‑1 (MQ=255)
                                                                                          |                                                                                         
AGCCATTCTATTAGTCAGTCTACAAATCATTTACTTTAAAATGATTAAAAAGCAACATGTCGCATAACACTTCAACATAATTAGAACCCTAGCAAAGATATCTATCTTTGTCAGGGTTCTTCTTTATGAATTATGAGATCGAATCTTCAACTAAAATTACGCCTTCATAGCAAGGACATT  >  USA300TCH1516_ALE/2559740‑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: