Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,105,994 T→C T12T (ACT→ACC pdhD → Dihydrolipoyl dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,105,9940TC99.2% 485.0 / ‑4.6 134T12T (ACT→ACCpdhDDihydrolipoyl dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (75/58);  total (76/58)
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.64e-01

Only 100 of 138 total aligned reads displayed.

ACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGG  >  USA300TCH1516_ALE/1105906‑1106080
                                                                                        |                                                                                      
accttAACCGTTTATTCGATAATCCAGAAACCATATTAACGGAGGGGCAAAACATGGTAGTTGGAGAATTTCCAAGGGACACAGACACtac                                                                                      >  2:2363242/4‑90 (MQ=255)
 caTTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACcata                                                                                     >  1:112026/1‑91 (MQ=255)
 caTTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACcata                                                                                     >  2:199328/1‑91 (MQ=255)
  aTTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAg                                                                                    <  1:1139630/91‑1 (MQ=255)
  aTTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAg                                                                                    <  2:1522602/91‑1 (MQ=255)
   ttAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGt                                                                                   <  1:2970671/91‑1 (MQ=255)
   ttAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGt                                                                                   <  2:2979392/91‑1 (MQ=255)
   ttAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGt                                                                                   <  2:2979351/91‑1 (MQ=255)
    tAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTa                                                                                  <  1:1082504/91‑1 (MQ=255)
    tAAACGTTTATTAAATAATCCAGAATTAGTATTAATGGAGGGGTAAAAAATGGTAGCTGGCGATTTCCCAATTCAAACAGATACCCTAGTa                                                                                  >  2:2058286/1‑91 (MQ=255)
      aaCGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAATACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAAt                                                                                <  2:2774871/91‑1 (MQ=255)
      aaCGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAAt                                                                                <  1:759581/91‑1 (MQ=255)
        cGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTCAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCg                                                                              >  2:2545823/1‑91 (MQ=255)
         gTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCgg                                                                             >  1:1655518/1‑91 (MQ=255)
         gTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCgg                                                                             >  1:1655850/1‑91 (MQ=255)
         gTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCgg                                                                             >  2:1776777/1‑91 (MQ=255)
           ttattaAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAg                                                                           >  2:1172925/1‑91 (MQ=255)
           ttattaAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAg                                                                           >  2:2179584/1‑91 (MQ=255)
           ttattaAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAg                                                                           >  1:1828942/1‑91 (MQ=255)
           ttattaAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAg                                                                           >  2:554462/1‑91 (MQ=255)
             attaAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCa                                                                         <  2:813129/91‑1 (MQ=255)
              ttaAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAg                                                                        >  2:2219427/1‑91 (MQ=255)
               taAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAgg                                                                       >  2:1884597/1‑91 (MQ=255)
                 aataatCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTc                                                                     <  1:1814489/91‑1 (MQ=255)
                    aatCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCtg                                                                  >  2:2199154/1‑91 (MQ=255)
                    aatCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCtg                                                                  >  2:1704137/1‑91 (MQ=255)
                        cAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCtggtgg                                                              >  2:1236289/1‑91 (MQ=255)
                        cAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTCGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCtggtgg                                                              >  1:2664947/1‑91 (MQ=255)
                         aGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGa                                                             >  2:2688582/1‑91 (MQ=255)
                         aGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGa                                                             <  2:1707140/91‑1 (MQ=255)
                          gAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGAt                                                            <  2:1486835/91‑1 (MQ=255)
                          gAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGAt                                                            <  2:1320089/91‑1 (MQ=255)
                          gAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGAt                                                            >  2:2070221/1‑91 (MQ=255)
                          gAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGAt                                                            <  1:1619115/91‑1 (MQ=255)
                          gAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGAt                                                            >  1:2503778/1‑91 (MQ=255)
                          gAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGAt                                                            >  1:2503769/1‑91 (MQ=255)
                             ttattattaATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACg                                                         >  2:942984/1‑91 (MQ=255)
                              tattattaATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGt                                                        <  1:2237503/91‑1 (MQ=255)
                                 tattaATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTgc                                                     <  1:1461264/91‑1 (MQ=255)
                                 tattaATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTgc                                                     <  1:1461374/91‑1 (MQ=255)
                                  attaATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTgca                                                    >  1:54470/1‑91 (MQ=255)
                                  attaATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTgca                                                    >  1:2943249/1‑91 (MQ=255)
                                   ttaATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTgcag                                                   <  1:1716273/91‑1 (MQ=255)
                                    taATGGGGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTgcagc                                                  <  1:822168/91‑1 (MQ=255)
                                    taATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTgcagc                                                  <  1:1871174/91‑1 (MQ=255)
                                    taATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTgcagc                                                  <  1:821875/91‑1 (MQ=255)
                                       tGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAAt                                               >  2:1236765/1‑91 (MQ=255)
                                       tGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAAt                                               >  2:1236314/1‑91 (MQ=255)
                                           ggggTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGt                                           >  2:457463/1‑91 (MQ=255)
                                           ggggTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGt                                           >  2:2224521/1‑91 (MQ=255)
                                            gggTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTg                                          >  1:980181/1‑91 (MQ=255)
                                            gggTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTg                                          >  1:93192/1‑91 (MQ=255)
                                            gggTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTg                                          >  1:560482/1‑91 (MQ=255)
                                             ggTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGc                                         <  1:527187/91‑1 (MQ=255)
                                               tAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAg                                       <  1:2843743/91‑1 (MQ=255)
                                                aaaaCATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGc                                      <  1:2437161/91‑1 (MQ=255)
                                                 aaaCATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCt                                     <  2:2429667/91‑1 (MQ=255)
                                                   aCATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCa                                   >  2:603446/1‑91 (MQ=255)
                                                     aTGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAAt                                 <  1:958110/91‑1 (MQ=255)
                                                     aTGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAAt                                 <  2:2291977/91‑1 (MQ=255)
                                                       ggTAGTTGGAGATTTCCCAATTGAAGCAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTa                               >  1:1927846/1‑91 (MQ=255)
                                                       ggTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTa                               >  2:415199/1‑91 (MQ=255)
                                                       ggTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTa                               >  2:345123/1‑91 (MQ=255)
                                                        gTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAg                              >  2:2200665/1‑91 (MQ=255)
                                                        gTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAg                              >  2:2482110/1‑91 (MQ=255)
                                                        gTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAg                              <  1:817532/91‑1 (MQ=255)
                                                         tAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAgg                             >  2:1145887/1‑91 (MQ=255)
                                                         tAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAgg                             >  2:2313588/1‑91 (MQ=255)
                                                         tAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAgg                             >  2:2313567/1‑91 (MQ=255)
                                                          aGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGa                            >  1:1262172/1‑91 (MQ=255)
                                                             tGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACaa                         >  1:796747/1‑91 (MQ=255)
                                                              ggaaaTTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACaaa                        >  2:2274659/5‑91 (MQ=255)
                                                              ggAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACaaa                        >  1:2812982/1‑91 (MQ=255)
                                                              ggAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACaaa                        >  2:2602542/1‑91 (MQ=255)
                                                                   tttCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTa                   >  2:1516581/1‑91 (MQ=255)
                                                                    ttCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTaa                  <  1:911525/91‑1 (MQ=255)
                                                                     tCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAAc                 <  1:2619224/91‑1 (MQ=255)
                                                                     tCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAAc                 <  2:646701/91‑1 (MQ=255)
                                                                     tCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAAc                 <  1:2382093/91‑1 (MQ=255)
                                                                     tCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAAc                 <  1:480507/91‑1 (MQ=255)
                                                                       ccAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACaa               >  1:1250919/1‑91 (MQ=255)
                                                                        cAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAAt              >  2:952009/1‑91 (MQ=255)
                                                                        cAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAAt              >  2:1136259/1‑91 (MQ=255)
                                                                        cAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAAt              <  1:1864452/91‑1 (MQ=255)
                                                                         aaTTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATc             <  2:906232/91‑1 (MQ=255)
                                                                          aTTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCg            >  1:108181/1‑91 (MQ=255)
                                                                          aTTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCg            <  2:898943/91‑1 (MQ=255)
                                                                          aTTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCg            <  1:1850019/91‑1 (MQ=255)
                                                                          aTTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATACGTGCAGCTCAATTAGGACAAAAAGTAACAATCg            <  1:13216/91‑1 (MQ=255)
                                                                            tGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGtt          >  1:27680/1‑91 (MQ=255)
                                                                             gAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTg         <  1:607385/91‑1 (MQ=255)
                                                                               aaCAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTgag       <  1:1753722/91‑1 (MQ=255)
                                                                               aaCAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTgag       >  1:1164542/1‑91 (MQ=255)
                                                                                aCAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTgaga      >  2:1355909/1‑91 (MQ=255)
                                                                                 cAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGaa     <  1:2961165/91‑1 (MQ=255)
                                                                                  aGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGaaa    >  2:1025102/1‑91 (MQ=255)
                                                                                  aGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGaaa    >  2:1865441/1‑91 (MQ=255)
                                                                                  aGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGaaa    >  2:2258790/1‑91 (MQ=255)
                                                                                   gATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAg   >  1:1178850/1‑91 (MQ=255)
                                                                                    aTACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCCATTAGGACAAAAAGTAACAATCGTTGAGAAAgg  >  2:2586969/1‑91 (MQ=255)
                                                                                        |                                                                                      
ACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGG  >  USA300TCH1516_ALE/1105906‑1106080

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: