Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 217,027 G→T 6.1% intergenic (‑221/+32) 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,0270GT 6.1% 686.6 / 25.0 198intergenic (‑221/+32)rocD2_1/brnQ_1Ornithine aminotransferase 2/Branched‑chain amino acid transport system 2 carrier protein
Reads supporting (aligned to +/- strand):  ref base G (96/90);  new base T (2/10);  total (98/100)
Fisher's exact test for biased strand distribution p-value = 3.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 204 total aligned reads displayed.

AAATGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAAC  >  USA300TCH1516_ALE/216937‑217117
                                                                                          |                                                                                          
aaaTGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAg                                                                                            <  1:1173638/91‑1 (MQ=255)
aaaTGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAg                                                                                            >  2:858969/1‑91 (MQ=255)
  aTGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCa                                                                                          <  1:3731544/91‑1 (MQ=255)
  aTGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCa                                                                                          >  2:2776689/1‑91 (MQ=255)
  aTGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCa                                                                                          <  2:1999047/91‑1 (MQ=255)
  aTGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCa                                                                                          >  2:1172373/1‑91 (MQ=255)
   tGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCat                                                                                         <  2:1899002/91‑1 (MQ=255)
    gAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCata                                                                                        >  2:3525429/1‑91 (MQ=255)
       aGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGt                                                                                     >  1:3411439/1‑91 (MQ=255)
         tACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGttt                                                                                   <  1:3502628/91‑1 (MQ=255)
         tACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGttt                                                                                   <  1:1601081/91‑1 (MQ=255)
          aCTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTa                                                                                  >  1:3799215/1‑91 (MQ=255)
           cTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTaa                                                                                 <  1:2224189/91‑1 (MQ=255)
            tAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAGTCATATGTTAAAc                                                                                <  1:507515/91‑1 (MQ=255)
            tAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAAc                                                                                >  1:2757179/1‑91 (MQ=255)
             aaTATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACa                                                                               >  1:2031418/1‑91 (MQ=255)
             aaTATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACa                                                                               <  2:2833042/91‑1 (MQ=255)
              atatTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACat                                                                              >  2:1139522/1‑91 (MQ=255)
              atatTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACat                                                                              <  1:391969/91‑1 (MQ=255)
              atatTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACat                                                                              <  1:1668769/91‑1 (MQ=255)
               tatTTTTTGTGGATAATCATTGAGATGTATGTATACGATTGCAGTTTAATGTAGTATGAGTAAAAAATGCCACAAGCATATGTTTAACata                                                                             <  2:2964638/91‑1 (MQ=255)
                 ttttttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATg                                                                           <  2:3718429/91‑1 (MQ=255)
                  tttttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAGTCATATGTTAAACATATgc                                                                          <  2:2018985/91‑2 (MQ=38)
                     ttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAGTCATATGTTAAACATATgcttg                                                                       <  2:1452847/91‑5 (MQ=37)
                     ttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACtg                                                                       <  1:1230138/91‑1 (MQ=255)
                     ttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACtg                                                                       >  2:3676394/1‑91 (MQ=255)
                     ttGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACtg                                                                       <  1:2708056/91‑1 (MQ=255)
                       gtgGATAATCATTGATCTGTATGTATCAGATTGAAGTTTACTGTAGTATGAATAAAAAATGCCACAAGCATATGTTTCACATATGCCtgtg                                                                     >  1:2523509/1‑91 (MQ=255)
                       gtgGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACtgtg                                                                     >  2:2547023/1‑91 (MQ=255)
                        tgGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTgg                                                                    <  2:587766/91‑1 (MQ=255)
                        tagATAATCATTGATAGGCATGTATACGATTGGAGGTGAATGTAGTATGTATAAAAAATGCCACAGTCATATGTTTAACATATGCTTGTgg                                                                    <  2:618323/89‑1 (MQ=21)
                           aTAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCAt                                                                 <  2:1322201/91‑1 (MQ=255)
                            tAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATa                                                                <  1:3256669/91‑1 (MQ=255)
                              aTCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGt                                                              <  2:2293918/91‑1 (MQ=255)
                                cATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATaatac                                                            <  1:1442577/91‑5 (MQ=255)
                                cATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGtat                                                            <  2:1670896/91‑1 (MQ=255)
                                  ttGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGtatta                                                          >  2:2093685/1‑91 (MQ=255)
                                   tGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGtattat                                                         >  1:377731/1‑91 (MQ=255)
                                      tatgtatgTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCt                                                      >  2:1488318/1‑91 (MQ=255)
                                         gtatgtatACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTtgt                                                   <  1:1676148/91‑1 (MQ=255)
                                         gtatgtatACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTtgt                                                   <  1:641637/91‑1 (MQ=255)
                                           atgtatACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTgg                                                 <  1:3582479/91‑1 (MQ=255)
                                           atgtatACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTgg                                                 >  1:3523507/1‑91 (MQ=255)
                                           atgtatACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTgg                                                 <  2:2457226/91‑1 (MQ=255)
                                            tgtatACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGa                                                >  2:1838003/1‑91 (MQ=255)
                                              tataCGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAAt                                              <  1:452946/91‑1 (MQ=255)
                                              tataCGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAAt                                              >  2:2151478/1‑91 (MQ=255)
                                               ataCGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAAtg                                             <  2:3536082/91‑2 (MQ=255)
                                                  cGATTGCAGTTTAATGTAGTATGAATAAAAAATTCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTa                                          >  2:3348837/1‑91 (MQ=255)
                                                  cGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAGCATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTa                                          >  1:3601499/1‑91 (MQ=255)
                                                    aTTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACg                                        <  2:3169060/91‑1 (MQ=255)
                                                       gCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATg                                     <  2:2882203/91‑1 (MQ=255)
                                                       gCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATg                                     <  1:3581008/91‑1 (MQ=255)
                                                       gCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATg                                     >  1:2742949/1‑91 (MQ=255)
                                                       gCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATg                                     >  2:1682572/1‑91 (MQ=255)
                                                       gCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATg                                     <  1:1044264/91‑1 (MQ=255)
                                                        cAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGc                                    <  2:682322/91‑1 (MQ=255)
                                                         aGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCt                                   >  2:1875654/1‑91 (MQ=255)
                                                         aGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCt                                   >  1:2961929/1‑91 (MQ=255)
                                                         aGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCt                                   <  2:1499879/91‑1 (MQ=255)
                                                         aGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCt                                   <  1:542613/91‑1 (MQ=255)
                                                          gTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCtt                                  <  1:2491458/91‑1 (MQ=255)
                                                            ttAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTaa                                >  1:1779630/1‑91 (MQ=255)
                                                            ttAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTaa                                >  1:918183/1‑91 (MQ=255)
                                                              aaTGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAAtt                              >  1:1315542/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:581835/1‑91 (MQ=255)
                                                                 gtagtaTGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTta                           >  2:257792/1‑91 (MQ=255)
                                                                   agtaTGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTtata                         >  1:871349/1‑91 (MQ=255)
                                                                        gAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTataaataa                    <  1:2547023/91‑1 (MQ=255)
                                                                         aaTAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAaataat                   <  2:1554702/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)
                                                                           tAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAAtttt                            <  2:1674336/80‑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:2639915/91‑1 (MQ=255)
                                                                                aaTGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTGATAAATAATGAAGCCg            >  2:2982688/1‑91 (MQ=255)
                                                                                 aTGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGa           <  2:3869464/91‑1 (MQ=255)
                                                                                   gCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATa         >  2:2526989/1‑91 (MQ=255)
                                                                                   gCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATa         >  1:1206145/1‑91 (MQ=255)
                                                                                   gCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATa         >  1:2957554/1‑91 (MQ=255)
                                                                                   gCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATa         >  1:3854400/1‑91 (MQ=255)
                                                                                    ccACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATaa        >  2:3500091/1‑91 (MQ=255)
                                                                                     cacaAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAAt       >  1:705243/1‑91 (MQ=255)
                                                                                     cacaAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAAt       <  1:3385763/91‑1 (MQ=255)
                                                                                      acaAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAg               >  2:2332132/1‑82 (MQ=255)
                                                                                      acaAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAg               <  1:2332132/82‑1 (MQ=255)
                                                                                      acaAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATg      <  2:1427348/91‑1 (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)
                                                                                         ggCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGaaa   <  2:479923/90‑1 (MQ=255)
                                                                                          gCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAAc  <  2:442382/91‑1 (MQ=255)
                                                                                          |                                                                                          
AAATGAAAGTACTAATATTTTTTGTGGATAATCATTGATATGTATGTATACGATTGCAGTTTAATGTAGTATGAATAAAAAATGCCACAAGCATATGTTTAACATATGACTGTGGCATAGTATTATGCTTGTGGAATTTTACGATGCTTAATTTTATAAATAATGAAGCCGATAATGAAAC  >  USA300TCH1516_ALE/216937‑217117

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: