Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 217,026 A→G 6.2% intergenic (‑220/+33) 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,0260AG 6.2% 602.1 / 19.0 194intergenic (‑220/+33)rocD2_1/brnQ_1Ornithine aminotransferase 2/Branched‑chain amino acid transport system 2 carrier protein
Reads supporting (aligned to +/- strand):  ref base A (93/88);  new base G (2/10);  total (95/99)
Fisher's exact test for biased strand distribution p-value = 3.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 203 total aligned reads displayed.

AAAATGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAA  >  USA300TCH1516_ALE/216936‑217116
                                                                                          |                                                                                          
 aaaTGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAg                                                                                           <  1:1173638/91‑1 (MQ=255)
 aaaTGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAg                                                                                           >  2:858969/1‑91 (MQ=255)
 aaaTGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAg                                                                                           <  1:3531820/91‑1 (MQ=255)
   aTGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCa                                                                                         >  2:2776689/1‑91 (MQ=255)
   aTGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCa                                                                                         <  1:3731544/91‑1 (MQ=255)
   aTGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCa                                                                                         <  2:1999047/91‑1 (MQ=255)
   aTGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCa                                                                                         <  2:2053396/91‑1 (MQ=255)
    tGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCat                                                                                        <  2:1437977/91‑1 (MQ=255)
    tGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCat                                                                                        <  2:1899002/91‑1 (MQ=255)
     gAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCata                                                                                       >  2:3525429/1‑91 (MQ=255)
        aGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACGAGCATATGt                                                                                    <  2:118434/91‑1 (MQ=255)
        aGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGt                                                                                    >  1:3411439/1‑91 (MQ=255)
          tACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGttt                                                                                  <  1:1601081/91‑1 (MQ=255)
          tACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGttt                                                                                  <  1:3502628/91‑1 (MQ=255)
           aCTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTa                                                                                 >  1:3799215/1‑91 (MQ=255)
             tAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAGTCATATGTTAAAc                                                                               <  1:507515/91‑1 (MQ=255)
              aaTATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACa                                                                              >  1:2031418/1‑91 (MQ=255)
              aaTATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACa                                                                              <  1:2805426/91‑1 (MQ=255)
              aaTATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACa                                                                              <  2:2833042/91‑1 (MQ=255)
              aaTATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACa                                                                              >  1:2549051/1‑91 (MQ=255)
               atatTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACat                                                                             <  2:2615459/91‑1 (MQ=255)
               atatTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACat                                                                             <  1:391969/91‑1 (MQ=255)
               atatTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACat                                                                             <  1:1668769/91‑1 (MQ=255)
               atatTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACat                                                                             >  2:1653409/1‑91 (MQ=255)
                tatTTTTTGTGGATAATCATTGAGATGTATGTATACGATTGCAGTTTAATGTAGTATGAGTAAAAAATGCCACAAGCATATGTTTAACata                                                                            <  2:2964638/91‑1 (MQ=255)
                  ttttttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATg                                                                          <  2:3718429/91‑1 (MQ=255)
                     tttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACt                                                                       >  1:2730109/1‑91 (MQ=255)
                      ttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATGAAAAATGCCTCAGTCATATGTTAAAAGTATgcttg                                                                      <  2:255330/91‑5 (MQ=255)
                      ttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACtg                                                                      <  1:1230138/91‑1 (MQ=255)
                      ttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACtg                                                                      >  2:3676394/1‑91 (MQ=255)
                        gtgGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACtgtg                                                                    >  2:2544843/1‑91 (MQ=255)
                        gtgGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACtgtg                                                                    >  1:2523506/1‑91 (MQ=255)
                         tgGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTgg                                                                   <  2:587766/91‑1 (MQ=255)
                         tagATAATCATTGATAGGCATGTATACGATTGGAGGTGAATGTAGTATGTATAAAAAATGCCACAGTCATATGTTTAACATATGCTTGTgg                                                                   <  2:618323/89‑1 (MQ=21)
                          ggATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGc                                                                  >  2:1774403/1‑91 (MQ=255)
                           gATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAGTCATATGTTAAACATATGCTTGTGGCa                                                                 <  2:1158342/91‑1 (MQ=25)
                             tAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATa                                                               <  1:3256669/91‑1 (MQ=255)
                               aTCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGt                                                             <  2:2295187/91‑1 (MQ=255)
                                 cATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATaatac                                                           <  1:1442577/91‑5 (MQ=255)
                                 cATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGtat                                                           <  1:274393/91‑1 (MQ=255)
                                   ttGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGtatta                                                         <  1:2489419/91‑1 (MQ=255)
                                    tGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGtattat                                                        >  1:377731/1‑91 (MQ=255)
                                    tGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGtattat                                                        >  2:1037224/1‑91 (MQ=255)
                                       tatgtatgTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCt                                                     >  2:1488437/1‑91 (MQ=255)
                                         tgtatgtaTACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTAGTATGCAtg                                                   >  2:2372893/1‑91 (MQ=255)
                                          gtatgtatACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTtgt                                                  <  1:1676148/91‑1 (MQ=255)
                                          gtatgtatACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTtgt                                                  <  1:641637/91‑1 (MQ=255)
                                             tgtatACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGa                                               >  2:1838003/1‑91 (MQ=255)
                                              gtatACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTggca                                              >  2:150760/1‑89 (MQ=255)
                                               tataCGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAAt                                             <  1:452946/91‑1 (MQ=255)
                                                ataCGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAAtt                                            <  2:2094498/91‑1 (MQ=255)
                                                ataCGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAAtg                                            <  2:3536082/91‑2 (MQ=255)
                                                   cGATTGCAGTTTAATGTAGTATGAATAAAAAATTCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTa                                         >  2:3348837/1‑91 (MQ=255)
                                                   cGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAGCATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTa                                         >  1:3601499/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)
                                                        gCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATg                                    <  2:2882203/91‑1 (MQ=255)
                                                        gCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATg                                    <  1:1044264/91‑1 (MQ=255)
                                                        gCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATg                                    <  1:3581275/91‑1 (MQ=255)
                                                         cAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGc                                   <  2:682322/91‑1 (MQ=255)
                                                          aGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCt                                  >  1:2961929/1‑91 (MQ=255)
                                                          aGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCt                                  >  2:1875654/1‑91 (MQ=255)
                                                          aGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCt                                  <  1:542613/91‑1 (MQ=255)
                                                           gTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCtt                                 >  1:2234854/1‑91 (MQ=255)
                                                             ttAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTaa                               >  1:918183/1‑91 (MQ=255)
                                                             ttAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTaa                               >  1:1779630/1‑91 (MQ=255)
                                                              tAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAAt                              >  1:2184598/1‑91 (MQ=255)
                                                               aaTGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAAtt                             >  1:1315542/1‑91 (MQ=255)
                                                               aaTGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAAtt                             >  2:3429562/1‑91 (MQ=255)
                                                                aTGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAAttt                            <  2:2495467/91‑1 (MQ=255)
                                                                 tGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAAtttt                           <  1:1653409/91‑1 (MQ=255)
                                                                  gtagtaTGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTta                          >  2:2247250/1‑91 (MQ=255)
                                                                  gtagtaTGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTta                          >  2:581835/1‑91 (MQ=255)
                                                                  gtagtaTGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTta                          <  2:1637017/91‑1 (MQ=255)
                                                                    agtaTGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTtata                        >  1:871349/1‑91 (MQ=255)
                                                                         gAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTataaataa                   <  1:2544843/91‑1 (MQ=255)
                                                                          aaTAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAaataat                  <  2:321963/91‑1 (MQ=255)
                                                                          aaTAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAaataat                  <  1:2093685/91‑1 (MQ=255)
                                                                           aTAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATg                 <  1:2863701/91‑1 (MQ=255)
                                                                             aaaaaaTGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGaa               <  2:377731/91‑1 (MQ=255)
                                                                                aaaTGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGcc            >  1:15117/1‑91 (MQ=255)
                                                                                 taTGCCACAGTCATATGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCg           >  1:3670441/2‑91 (MQ=25)
                                                                                 aaTGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCg           <  1:2640877/91‑1 (MQ=255)
                                                                                 aaTGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTGATAAATAATGAAGCCg           >  2:2982688/1‑91 (MQ=255)
                                                                                  aTGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGa          <  2:3869464/91‑1 (MQ=255)
                                                                                   tGCCACAGTCATATGTTAAACATATGCTTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGAt         >  2:1384492/1‑91 (MQ=25)
                                                                                   tGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGAt         >  2:1681355/1‑91 (MQ=255)
                                                                                    gCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATa        >  1:1206145/1‑91 (MQ=255)
                                                                                    gCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATa        >  1:3854400/1‑91 (MQ=255)
                                                                                    gCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATa        >  1:2957554/1‑91 (MQ=255)
                                                                                     ccACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATaa       >  2:3500091/1‑91 (MQ=255)
                                                                                      cacaAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAAt      <  1:3385763/91‑1 (MQ=255)
                                                                                      cacaAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAAt      >  1:705243/1‑91 (MQ=255)
                                                                                       acaAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATg     >  2:493129/1‑91 (MQ=255)
                                                                                       acaAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATg     >  1:1457849/1‑91 (MQ=255)
                                                                                        caAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGa    >  1:1859987/1‑91 (MQ=255)
                                                                                        caAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGa    >  1:3165402/1‑91 (MQ=255)
                                                                                         aaGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGaa   >  2:3287055/1‑91 (MQ=255)
                                                                                          aGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGaaa  <  2:479901/91‑1 (MQ=255)
                                                                                          |                                                                                          
AAAATGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAA  >  USA300TCH1516_ALE/216936‑217116

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: