Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,398,086 A→T 5.2% intergenic (‑247/‑43) modA ← / → fdhD Molybdate‑binding periplasmic protein/Sulfurtransferase FdhD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,398,0860AT 5.2% 858.5 / 29.3 251intergenic (‑247/‑43)modA/fdhDMolybdate‑binding periplasmic protein/Sulfurtransferase FdhD
Reads supporting (aligned to +/- strand):  ref base A (105/133);  new base T (8/5);  total (113/138)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

Only 100 of 257 total aligned reads displayed.

CGAATTCAATTTAAATCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCAACCAATTGTTAGATATGAAGAT  >  USA300TCH1516_ALE/2397996‑2398176
                                                                                           |                                                                                          
cGAATTCAATTTAAATCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCaaa                                                                                            >  2:433662/1‑91 (MQ=255)
 gAATTCAATTTAAATCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAAt                                                                                           >  1:2594856/1‑91 (MQ=255)
 gAATTCAATTTAAATCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAAt                                                                                           <  1:1798413/91‑1 (MQ=255)
 aaaTTCAATTTAAATCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAAt                                                                                           <  2:3264462/90‑1 (MQ=255)
  aaTTCAATTTAAATCAAACATATTACTCTCTATTAATTTTTATCAATCAAATTTCAAACTTAATCAAATATATTAAAATTTTAATCAAAtt                                                                                          <  2:367058/91‑1 (MQ=255)
  aaTTCAATTTAAATCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAAtt                                                                                          <  2:2545651/91‑1 (MQ=255)
    ttCAATTTAAATCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGc                                                                                        >  1:2093359/1‑91 (MQ=255)
     tCAATTTAAATCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCt                                                                                       <  2:3740733/91‑1 (MQ=255)
      cAATTTAAATCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAGCTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCtt                                                                                      <  1:3352512/91‑1 (MQ=255)
      cAATTTAAATCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCtt                                                                                      <  1:799399/91‑1 (MQ=255)
       aaTTTAAATCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCttt                                                                                     >  2:3402791/1‑91 (MQ=255)
        aTTTAAATCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTa                                                                                    >  1:1538391/1‑91 (MQ=255)
         tttAAATCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATATTAAAATTTTAAAGCAATTTGCttaa                                                                                    <  2:2248564/91‑3 (MQ=255)
          ttAAATCAAACAGATTGCTCGCTATTAATTTTTATCTATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTaaa                                                                                  <  2:3904431/91‑1 (MQ=255)
          ttAAATCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAAGCAATTTGCttaaa                                                                                   <  2:3571855/91‑4 (MQ=255)
           tAAATCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTaaaa                                                                                 <  1:3438800/91‑1 (MQ=255)
            aaaTCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAAt                                                                                <  1:1989788/91‑1 (MQ=255)
            aaaTCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAAt                                                                                <  1:1697784/91‑1 (MQ=255)
             aaTCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATa                                                                               >  1:564766/1‑91 (MQ=255)
              aTCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATaa                                                                              <  2:277535/91‑1 (MQ=255)
              aTCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATaa                                                                              <  1:2928898/91‑1 (MQ=255)
                cAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAg                                                                            <  2:691879/91‑1 (MQ=255)
                  aaCAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGtt                                                                          >  1:469420/1‑91 (MQ=255)
                  aaCAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGtt                                                                          <  2:2497115/91‑1 (MQ=255)
                   aCAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGttt                                                                         <  2:32316/91‑1 (MQ=255)
                     aGATTGCTCGCTATTAATTTTTATCAATCAACTTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTaa                                                                       >  2:2598209/1‑91 (MQ=255)
                     aGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTaa                                                                       >  1:2999862/1‑91 (MQ=255)
                      gATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAg                                                                      <  1:1231588/91‑1 (MQ=255)
                      gATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAg                                                                      <  1:3774890/91‑1 (MQ=255)
                       aTTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGc                                                                     <  2:3135985/91‑1 (MQ=255)
                       aTTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGc                                                                     <  1:3282000/91‑1 (MQ=255)
                       aTTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGc                                                                     >  1:1891289/1‑91 (MQ=255)
                         tGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTa                                                                   <  2:270596/91‑1 (MQ=255)
                         tGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTa                                                                   <  1:1674492/91‑1 (MQ=255)
                          gCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTaa                                                                  >  1:1643254/1‑91 (MQ=255)
                           cTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAAt                                                                 <  2:2386174/91‑1 (MQ=255)
                           cTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAAt                                                                 >  2:1234234/1‑91 (MQ=255)
                           cTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAAt                                                                 <  2:2774857/91‑1 (MQ=255)
                            tCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAAtt                                                                >  2:3845636/1‑91 (MQ=255)
                            tCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAAtt                                                                <  1:863335/91‑1 (MQ=255)
                             cGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTa                                                               <  2:1413352/91‑1 (MQ=255)
                              gCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTaa                                                              >  2:137923/1‑91 (MQ=255)
                                tATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAg                                                            <  1:2290550/91‑1 (MQ=255)
                                tATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAg                                                            >  1:212669/1‑91 (MQ=255)
                                   tAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAgggg                                                         <  1:603991/91‑1 (MQ=255)
                                    aaTTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGt                                                        <  1:1649797/91‑1 (MQ=255)
                                    aaTTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGt                                                        <  1:704456/91‑1 (MQ=255)
                                    aaTTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGt                                                        <  1:3621905/91‑1 (MQ=255)
                                    aaTTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGt                                                        >  1:1016424/1‑91 (MQ=255)
                                     aTTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGtt                                                       >  2:1761284/1‑91 (MQ=255)
                                       ttttATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAg                                                     <  2:286826/91‑1 (MQ=255)
                                        tttATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGa                                                    >  1:2876248/1‑91 (MQ=255)
                                        tttATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGa                                                    >  2:496638/1‑91 (MQ=255)
                                           atcaatcaAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATaa                                                 >  1:2792656/1‑91 (MQ=255)
                                             caatcaaATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATg                                               <  2:113433/91‑1 (MQ=255)
                                              aatcaaATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGa                                              >  1:270201/1‑91 (MQ=255)
                                               atcaaATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGaa                                             >  1:1798248/1‑91 (MQ=255)
                                               atcaaATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGaa                                             <  2:2902220/91‑1 (MQ=255)
                                                 aaaaTTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACa                                           <  1:2515319/90‑1 (MQ=255)
                                                  aaaTTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACaa                                          <  1:1106691/91‑1 (MQ=255)
                                                  aaaTTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACaa                                          >  1:1483172/1‑91 (MQ=255)
                                                  aaaTTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACaa                                          >  2:1651013/1‑91 (MQ=255)
                                                    aTTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAg                                        >  1:4122/1‑91 (MQ=255)
                                                     tttCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGa                                       <  2:3770122/91‑1 (MQ=255)
                                                      ttCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGAt                                      >  2:2922857/1‑91 (MQ=255)
                                                         aaaCTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTa                                   >  1:192318/1‑91 (MQ=255)
                                                         aaaCTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTa                                   <  2:2325305/91‑1 (MQ=255)
                                                          aaCTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTAt                                  <  2:2433395/91‑1 (MQ=255)
                                                            cTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAACGATGTATCt                                >  2:3956082/1‑91 (MQ=255)
                                                                aGCAAATATGTTAAAATTTTAAGCCAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAg                            >  2:881967/1‑91 (MQ=255)
                                                                aGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAg                            >  2:106359/1‑91 (MQ=255)
                                                                 gCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTATAGGGGTTAGATAATGAACAAAGATGTATCTTTAgg                           <  2:837078/91‑1 (MQ=255)
                                                                 gCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAgg                           >  1:95898/1‑91 (MQ=255)
                                                                 gCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAgg                           >  2:2307638/1‑91 (MQ=255)
                                                                 gCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAgg                           <  1:3930487/91‑1 (MQ=255)
                                                                  cAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGt                          <  2:1933477/91‑1 (MQ=255)
                                                                   aaaTATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTc                         >  2:335175/1‑91 (MQ=255)
                                                                    aaTATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCa                        <  1:3470169/91‑1 (MQ=255)
                                                                    aaTATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCa                        >  2:182197/1‑91 (MQ=255)
                                                                     atatGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGATAGATAATGAACAAAGACGTATCTTTAGGTCaa                       >  2:1571972/1‑91 (MQ=255)
                                                                      tatGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCAAc                      >  2:953992/1‑91 (MQ=255)
                                                                      tatGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCAAc                      <  2:3528727/91‑1 (MQ=255)
                                                                      tatGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCAAc                      <  1:2357332/91‑1 (MQ=255)
                                                                      tatGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAANTAAAGGGGTTAGATAATGANCAAAGATGTATCTTTAGGTCAAc                      >  2:1456917/1‑91 (MQ=255)
                                                                      tatGTTAAAAGTTTAAGCAAATTGCTTTAAAATAAAGTTTCAGCTAATTGAGGGGCTTAGATATGGAACAAAGATGTATCTTTAGGGCAAc                      <  2:3313896/91‑1 (MQ=255)
                                                                       atGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCAAcc                     <  2:3275834/91‑1 (MQ=255)
                                                                           tAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCAACCAAtt                 <  1:750454/91‑1 (MQ=255)
                                                                             aaaTTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCAACCAATTGt               >  2:1807545/1‑91 (MQ=255)
                                                                              aaTTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCAACCAATTGtt              <  1:2909160/91‑1 (MQ=255)
                                                                              aaTTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCAACCAATTGtt              >  2:2334515/1‑91 (MQ=255)
                                                                               aTTTTAAAGCAATTTGCTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCAACCAATTGTTa             >  1:437354/1‑91 (MQ=38)
                                                                                   tAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAGGATGTATCTTTAGGTCAACCAATTGTTAGata         >  2:909668/1‑91 (MQ=255)
                                                                                   tAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAGGATGTATCTTTAGGTCAACCAATTGTTAGata         >  2:1673610/1‑91 (MQ=255)
                                                                                   tAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCAACCAATTGTTAGata         >  1:1709727/1‑91 (MQ=255)
                                                                                     aaGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCAACCAATTGTTAGatat        >  1:2424704/1‑91 (MQ=255)
                                                                                      aGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCAACCAATTGTTAGATATg       <  1:3577758/91‑1 (MQ=255)
                                                                                      aGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCAACCAATTGTTAGATATg       >  1:11692/1‑91 (MQ=255)
                                                                                       gCAATTTGCTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCAACCAATTGTTAGATATGaa     >  1:3267876/1‑91 (MQ=255)
                                                                                         aaaTTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCAACCAATTGTTAGATATGAAg    >  2:2004034/1‑91 (MQ=255)
                                                                                           aTTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCAACCAATTGTTAGATATGAAGAt  >  2:2872984/1‑91 (MQ=255)
                                                                                           |                                                                                          
CGAATTCAATTTAAATCAAACAGATTGCTCGCTATTAATTTTTATCAATCAAATTTCAAACTTAAGCAAATATGTTAAAATTTTAAGCAAATTGCTTTAAAATAAAGTTTAAGCTAATTAAAGGGGTTAGATAATGAACAAAGATGTATCTTTAGGTCAACCAATTGTTAGATATGAAGAT  >  USA300TCH1516_ALE/2397996‑2398176

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: