Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 79,612 A→T 11.5% intergenic (+221/‑39) USA300HOU_RS00355 → / → cmoB hypothetical protein/tRNA U34 carboxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE79,6120AT11.5% 653.4 / 59.9 218intergenic (+221/‑39)USA300HOU_RS00355/cmoBhypothetical protein/tRNA U34 carboxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (89/104);  new base T (9/16);  total (98/120)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

Only 100 of 231 total aligned reads displayed.

TCTCGCGACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATG  >  USA300TCH1516_ALE/79523‑79702
                                                                                         |                                                                                          
tctcGCGACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAAc                                                                                           >  1:807899/1‑91 (MQ=255)
tctcGCGACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAAc                                                                                           <  1:2130959/91‑1 (MQ=255)
tctcGCGACCACAGCTTAATTTTTGGTGTTCAATCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAAc                                                                                           <  1:3296168/91‑1 (MQ=255)
 ctcGCGACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACg                                                                                          >  2:3018427/1‑91 (MQ=255)
  tcGCGACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGt                                                                                         >  1:2654165/1‑91 (MQ=255)
  tcGCGACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGt                                                                                         <  2:3873138/91‑1 (MQ=255)
    gcgACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTaa                                                                                       >  1:2835617/1‑91 (MQ=255)
    gcgACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTaa                                                                                       >  2:1141191/1‑91 (MQ=255)
    gcgACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTaa                                                                                       >  1:912307/1‑91 (MQ=255)
       aCCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAatta                                                                                    <  2:2113133/91‑1 (MQ=255)
        ccACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAataat                                                                                   <  2:3695122/91‑4 (MQ=255)
        ccACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAattat                                                                                   <  1:1388705/91‑1 (MQ=255)
         cacaGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAattatt                                                                                  >  2:1074913/1‑91 (MQ=255)
         cacaGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAattatt                                                                                  <  2:2273891/91‑1 (MQ=255)
          acaGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAAtaatta                                                                                 <  2:574487/91‑1 (MQ=255)
          acaGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAAttatta                                                                                 >  1:3472907/1‑91 (MQ=255)
             gCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACgt                                                                              <  2:3530264/91‑2 (MQ=255)
             gCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTACTTATTACGa                                                                              >  1:3719761/1‑91 (MQ=255)
             gCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGa                                                                              >  2:3488554/1‑91 (MQ=255)
             gCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGa                                                                              <  1:3214195/91‑1 (MQ=255)
              cTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGAt                                                                             <  1:1919229/91‑1 (MQ=255)
               ttAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACgttt                                                                            <  1:3776811/91‑4 (MQ=39)
               ttAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGAtt                                                                            <  1:2392677/91‑1 (MQ=255)
               ttAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTAATATTCAAAACGTAATTATTACGAtt                                                                            >  1:3615031/1‑91 (MQ=255)
                 aaTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGTTTTg                                                                          >  2:2524918/1‑91 (MQ=39)
                 aaTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTg                                                                          >  2:2961614/1‑91 (MQ=255)
                  aTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGTTTTGa                                                                         <  2:3437537/91‑1 (MQ=38)
                  aTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGa                                                                         >  1:654140/1‑91 (MQ=255)
                   tttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGTTTTGaa                                                                        <  1:835402/91‑1 (MQ=38)
                    ttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGTTTTGAAt                                                                       <  1:2524918/91‑1 (MQ=38)
                    ttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGTTTTGAAt                                                                       <  1:3248763/91‑1 (MQ=38)
                    ttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGTTTTGAAt                                                                       <  1:3247334/91‑1 (MQ=38)
                    ttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAAt                                                                       <  2:2948560/91‑1 (MQ=255)
                     tttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAAtt                                                                      >  2:3663011/1‑91 (MQ=255)
                     tttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAAtt                                                                      <  2:2948151/91‑1 (MQ=255)
                       tGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTaa                                                                    >  1:2313726/1‑91 (MQ=255)
                       tGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTaa                                                                    >  1:947588/1‑91 (MQ=255)
                          tgtTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCg                                                                 >  2:1265984/1‑91 (MQ=255)
                            ttcattcaTGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgag                                                               >  1:3318900/1‑91 (MQ=255)
                            ttcattcaTGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgag                                                               <  2:362270/91‑1 (MQ=255)
                             tcattcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATAAAGCgagg                                                              <  2:507449/91‑1 (MQ=255)
                              cattcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgagga                                                             <  2:1399690/91‑1 (MQ=255)
                                ttcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGgaga                                                           <  2:1920644/91‑1 (MQ=255)
                                 tcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaa                                                          <  2:2707728/91‑1 (MQ=255)
                                 tcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaa                                                          <  2:633190/91‑1 (MQ=255)
                                  catGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaat                                                         <  1:1427859/91‑1 (MQ=255)
                                     gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                      >  2:2662574/1‑91 (MQ=255)
                                     gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                      >  2:1129561/1‑91 (MQ=255)
                                     gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                      >  1:1664991/1‑91 (MQ=255)
                                      aCTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaa                                                     >  1:2882064/1‑91 (MQ=255)
                                      aCTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaa                                                     <  2:639371/91‑1 (MQ=255)
                                       cTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaat                                                    <  2:784949/91‑1 (MQ=255)
                                         ttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgcaatga                                                  >  2:3568775/1‑91 (MQ=255)
                                         ttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatga                                                  <  1:1666638/91‑1 (MQ=255)
                                           aaaaaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACt                                                <  2:614552/91‑1 (MQ=255)
                                           aaaaaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACt                                                <  1:1074913/91‑1 (MQ=255)
                                           aaaaaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACt                                                <  1:3140758/91‑1 (MQ=255)
                                           aaaaaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACt                                                <  1:1722420/91‑1 (MQ=255)
                                            aaaaaTCCATATTGCATAAATGTACGTAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTa                                               <  2:3366838/91‑1 (MQ=255)
                                             aaaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTaa                                              <  1:1005179/91‑1 (MQ=255)
                                              aaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAg                                             >  2:160838/1‑91 (MQ=255)
                                              aaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAg                                             >  2:416274/1‑91 (MQ=255)
                                              aaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAg                                             >  2:3276113/1‑91 (MQ=255)
                                               aaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGa                                            <  1:2846924/91‑1 (MQ=255)
                                                 tCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACt                                          >  2:145744/1‑91 (MQ=255)
                                                  ccTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACtt                                         >  1:1896137/1‑91 (MQ=255)
                                                   cTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTa                                        <  1:797259/91‑1 (MQ=255)
                                                   cTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTa                                        <  1:2212933/91‑1 (MQ=255)
                                                     tATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATg                                      >  1:1778817/1‑91 (MQ=255)
                                                     tATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATg                                      <  1:1325923/91‑1 (MQ=255)
                                                     tATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATg                                      <  1:3705208/91‑1 (MQ=255)
                                                      aTTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGa                                     >  2:183249/1‑91 (MQ=255)
                                                      aTTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGa                                     <  1:614087/91‑1 (MQ=255)
                                                       ttGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGAc                                    >  1:2440577/1‑91 (MQ=255)
                                                        tGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACg                                   <  2:766525/91‑1 (MQ=255)
                                                             aaaTGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTtgg                              <  2:692189/91‑1 (MQ=255)
                                                             aaaTGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTtgg                              >  1:2266479/1‑91 (MQ=255)
                                                             aaaTGTACATAGTGTAGGACTATTCAAAACTTAATTATTACGATTTGAATTAAGCGAGGAGCATGAAATGACTAAGACTTATGACGTTtgg                              <  1:2080380/91‑1 (MQ=255)
                                                              aaTGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTtggt                             <  1:2580410/91‑1 (MQ=255)
                                                                tGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTtggtgg                           <  2:477806/91‑1 (MQ=255)
                                                                   aCATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCaa                        >  1:1994523/1‑91 (MQ=255)
                                                                        gtgtAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAgg                   >  1:3088060/1‑91 (MQ=255)
                                                                         tgtAGTAATATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGt                  <  2:3457478/91‑1 (MQ=255)
                                                                            agtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCaa               >  1:1871401/1‑91 (MQ=255)
                                                                            agtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCaa               <  2:3088060/91‑1 (MQ=255)
                                                                              taCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGa             <  2:2544380/91‑1 (MQ=255)
                                                                              taCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGa             <  2:1548970/91‑1 (MQ=255)
                                                                                cTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAAt           >  1:2867653/1‑91 (MQ=255)
                                                                                cTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAAt           <  2:3791157/91‑1 (MQ=255)
                                                                                  aTTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCa         <  1:1177998/91‑1 (MQ=255)
                                                                                    tCAAATCGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAga       >  2:3130023/1‑91 (MQ=38)
                                                                                    tCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAga       <  2:1651952/91‑1 (MQ=255)
                                                                                    tCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAga       <  2:3159345/91‑1 (MQ=255)
                                                                                     cAAATCGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgat      >  2:3284788/1‑91 (MQ=38)
                                                                                     cAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgat      >  2:1329844/1‑91 (MQ=255)
                                                                                      aaaaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatg     >  1:1517744/1‑91 (MQ=255)
                                                                                      aaaaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatg     >  2:3616427/1‑91 (MQ=255)
                                                                                      aaaaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatg     <  2:2800319/91‑1 (MQ=255)
                                                                                        aaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgat   <  2:3579350/91‑1 (MQ=255)
                                                                                         aCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgatg  >  1:1777151/1‑91 (MQ=255)
                                                                                         |                                                                                          
TCTCGCGACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATG  >  USA300TCH1516_ALE/79523‑79702

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: