Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 217,036 T→A 12.0% intergenic (‑230/+23) 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,0360TA12.0% 645.3 / 64.7 215intergenic (‑230/+23)rocD2_1/brnQ_1Ornithine aminotransferase 2/Branched‑chain amino acid transport system 2 carrier protein
Reads supporting (aligned to +/- strand):  ref base T (96/93);  new base A (17/9);  total (113/102)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

Only 100 of 221 total aligned reads displayed.

TACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAAC  >  USA300TCH1516_ALE/216946‑217126
                                                                                          |                                                                                          
 aCTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTa                                                                                           >  1:1744380/1‑91 (MQ=255)
    aaTATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACa                                                                                        <  1:2805426/91‑1 (MQ=255)
    aaTATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACa                                                                                        <  2:2833042/91‑1 (MQ=255)
    aaTATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACa                                                                                        >  1:2031418/1‑91 (MQ=255)
     atatTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACat                                                                                       <  1:391969/91‑1 (MQ=255)
     atatTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACat                                                                                       <  2:2615459/91‑1 (MQ=255)
     atatTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACat                                                                                       >  2:1653409/1‑91 (MQ=255)
      tatTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACata                                                                                      <  1:557997/91‑1 (MQ=255)
      tatTTTTTGTGGATAATCATTGAGATGTATGTATACGATTGCAGTTTAATGTAGTATGAGTAAAAAATGCCACAAGCATATGTTTAACata                                                                                      <  2:2964638/91‑1 (MQ=255)
         tttttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAGTCATATGTTAAACATATgc                                                                                   <  2:2018985/91‑2 (MQ=38)
           tttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAGTCATATGTTAAACATATgctt                                                                                 <  1:1418551/91‑4 (MQ=37)
           tttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACt                                                                                 >  1:2730109/1‑91 (MQ=255)
            ttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATGAAAAATGCCTCAGTCATATGTTAAAAGTATgcttg                                                                                <  2:255330/91‑5 (MQ=255)
            ttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAGTCATATGTTAAACATATgcttg                                                                                <  1:1831706/91‑5 (MQ=37)
            ttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAGTCATATGTTAAACATATgcttg                                                                                <  2:1452847/91‑5 (MQ=37)
            ttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAGTCATATGTTAAACATATgcttg                                                                                <  2:932336/91‑5 (MQ=37)
            ttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACtg                                                                                <  1:1230138/91‑1 (MQ=255)
             tgtgGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACtgt                                                                               <  1:788029/91‑1 (MQ=255)
              gtgGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACtgtg                                                                              >  2:2544843/1‑91 (MQ=255)
              gtgGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACtgtg                                                                              >  1:2523506/1‑91 (MQ=255)
              gtgGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACtgtg                                                                              >  1:2876741/1‑91 (MQ=255)
               tgGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTgg                                                                             <  2:587805/91‑1 (MQ=255)
                ggATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGc                                                                            >  2:1774403/1‑91 (MQ=255)
                 gATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAGTCATATGTTAAACATATGCTTGTGGCa                                                                           <  2:1158342/91‑1 (MQ=25)
                   tAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATa                                                                         <  1:3256669/91‑1 (MQ=255)
                   tAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATa                                                                         <  1:540787/91‑1 (MQ=255)
                     aTCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGt                                                                       <  2:2295187/91‑1 (MQ=255)
                      tCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGta                                                                      <  1:3805003/91‑1 (MQ=255)
                        aTTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGtatt                                                                    >  1:652288/1‑91 (MQ=255)
                         ttGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGtatta                                                                   <  1:2489419/91‑1 (MQ=255)
                          tGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGtattat                                                                  >  2:1037224/1‑91 (MQ=255)
                             tatgtatgTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCt                                                               >  2:1488437/1‑91 (MQ=255)
                              atgtatgtATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCtt                                                              >  2:186911/1‑91 (MQ=255)
                               tgtatgtaTACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTAGTATGCAtg                                                             >  2:2372893/1‑91 (MQ=255)
                                gtatgtatACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTtgt                                                            >  2:3043236/1‑91 (MQ=255)
                                  atgtatACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTgg                                                          >  1:3523507/1‑91 (MQ=255)
                                    gtatACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTggca                                                        >  2:150760/1‑89 (MQ=255)
                                     tataCGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAAt                                                       <  1:1863258/91‑1 (MQ=255)
                                      ataCGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAAtt                                                      <  2:2094498/91‑1 (MQ=255)
                                      ataCGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAAtt                                                      >  1:3453039/1‑91 (MQ=255)
                                         cGATTGCAGTTTAATGTAGTATGAATAAAAAATTCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTa                                                   >  2:3348837/1‑91 (MQ=255)
                                         cGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTa                                                   >  1:3793168/1‑91 (MQ=255)
                                         cGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTa                                                   >  2:1672105/1‑91 (MQ=255)
                                           aTTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACg                                                 <  2:3169060/91‑1 (MQ=255)
                                              gCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATg                                              <  1:1044264/91‑1 (MQ=255)
                                              gCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATg                                              >  2:1682572/1‑91 (MQ=255)
                                              gCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATg                                              <  1:3581275/91‑1 (MQ=255)
                                               cAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGc                                             <  2:682322/91‑1 (MQ=255)
                                                 gTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCtt                                           >  1:2234854/1‑91 (MQ=255)
                                                  tttGATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTa                                          <  1:3060651/91‑1 (MQ=255)
                                                   ttAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTaa                                         >  1:918183/1‑91 (MQ=255)
                                                    tAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAAt                                        >  2:2980531/1‑91 (MQ=255)
                                                     aaTGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAAtt                                       <  2:3721768/91‑1 (MQ=255)
                                                      aTGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAAttt                                      <  2:2495467/91‑1 (MQ=255)
                                                       tGTAGTATGAATGAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAAtttt                                     <  2:3226538/91‑1 (MQ=255)
                                                       tGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAAtttt                                     <  1:1653409/91‑1 (MQ=255)
                                                        gtagtaTGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTta                                    <  2:1637017/91‑1 (MQ=255)
                                                        gtagtaTGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTta                                    >  2:49372/1‑91 (MQ=255)
                                                        gtagtaTGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTta                                    >  1:1554702/1‑91 (MQ=255)
                                                           gtaTGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTataa                                 >  2:2863129/1‑91 (MQ=255)
                                                               gAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTataaataa                             <  1:2547023/91‑1 (MQ=255)
                                                                  tAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAAtttt                                     >  1:1674336/1‑80 (MQ=255)
                                                                     aataTGCAACACGCATATGTTTATCTGATCACTGTGCCATATTATTATGCATGTGGAAGTTTACGATGCTTAATTTTATGAATAATGAAGc                       <  2:586407/88‑1 (MQ=255)
                                                                     aaaaTGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGc                       >  1:2844122/1‑91 (MQ=255)
                                                                       taTGCCACAGTCATATGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCg                     >  1:3670441/2‑91 (MQ=25)
                                                                       aaTGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCg                     <  1:2639915/91‑1 (MQ=255)
                                                                        aTGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGa                    <  2:3869464/91‑1 (MQ=255)
                                                                         tGCCACAGTCATATGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGAt                   >  2:1384492/1‑91 (MQ=25)
                                                                         tGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGAt                   >  2:3595361/1‑91 (MQ=255)
                                                                          gCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATa                  >  2:1878398/1‑91 (MQ=255)
                                                                          gCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATa                  >  1:1206145/1‑91 (MQ=255)
                                                                           ccACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATaa                 >  2:3500091/1‑91 (MQ=255)
                                                                           ccACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATaa                 >  1:2961490/1‑91 (MQ=255)
                                                                            cacaAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAAt                <  1:3385763/91‑1 (MQ=255)
                                                                             acaAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATg               >  1:1457850/1‑91 (MQ=255)
                                                                             acaAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATg               <  2:3702993/91‑1 (MQ=255)
                                                                              caggCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGa              <  1:3597925/88‑1 (MQ=255)
                                                                              caAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGa              >  1:3165402/1‑91 (MQ=255)
                                                                               aaGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGaa             >  2:3287710/1‑91 (MQ=255)
                                                                                gtcATATGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGaaa            >  1:3852757/3‑91 (MQ=38)
                                                                                gtcATATGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGaaa            >  1:1141397/3‑91 (MQ=38)
                                                                                gtcATATGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGaaa            >  2:2886059/3‑91 (MQ=38)
                                                                                gtcATATGTTAAACATATGCTTGTGGCATAGCATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGaaa            >  1:119353/3‑91 (MQ=37)
                                                                                aGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGaaa            <  2:479994/91‑1 (MQ=255)
                                                                                aGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGaaa            <  2:479901/91‑1 (MQ=255)
                                                                                 tcATATGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAAc           >  1:2749268/2‑91 (MQ=38)
                                                                                  cATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAAcc          <  1:2675466/91‑1 (MQ=255)
                                                                                  cATATGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAAcc          >  1:1452162/1‑91 (MQ=38)
                                                                                   atatGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCa         <  2:3793168/91‑1 (MQ=255)
                                                                                   atatGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCa         >  2:2243993/1‑91 (MQ=38)
                                                                                    tatGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCaa        <  2:342373/91‑1 (MQ=255)
                                                                                    tatGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCaa        >  2:1975306/1‑91 (MQ=255)
                                                                                    tatGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCaa        <  2:705243/91‑1 (MQ=255)
                                                                                    tatGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCaa        >  1:2114448/1‑91 (MQ=39)
                                                                                     atGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAAt       <  1:49372/91‑1 (MQ=255)
                                                                                      tGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATc      >  1:360212/1‑91 (MQ=255)
                                                                                       gTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCg                     <  2:2186368/75‑1 (MQ=255)
                                                                                       gTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCg                     >  1:2186368/1‑75 (MQ=255)
                                                                                          tAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAAc  >  2:2701085/1‑91 (MQ=255)
                                                                                          tAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAAc  >  2:775473/1‑91 (MQ=255)
                                                                                          |                                                                                          
TACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAACCAATCAAAC  >  USA300TCH1516_ALE/216946‑217126

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