Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,231,474 G→C 11.0% G40A (GGT→GCT)  stp → Serine/threonine phosphatase stp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,231,4740GC11.0% 317.3 / 24.6 102G40A (GGT→GCT) stpSerine/threonine phosphatase stp
Reads supporting (aligned to +/- strand):  ref base G (50/37);  new base C (8/3);  total (60/42)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

CTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCA  >  USA300TCH1516_ALE/1231384‑1231564
                                                                                          |                                                                                          
cTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATggg                                                                                            <  1:2434914/91‑1 (MQ=255)
cTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATggg                                                                                            <  1:1865448/91‑1 (MQ=255)
 tGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGt                                                                                           <  1:3688078/91‑1 (MQ=255)
  ggACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTg                                                                                          >  1:1045614/1‑91 (MQ=255)
  ggACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTg                                                                                          >  1:306401/1‑91 (MQ=255)
  ggACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTg                                                                                          >  1:548802/1‑91 (MQ=255)
  ggACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGCTATGGGTg                                                                                          <  1:4078012/91‑1 (MQ=255)
    acaacaTAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCCAACTAATCAACAACTTTTAGTTCTGTGTGATGGTAtggctggc                                                                                        >  1:4148267/1‑91 (MQ=255)
     caaccTAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGcc                                                                                       <  1:3455698/91‑1 (MQ=255)
     caacaTAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGcc                                                                                       >  1:1843386/1‑91 (MQ=255)
        caTAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATa                                                                                    >  1:2564220/1‑91 (MQ=255)
          tAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATaaa                                                                                  <  1:789248/91‑1 (MQ=255)
            gagaTAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGc                                                                                >  1:3493448/1‑91 (MQ=255)
              gaTAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGCTGGCCATAAAGCAg                                                                              >  1:2706364/1‑91 (MQ=255)
               aTAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGCTGGCCATAAAGCAgg                                                                             >  1:1920915/1‑91 (MQ=255)
                 aaGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGgag                                                                           <  1:4002736/91‑1 (MQ=255)
                 aaGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGgag                                                                           <  1:2010213/91‑1 (MQ=255)
                    aaTGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAg                                                                        >  1:611246/1‑91 (MQ=255)
                     aTGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGt                                                                       >  1:1027934/1‑91 (MQ=255)
                       gAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGCTGGCCATAAAGCAGGAGAAGTTg                                                                     <  1:3093901/91‑1 (MQ=255)
                         aGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCa                                                                   >  1:1052154/1‑91 (MQ=255)
                          gATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCaa                                                                  >  1:725329/1‑91 (MQ=255)
                          gATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCaa                                                                  >  1:726225/1‑91 (MQ=255)
                           aTGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAg                                                                 <  1:2870780/91‑1 (MQ=255)
                           aTGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGCTGGCCATAAAGCAGGAGAAGTTGCAAg                                                                 <  1:1818606/91‑1 (MQ=255)
                             gCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTa                                                               <  1:2539189/91‑1 (MQ=255)
                                 gtggtATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGCTGGCCATAAAGCAGGAGAAGTTGCAAGTAAAtt                                                           >  1:2310836/1‑91 (MQ=255)
                                  tggtATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAAttt                                                          <  1:3881159/91‑1 (MQ=255)
                                      aTTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTa                                                      >  1:289956/1‑91 (MQ=255)
                                      aTTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTa                                                      >  1:1738994/1‑91 (MQ=255)
                                       ttttttATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTAc                                                     >  1:1409578/1‑91 (MQ=255)
                                        tttttATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACa                                                    <  1:2237588/91‑1 (MQ=255)
                                         ttttATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAg                                                   >  1:382775/1‑91 (MQ=255)
                                             ataATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGa                                               <  1:202836/91‑1 (MQ=255)
                                                aTCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGTCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGtt                                            >  1:3591669/1‑91 (MQ=255)
                                                aTCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGtt                                            >  1:3479991/1‑91 (MQ=255)
                                                aTCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGtt                                            >  1:3480186/1‑91 (MQ=255)
                                                aTCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGtt                                            >  1:3480328/1‑91 (MQ=255)
                                                aTCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGtt                                            >  1:271546/1‑91 (MQ=255)
                                                aTCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGATGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGtt                                            >  1:129877/1‑91 (MQ=255)
                                                 tCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTg                                           <  1:607498/91‑1 (MQ=255)
                                                 tCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTg                                           <  1:2662236/91‑1 (MQ=255)
                                                  cAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGa                                          >  1:2768281/1‑91 (MQ=255)
                                                  cAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGa                                          >  1:3816326/1‑91 (MQ=255)
                                                  cAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGa                                          <  1:2377362/91‑1 (MQ=255)
                                                  cAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGa                                          <  1:1223613/91‑1 (MQ=255)
                                                  cAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGa                                          >  1:1945769/1‑91 (MQ=255)
                                                  cAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGCTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGa                                          >  1:4119345/1‑91 (MQ=255)
                                                   aaaCTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGTCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGaa                                         <  1:631328/91‑1 (MQ=255)
                                                     aCTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAAt                                       <  1:3995654/91‑1 (MQ=255)
                                                     aCTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCAGAAAGCAGGCGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAAt                                       <  1:3311643/91‑1 (MQ=255)
                                                         aTCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGTCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCg                                   <  1:401307/91‑1 (MQ=255)
                                                         aTCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCg                                   >  1:2404198/1‑91 (MQ=255)
                                                         aTCAACAACTTTTAGTTCTGTGTGATGGTATGGATGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCg                                   <  1:2196517/91‑1 (MQ=255)
                                                          tCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGt                                  >  1:2029576/1‑91 (MQ=255)
                                                          tCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGt                                  >  1:1933184/1‑91 (MQ=255)
                                                          tCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGt                                  >  1:390169/1‑91 (MQ=255)
                                                          tCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGt                                  >  1:3046301/1‑91 (MQ=255)
                                                          tCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGt                                  <  1:2400899/91‑1 (MQ=255)
                                                          tCAACAACTTTTAGTTCTGTGTGATGGTATGGCTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGt                                  >  1:2998319/1‑91 (MQ=255)
                                                            aacaacTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGttt                                >  1:197086/1‑91 (MQ=255)
                                                            aacaacTTTTAGTTCTGTGTGATGGTATGGGTGGCCAAAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGttt                                >  1:3483182/1‑91 (MQ=255)
                                                              caacTTTTAGTTCTGTGTGATGGTATGGGTGTCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTg                              >  1:70775/1‑91 (MQ=255)
                                                              caacTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTg                              >  1:2689011/1‑91 (MQ=255)
                                                                acTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGaa                            >  1:213210/1‑91 (MQ=255)
                                                                acTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGaa                            >  1:55439/1‑91 (MQ=255)
                                                                 cTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAg                           <  1:1785243/91‑1 (MQ=255)
                                                                 cTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAg                           <  1:3783287/91‑1 (MQ=255)
                                                                  ttttAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGc                          <  1:2975076/91‑1 (MQ=255)
                                                                  ttttAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGc                          >  1:2503801/1‑91 (MQ=255)
                                                                  ttttAGTTCTGTGTGATGGTATGGATGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGc                          <  1:454747/91‑1 (MQ=255)
                                                                   tttAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCg                         >  1:708285/1‑91 (MQ=255)
                                                                    ttAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCgg                        <  1:898544/91‑1 (MQ=255)
                                                                    ttAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCgg                        <  1:1053356/91‑1 (MQ=255)
                                                                    ttAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCgg                        <  1:47500/91‑1 (MQ=255)
                                                                      aGTTCTGTGTGATGGTATGTGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGaa                      <  1:3150164/91‑1 (MQ=255)
                                                                      aGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGaa                      >  1:1066067/1‑91 (MQ=255)
                                                                      aGTTCTGTGTGATGCTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGaa                      >  1:791768/1‑91 (MQ=255)
                                                                       gTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGaaa                     >  1:2380911/1‑91 (MQ=255)
                                                                        ttCTGTGTGATGGTATGGCTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGaaaa                    >  1:3586194/1‑91 (MQ=255)
                                                                            gtgtgATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCtt                <  1:1930327/91‑1 (MQ=255)
                                                                            gtgtgATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCtt                >  1:449039/1‑91 (MQ=255)
                                                                            gtgtgATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCtt                >  1:2098146/1‑91 (MQ=255)
                                                                            gtgtgATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCtt                <  1:1930582/91‑1 (MQ=255)
                                                                             tgtgATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTta               <  1:3694526/91‑1 (MQ=255)
                                                                             tgtgATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTta               >  1:3838409/1‑91 (MQ=255)
                                                                              gtgATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTtat              <  1:522710/91‑1 (MQ=255)
                                                                               tgATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTtata             <  1:1567432/91‑1 (MQ=255)
                                                                               tgATGGTATGGATGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTtata             >  1:1866965/1‑91 (MQ=255)
                                                                                  tGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGaa          <  1:2679935/91‑1 (MQ=255)
                                                                                   ggTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGaac         <  1:1473144/91‑1 (MQ=255)
                                                                                   ggTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGaac         <  1:3358856/91‑1 (MQ=255)
                                                                                    taccgggTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACATATGAGTTGAAATCCCGTTTTGAAGCGGAATATCTTATAGaaca        >  1:4147006/5‑91 (MQ=255)
                                                                                    taccgggTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAGATCTTATAGaaca        >  1:4146935/5‑91 (MQ=255)
                                                                                    gTATGGCTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGaaca        >  1:1529946/1‑91 (MQ=255)
                                                                                    gTATGGCTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGaaca        >  1:349241/1‑91 (MQ=255)
                                                                                      aTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGaacaac      >  1:606474/1‑91 (MQ=255)
                                                                                      aTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGaacaac      >  1:2018102/1‑91 (MQ=255)
                                                                                        gggTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACAt    >  1:189960/1‑91 (MQ=255)
                                                                                          gTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCa  >  1:1536304/1‑91 (MQ=255)
                                                                                          |                                                                                          
CTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCA  >  USA300TCH1516_ALE/1231384‑1231564

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