Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,744,205 A→G 10.8% intergenic (+142/+51) fda → / ← mqo2 Fructose‑bisphosphate aldolase class 1/putative malate:quinone oxidoreductase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,744,2050AG10.8% 627.2 / 52.5 223intergenic (+142/+51)fda/mqo2Fructose‑bisphosphate aldolase class 1/putative malate:quinone oxidoreductase 2
Reads supporting (aligned to +/- strand):  ref base A (115/84);  new base G (7/17);  total (122/101)
Fisher's exact test for biased strand distribution p-value = 9.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

Only 100 of 238 total aligned reads displayed.

TGTGCTTTAAAGGCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAACCTAATTCTAAGTCTTTTGAAGTT  >  USA300TCH1516_ALE/2744115‑2744295
                                                                                          |                                                                                          
tgtgCTTTAAAGGCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGaa                                                                                            >  1:791265/1‑91 (MQ=255)
tgtgCTTTAAAGGCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGaa                                                                                            >  1:2890352/1‑91 (MQ=255)
tgtgCTTTAAAGGCATTGTGTATGAGACTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGaa                                                                                            >  2:2219564/1‑91 (MQ=255)
 gtgCTTTAAAGGCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGaaa                                                                                           <  1:82961/91‑1 (MQ=255)
 gtgCTTTAAAGGCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGaaa                                                                                           >  1:2237962/1‑91 (MQ=255)
    cTTTAAAGGCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAAt                                                                                        <  2:2507267/91‑1 (MQ=255)
    cTTTAAAGGCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAAt                                                                                        >  2:2142475/1‑91 (MQ=255)
    cTTTAAAGGCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAAt                                                                                        >  2:1702074/1‑91 (MQ=255)
    cTTTAAAGGCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAAt                                                                                        >  2:1233053/1‑91 (MQ=255)
    cTTTAAAGGCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAAt                                                                                        >  1:3147778/1‑91 (MQ=255)
    cTTTAAAGGCATTGTGTATGAGTCTTTTTTGATTAAACTTGACTTTGCCGCATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAAt                                                                                        >  1:1110173/1‑91 (MQ=255)
     tttAAAGGCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATc                                                                                       >  1:3648219/1‑91 (MQ=255)
      ttAAAGGCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCt                                                                                      <  2:1431290/91‑1 (MQ=255)
      ttAAAGGCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCt                                                                                      <  2:2131100/91‑1 (MQ=255)
      ttAAAGGCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCt                                                                                      >  2:1826356/1‑91 (MQ=255)
      ttAAAGGCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCt                                                                                      >  2:597756/1‑91 (MQ=255)
        aaaGGCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCttt                                                                                    >  2:3835000/1‑91 (MQ=255)
          aGGCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTaa                                                                                  <  1:3354083/91‑1 (MQ=255)
          aGGCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTaa                                                                                  <  1:432823/91‑1 (MQ=255)
           ggCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAg                                                                                 >  2:1539098/1‑91 (MQ=255)
           ggCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAg                                                                                 <  1:2970905/91‑1 (MQ=255)
           ggCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAACTGAGAAATCTTAAAg                                                                                 <  1:2373347/91‑1 (MQ=38)
             cATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGAt                                                                               >  1:2975399/1‑91 (MQ=255)
              aTTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGAtt                                                                              >  2:528834/1‑91 (MQ=255)
               ttGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGAttt                                                                             <  1:1929953/91‑1 (MQ=255)
               ttGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGAttt                                                                             >  2:3069709/1‑91 (MQ=255)
                tgtgTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTtc                                                                            >  2:24416/1‑91 (MQ=255)
                tgtgTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTtc                                                                            <  2:1897646/91‑1 (MQ=255)
                 gtgtATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTtct                                                                           >  1:1213641/1‑91 (MQ=255)
                  tgtATGAGTCTTTTTTGTTTAATGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTtctc                                                                          >  2:504201/1‑91 (MQ=255)
                  tgtATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTCAGTAAAAAAATCGGACAGCTGAAAAATCTTAAAGATTtctc                                                                          >  2:3500136/1‑91 (MQ=255)
                  tgtATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAGAAATCTTAAAGAtttttc                                                                          <  1:3162824/91‑4 (MQ=38)
                   gtATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCa                                                                         >  2:1695794/1‑91 (MQ=255)
                   gtATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCa                                                                         >  2:2877354/1‑91 (MQ=255)
                    tATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAACTGAGAAATCTTAAAGATTTTTCAg                                                                        <  2:639907/91‑1 (MQ=25)
                     aTGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGt                                                                       >  2:3330557/1‑91 (MQ=255)
                     aTGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGt                                                                       >  2:61225/1‑91 (MQ=255)
                     aTGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGt                                                                       <  1:3468563/91‑1 (MQ=255)
                       gAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTg                                                                     <  1:1706730/91‑1 (MQ=255)
                       gAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAACTGAGAAATCTTAAAGATTTTTCAgctg                                                                     <  1:3549081/91‑4 (MQ=21)
                        aGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAACTGAGAAATCTTAAAGATTTTTCAGCTGt                                                                    <  2:2089613/91‑1 (MQ=21)
                            ttttttGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGa                                                                >  2:2564342/1‑91 (MQ=255)
                             tttttGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGAt                                                               >  1:3685435/1‑91 (MQ=255)
                             tttttGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGAt                                                               >  1:560501/1‑91 (MQ=255)
                             tttttGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGAt                                                               >  1:2103232/1‑91 (MQ=255)
                             tttttGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGAt                                                               >  1:1652849/1‑91 (MQ=255)
                              ttttGATTATAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGAtt                                                              >  2:737227/1‑91 (MQ=255)
                              ttttGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAACTGAGAAATCTTAAAGATTTTTCAGCTGTCCGAtt                                                              <  2:704272/91‑1 (MQ=18)
                               tttGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGAttt                                                             >  1:739581/1‑91 (MQ=255)
                                 tGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTa                                                           >  1:3222292/1‑91 (MQ=255)
                                 tGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTa                                                           <  2:679432/91‑1 (MQ=255)
                                 tGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTa                                                           >  2:1135939/1‑91 (MQ=255)
                                 tGATTAAAGTTGAATTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTa                                                           >  1:245028/1‑91 (MQ=255)
                                  gATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGAtttt                                                            <  2:1316985/89‑1 (MQ=255)
                                  gATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTaa                                                          <  1:588775/91‑1 (MQ=255)
                                   aTTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTaaa                                                         >  1:1736447/1‑91 (MQ=255)
                                   aTTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTaaa                                                         >  2:1502725/1‑91 (MQ=255)
                                   aTTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTaaa                                                         >  1:1503269/1‑91 (MQ=255)
                                   aTTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTaaa                                                         >  1:1501725/1‑91 (MQ=255)
                                   aTTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTaaa                                                         >  1:1501832/1‑91 (MQ=255)
                                     tAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAata                                                       >  2:2383226/1‑91 (MQ=255)
                                         gTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGc                                                   >  2:3424666/1‑91 (MQ=255)
                                           tGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTa                                                 <  2:3152529/91‑1 (MQ=255)
                                             aTTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTAtt                                               <  2:895646/91‑1 (MQ=255)
                                             aTTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTAtt                                               <  1:2310088/91‑1 (MQ=255)
                                              ttttGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTa                                              >  1:2748299/1‑91 (MQ=255)
                                              ttttGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTa                                              >  1:3883505/1‑91 (MQ=255)
                                              ttttGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTa                                              >  1:1854020/1‑91 (MQ=255)
                                               tttGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAg                                             >  1:1053800/1‑91 (MQ=255)
                                                 tGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGaa                                           >  1:3377865/1‑91 (MQ=255)
                                                   ccGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAAtt                                         >  1:2766960/1‑91 (MQ=255)
                                                   ccGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAAtt                                         <  1:3828848/91‑1 (MQ=255)
                                                   ccGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAAtt                                         <  2:1263584/91‑1 (MQ=255)
                                                      tatgtatTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAg                                      <  2:3014865/91‑1 (MQ=255)
                                                        tgtatTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGtt                                    <  2:709391/91‑1 (MQ=255)
                                                          tatTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGtttt                                  <  1:673893/91‑1 (MQ=255)
                                                           atTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGCTTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTc                                 >  2:2996498/1‑91 (MQ=255)
                                                           atTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTc                                 <  1:205018/91‑1 (MQ=255)
                                                             tAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCgt                               <  1:1777920/91‑1 (MQ=255)
                                                                aTTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCgtagt                            >  1:3556910/1‑91 (MQ=255)
                                                                  tACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTATCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTaa                          >  2:3195540/1‑91 (MQ=255)
                                                                   aCGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAAc                         >  1:129635/1‑91 (MQ=255)
                                                                    cGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAAcc                        <  2:2693192/91‑1 (MQ=255)
                                                                          aaaaTCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAACCTAATTc                  <  1:94027/91‑1 (MQ=255)
                                                                          aaaaTCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAACCTAATTc                  <  1:1233053/91‑1 (MQ=255)
                                                                           aaaTCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAACCTAATTCt                 <  2:359646/91‑1 (MQ=255)
                                                                           aaaTCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAACCTAATTCt                 <  1:2546639/91‑1 (MQ=255)
                                                                           aaaTCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATAGAGCTATTAGAATTTAGTTTTCGTAGTAACCTAATTCt                 <  1:3274184/91‑1 (MQ=255)
                                                                             aTCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAACCTAATTCTaa               <  2:1699694/91‑1 (MQ=255)
                                                                              tCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAACCTAATTCTAAg              >  1:1151760/1‑91 (MQ=255)
                                                                              tCGGACAACTGAGAAATCTTAAAGATTTTTCAGCTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAACCTAATTCTAAg              >  1:1381030/1‑91 (MQ=18)
                                                                                ggACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAACCTAATTCTAAGTc            >  2:3361356/1‑91 (MQ=255)
                                                                                ggACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAACCTAATTCTAAGTc            >  2:3312707/1‑91 (MQ=255)
                                                                                ggACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAACCTAATTCTAAGTc            >  2:764744/1‑91 (MQ=255)
                                                                                    aGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAACCTAATTCTAAGTCtttt        >  1:1239148/1‑91 (MQ=255)
                                                                                    aGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAACCTAATTCTAAGTCtttt        <  1:3846580/91‑1 (MQ=255)
                                                                                     gCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAACCTAATTCTAAGTCTTTTg       <  2:3240430/91‑1 (MQ=255)
                                                                                     acTGAGAAATCTTAAAGATTTTTCAGCTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAACCTAATTCTAAGTCTTTTg       >  2:1404625/2‑91 (MQ=25)
                                                                                      cTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAACCTAATTCTAAGTCTTTTGa      >  1:3419576/1‑91 (MQ=255)
                                                                                       tGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAACCTAATTCTAAGTCTTTTGaa     <  2:1555418/91‑1 (MQ=255)
                                                                                          |                                                                                          
TGTGCTTTAAAGGCATTGTGTATGAGTCTTTTTTGATTAAAGTTGATTTTGCCGTATGTATTAGATTACGTAAAAAAATCGGACAGCTGAAAAATCTTTAAGATTTCTCAGTTGTCCGATTTTAAATATAGCTATTAGAATTTAGTTTTCGTAGTAACCTAATTCTAAGTCTTTTGAAGTT  >  USA300TCH1516_ALE/2744115‑2744295

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: