Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,172,075 G→C 11.0% E42Q (GAA→CAA)  leuB → 3‑isopropylmalate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,172,0750GC11.0% 629.5 / 53.3 193E42Q (GAA→CAA) leuB3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (91/79);  new base C (11/10);  total (103/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 195 total aligned reads displayed.

GGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTG  >  USA300TCH1516_ALE/2171985‑2172164
                                                                                           |                                                                                          
ggAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACg                                                                                            <  2:1669316/91‑1 (MQ=255)
  aaTCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGaa                                                                                          <  1:984550/91‑1 (MQ=255)
   aTCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGCGCACCACGAAt                                                                                         >  2:3842933/1‑91 (MQ=255)
   aTCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAAt                                                                                         >  2:2760197/1‑91 (MQ=255)
   aTCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAAt                                                                                         >  2:1140666/1‑91 (MQ=255)
     cGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAAttt                                                                                       >  2:993471/1‑91 (MQ=255)
       gTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAAtttcg                                                                                     >  2:3952023/1‑89 (MQ=255)
       gTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtgg                                                                                     >  1:3815809/1‑91 (MQ=255)
       gTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtgg                                                                                     <  1:2439174/91‑1 (MQ=255)
        tCCAGAAATTTTGAACGGATCGCGATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtggt                                                                                    <  2:306170/91‑1 (MQ=255)
         ccAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtggtg                                                                                   >  2:2949171/1‑91 (MQ=255)
          cAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtggtgg                                                                                  <  1:1022387/91‑1 (MQ=255)
          cAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtggtgg                                                                                  <  2:2016868/91‑1 (MQ=255)
          cAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtggtgg                                                                                  <  1:1506857/91‑1 (MQ=255)
           aGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTggtggt                                                                                 >  1:1708944/1‑91 (MQ=255)
            gAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTcgtggt                                                                                 <  2:3825896/91‑1 (MQ=38)
            gAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGgtggtg                                                                                <  2:905593/91‑1 (MQ=255)
            gAAATTTTGAACGGATCTCTATCATTGCCTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTcgtggt                                                                                 <  2:510795/91‑1 (MQ=37)
              aaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTCGTGGTGc                                                                               <  1:175118/91‑1 (MQ=37)
              aaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGcc                                                                              <  2:3279592/91‑1 (MQ=255)
              aaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGcc                                                                              <  2:2253535/91‑1 (MQ=255)
              aaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGcc                                                                              >  2:193617/1‑91 (MQ=255)
               aTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCct                                                                             <  2:1474128/91‑1 (MQ=255)
                ttttGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCctc                                                                            >  1:3378616/1‑91 (MQ=255)
                   tGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTAt                                                                        >  2:3180745/1‑91 (MQ=255)
                    gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTAtt                                                                       >  1:2581322/1‑91 (MQ=255)
                     aaCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTg                                                                      >  2:2197157/1‑91 (MQ=255)
                     aaCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTg                                                                      >  1:3508861/1‑91 (MQ=255)
                      aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGa                                                                     >  2:2717917/1‑91 (MQ=255)
                      aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGa                                                                     <  1:728662/91‑1 (MQ=255)
                      aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGa                                                                     >  1:1444419/1‑91 (MQ=255)
                      aCGGATCTCTATCACTGCTTGAAGCTATAAGTAATAAATATAAGTTTAATTATCGAATAGAGCACCACCAAATTCGTGGTGCTCTAtttga                                                                    <  2:3677822/91‑4 (MQ=16)
                         gATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATAc                                                                  >  1:125484/1‑91 (MQ=255)
                           tctcTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACAt                                                                <  1:996770/91‑1 (MQ=255)
                           tctcTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACAt                                                                <  1:3598111/91‑1 (MQ=255)
                           tctcTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACAt                                                                >  2:409652/1‑91 (MQ=255)
                              ctATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCg                                                             <  1:801916/91‑1 (MQ=255)
                               tATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCgg                                                            >  2:2928689/1‑91 (MQ=255)
                                aTCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGc                                                           >  2:389062/1‑91 (MQ=255)
                                aTCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGc                                                           >  2:2319687/1‑91 (MQ=255)
                                  cATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGa                                                         >  1:329090/1‑91 (MQ=255)
                                   aTTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAg                                                        >  2:3387012/1‑91 (MQ=255)
                                    ttGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGc                                                       >  2:2309772/1‑91 (MQ=255)
                                        ttGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCttt                                                   >  1:463257/1‑91 (MQ=255)
                                         tGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTa                                                  >  1:2305761/1‑91 (MQ=255)
                                          gAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTaa                                                 <  1:461694/91‑1 (MQ=255)
                                           aaaTTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAAc                                                >  2:1749700/1‑91 (MQ=255)
                                           aaaTTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAAc                                                <  1:925849/91‑1 (MQ=255)
                                            aaTTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACt                                               <  1:2929064/91‑1 (MQ=255)
                                            aaTTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACt                                               <  2:862877/91‑1 (MQ=255)
                                              ttATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTga                                             <  1:1657558/91‑1 (MQ=255)
                                                  aaGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGaaa                                         >  1:2479627/1‑91 (MQ=255)
                                                  aaGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGaaa                                         >  1:2160224/1‑91 (MQ=255)
                                                   aGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGcaaa                                        >  2:3848706/1‑91 (MQ=255)
                                                   aGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGaaaa                                        >  1:2492615/1‑91 (MQ=255)
                                                   aGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGaaaa                                        >  1:895480/1‑91 (MQ=255)
                                                     taataaATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAAcc                                      >  1:2005855/1‑91 (MQ=255)
                                                      aataaataTAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCt                                     <  1:2216033/91‑1 (MQ=255)
                                                         aaataTAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTaa                                  <  1:1060265/91‑1 (MQ=255)
                                                            tataACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATg                               <  2:1359328/91‑1 (MQ=255)
                                                              taACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCg                             <  2:3039409/91‑1 (MQ=255)
                                                                 cTTTAATTATTAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCgtgt                          >  2:1055697/1‑91 (MQ=255)
                                                                 cTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCgtgt                          >  1:1276026/1‑91 (MQ=255)
                                                                 cTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCCCGCCCTTAACTGAGAAAACCTTAACTGCgtgt                          >  2:1611520/1‑91 (MQ=255)
                                                                  tttAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTa                         >  1:137995/1‑91 (MQ=255)
                                                                   ttAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaa                        >  1:1330236/1‑91 (MQ=255)
                                                                    tAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaaa                       >  1:3923710/1‑91 (MQ=255)
                                                                     aaTTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTaaaa                      >  1:1145181/1‑91 (MQ=255)
                                                                      aTTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAag                     <  1:362775/91‑1 (MQ=255)
                                                                      aTTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGGAAAag                     <  2:2434247/91‑1 (MQ=255)
                                                                          atcaaTAGAGGCACCACCAAATTCGTGGTGCTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGc                  >  1:3751531/4‑91 (MQ=17)
                                                                           tcaaTAGAGGCACCACCAAATTCGTGGTGCTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCa                 >  1:711013/3‑91 (MQ=17)
                                                                           cAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGATCAg                >  1:3194616/1‑91 (MQ=255)
                                                                           cAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAg                >  1:49484/1‑91 (MQ=255)
                                                                           cAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAg                >  1:871095/1‑91 (MQ=255)
                                                                           cAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAg                >  2:1114447/1‑91 (MQ=255)
                                                                            caaTAGAGGCACCACCAAATTCGAGGTGCTCTATTTGATACATCCGGCGAGCCTTTAACTGAGAAACCCTTAAATGCCTGTAAAAGAGCAg                >  2:2145209/2‑91 (MQ=9)
                                                                            aaaTAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGa               >  1:2806684/1‑91 (MQ=255)
                                                                             aaTAGAGGCACCACCAAATTCGTGGTGCTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGa               >  2:329236/1‑91 (MQ=18)
                                                                              aTAGAGGCACCACCAAATTCGTGGTGCTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGAt              >  2:3107851/1‑91 (MQ=18)
                                                                              aTAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAATAGCatctg             <  1:571683/91‑5 (MQ=255)
                                                                              aTAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATg             >  1:3905014/1‑91 (MQ=255)
                                                                               tAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGc            <  2:3771425/91‑1 (MQ=255)
                                                                                agagCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAACCCCTTAAATGCGTGTAAAATAGCATCTGCt           <  2:2651595/91‑1 (MQ=255)
                                                                                 gaggCACCACCAAATTCGTGGTGCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTa          >  1:2981733/4‑91 (MQ=37)
                                                                                 gagCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCATTTAACTGAGAAAGCCTTAAATGCGTGTAAAAGAGCAGATGCTa          <  2:3873192/91‑1 (MQ=255)
                                                                                  agCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAt         <  1:2338974/91‑1 (MQ=255)
                                                                                   ggCACCACCAAATTCGTGGTGCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtt        >  1:1842209/2‑91 (MQ=37)
                                                                                    gCACCACGAATTTGGTGGTGCCTCTATTGCTACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtt        >  2:3446694/1‑91 (MQ=255)
                                                                                    gCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtt        <  2:626283/91‑1 (MQ=255)
                                                                                    gCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtt        <  1:3553625/91‑1 (MQ=255)
                                                                                    gCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtt        <  1:3542359/91‑1 (MQ=255)
                                                                                     caccacGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAttt       <  2:2373912/91‑1 (MQ=255)
                                                                                     caccacGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAttt       <  1:2568621/91‑1 (MQ=255)
                                                                                      accacGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGACAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtttt      >  2:3118295/1‑91 (MQ=255)
                                                                                      accacGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtttt      >  2:2972482/1‑91 (MQ=255)
                                                                                       ccacGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTa     >  1:2289436/1‑91 (MQ=255)
                                                                                       ccacGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTa     >  2:1097410/1‑91 (MQ=255)
                                                                                       ccacGAATTGGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTa     >  1:1977069/1‑91 (MQ=255)
                                                                                          cGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTg  >  2:688512/1‑91 (MQ=255)
                                                                                           |                                                                                          
GGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTG  >  USA300TCH1516_ALE/2171985‑2172164

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: