Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 217,046 C→T 8.4% intergenic (‑240/+13) rocD2_1 ← / ← brnQ_1 Ornithine aminotransferase 2/Branched‑chain amino acid transport system 2 carrier protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE217,0460CT 8.4% 643.2 / 29.1 206intergenic (‑240/+13)rocD2_1/brnQ_1Ornithine aminotransferase 2/Branched‑chain amino acid transport system 2 carrier protein
Reads supporting (aligned to +/- strand):  ref base C (101/87);  new base T (16/2);  total (117/89)
Fisher's exact test for biased strand distribution p-value = 4.87e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 208 total aligned reads displayed.

TTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGAACAAC  >  USA300TCH1516_ALE/216957‑217136
                                                                                         |                                                                                          
 ttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACtg                                                                                          <  1:2708056/91‑1 (MQ=255)
 ttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACtg                                                                                          <  1:1230138/91‑1 (MQ=255)
  tgtgGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACtgt                                                                                         >  1:3395094/1‑91 (MQ=255)
  tgtgGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACtgt                                                                                         >  2:1946955/1‑91 (MQ=255)
  tgtgGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACtgt                                                                                         <  1:788029/91‑1 (MQ=255)
   gtgGATAATCATTGATCTGTATGTATCAGATTGAAGTTTACTGTAGTATGAATAAAAAATGCCACAAGCATATGTTTCACATATGCCtgtg                                                                                        >  1:2523509/1‑91 (MQ=255)
   gtgGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACtgtg                                                                                        >  2:2547023/1‑91 (MQ=255)
    tgGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTgg                                                                                       <  2:587805/91‑1 (MQ=255)
        tAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATa                                                                                   <  1:3256728/91‑1 (MQ=255)
          aTCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGt                                                                                 <  2:2293918/91‑1 (MQ=255)
            cATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGtat                                                                               <  2:3865179/91‑1 (MQ=255)
             aTTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGtatt                                                                              >  1:652288/1‑91 (MQ=255)
              ttGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGtatta                                                                             >  2:2093685/1‑91 (MQ=255)
               tGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGtattat                                                                            >  2:626141/1‑91 (MQ=255)
               tGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGtattat                                                                            >  2:1037224/1‑91 (MQ=255)
                  tatgtatgTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCt                                                                         >  2:1488437/1‑91 (MQ=255)
                    tgtatgtaTACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTtg                                                                       >  2:3877486/1‑91 (MQ=255)
                     gtatgtatACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTtgt                                                                      >  2:3043236/1‑91 (MQ=255)
                     gtatgtatACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTtgt                                                                      <  1:1676148/91‑1 (MQ=255)
                       atgtatACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTgg                                                                    <  2:2457226/91‑1 (MQ=255)
                       atgtatACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTgg                                                                    <  1:3582479/91‑1 (MQ=255)
                        tgtatACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGa                                                                   >  2:1838003/1‑91 (MQ=255)
                         gtatACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTggca                                                                  >  2:150760/1‑89 (MQ=255)
                          tataCGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAAt                                                                 >  2:2151478/1‑91 (MQ=255)
                           ataCGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAAtg                                                                <  2:3536005/91‑2 (MQ=255)
                              cGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAGCATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTa                                                             >  1:3601499/1‑91 (MQ=255)
                              cGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTa                                                             >  1:3793168/1‑91 (MQ=255)
                              cGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTa                                                             >  2:1672105/1‑91 (MQ=255)
                              cGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCCTGTGGAATTTTa                                                             >  2:1481392/1‑91 (MQ=255)
                                aTTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACg                                                           <  2:3169060/91‑1 (MQ=255)
                                aTTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACg                                                           >  1:3741534/1‑91 (MQ=255)
                                   gCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATg                                                        <  1:1044264/91‑1 (MQ=255)
                                   gCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATg                                                        <  1:3581008/91‑1 (MQ=255)
                                   gCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATg                                                        <  2:915410/91‑1 (MQ=255)
                                   gCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATg                                                        >  1:2742949/1‑91 (MQ=255)
                                     aGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCt                                                      >  1:2961929/1‑91 (MQ=255)
                                     aGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCt                                                      <  1:542613/91‑1 (MQ=255)
                                     aGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCt                                                      >  2:1875654/1‑91 (MQ=255)
                                      gTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCtt                                                     <  1:2491458/91‑1 (MQ=255)
                                       tttGATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTa                                                    <  1:3060651/91‑1 (MQ=255)
                                        ttAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTaa                                                   >  1:918183/1‑91 (MQ=255)
                                          aaTGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAAtt                                                 >  2:3429562/1‑91 (MQ=255)
                                            tGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAAtttt                                               <  1:1653409/91‑1 (MQ=255)
                                             gtagtaTGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTta                                              <  2:1637177/91‑1 (MQ=255)
                                             gtagtaTGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTta                                              >  2:49372/1‑91 (MQ=255)
                                             gtagtaTGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGACGCTTAATTTta                                              >  1:3512596/1‑91 (MQ=255)
                                                gtaTGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTataa                                           >  2:2863129/1‑91 (MQ=255)
                                                    gAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTataaataa                                       <  1:2547023/91‑1 (MQ=255)
                                                     aaTAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAcataat                                      >  2:3421820/1‑91 (MQ=255)
                                                     aaTAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAaataat                                      <  2:321963/91‑1 (MQ=255)
                                                      aTAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATg                                     <  1:2863701/91‑1 (MQ=255)
                                                        aaaaaaTGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGaa                                   <  2:377731/91‑1 (MQ=255)
                                                          aataTGCAACACGCATATGTTTATCTGATCACTGTGCCATATTATTATGCATGTGGAAGTTTACGATGCTTAATTTTATGAATAATGAAGc                                 <  2:586407/88‑1 (MQ=255)
                                                           aaaTGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGcc                                >  1:15117/1‑91 (MQ=255)
                                                            aaTGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCg                               <  1:2639915/91‑1 (MQ=255)
                                                            aaTGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTGATAAATAATGAAGCCg                               >  2:2982688/1‑91 (MQ=255)
                                                              tGCCACAGTCATATGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGAt                             >  2:1384492/1‑91 (MQ=25)
                                                              tGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGAt                             >  2:1681355/1‑91 (MQ=255)
                                                              tGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGAt                             >  2:3595361/1‑91 (MQ=255)
                                                               gCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATa                            >  1:2957554/1‑91 (MQ=255)
                                                               gCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATa                            >  1:1206145/1‑91 (MQ=255)
                                                               gCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATa                            >  1:3854400/1‑91 (MQ=255)
                                                               gCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATa                            >  2:2526989/1‑91 (MQ=255)
                                                                  acaAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAg                                  <  1:2332132/82‑1 (MQ=255)
                                                                  acaAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAg                                  >  2:2332132/1‑82 (MQ=255)
                                                                  acaAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATg                         >  1:1457849/1‑91 (MQ=255)
                                                                  acaAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATg                         <  2:3702993/91‑1 (MQ=255)
                                                                   caAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGa                        >  1:3165402/1‑91 (MQ=255)
                                                                   caAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGa                        >  1:1859987/1‑91 (MQ=255)
                                                                    agtcATATGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGaa                       >  1:382973/4‑91 (MQ=38)
                                                                    aaGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGaa                       >  2:3287055/1‑91 (MQ=255)
                                                                     gtcATATGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGaaa                      >  1:1141397/3‑91 (MQ=38)
                                                                     gtcATATGTTAAACATATGCTTGTGGCATAGCATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGaaa                      >  1:119353/3‑91 (MQ=37)
                                                                     aGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGaaa                      <  2:479901/91‑1 (MQ=255)
                                                                     aGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGaaa                      <  2:479994/91‑1 (MQ=255)
                                                                     aGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGaaa                      >  2:683063/1‑91 (MQ=255)
                                                                      tcATATGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAAc                     >  1:2749268/2‑91 (MQ=38)
                                                                       cATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAAcc                    <  1:2675466/91‑1 (MQ=255)
                                                                        atatGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCa                   >  2:2243993/1‑91 (MQ=38)
                                                                         tatGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCaa                  <  2:293338/91‑1 (MQ=255)
                                                                         tatGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCaa                  >  2:1974971/1‑91 (MQ=255)
                                                                         tatGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCaa                  >  1:2114448/1‑91 (MQ=39)
                                                                          atGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAAt                 <  1:49372/91‑1 (MQ=255)
                                                                          atGTTTAACATATGACTGTGGCATAGTATTATGCTTATGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAAt                 >  2:3334315/1‑91 (MQ=255)
                                                                          atGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATGATGAAGCCGATAATGAAACCAAt                 >  2:2969928/1‑91 (MQ=38)
                                                                           tGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATc                >  2:752435/1‑91 (MQ=255)
                                                                           tGTTAAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATc                >  2:260554/1‑91 (MQ=255)
                                                                            gTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCg                               >  1:2186368/1‑75 (MQ=255)
                                                                              ttAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCaaa             >  2:1686017/1‑91 (MQ=255)
                                                                               tAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAAc            >  2:2701085/1‑91 (MQ=255)
                                                                               tAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAAc            >  1:1326687/1‑91 (MQ=255)
                                                                                 aCATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTg          >  1:3544351/1‑91 (MQ=255)
                                                                                  cATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTga         >  1:87183/1‑91 (MQ=255)
                                                                                   atatGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTgag        <  1:1878398/91‑1 (MQ=255)
                                                                                     atGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGaa      <  1:1774899/91‑1 (MQ=255)
                                                                                     atGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGaa      >  1:1442999/1‑91 (MQ=255)
                                                                                       gACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGaaca    >  2:1627396/1‑91 (MQ=255)
                                                                                       gACTGTGGCATAGTATTATGCTTGTCGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGaaca    <  1:3622109/91‑1 (MQ=255)
                                                                                        aCTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGaacaa   <  1:430952/91‑1 (MQ=255)
                                                                                         cTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGaacaac  >  2:2198922/1‑91 (MQ=255)
                                                                                         |                                                                                          
TTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAACTGAGAACAAC  >  USA300TCH1516_ALE/216957‑217136

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: