Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 109,964 T→A 7.2% intergenic (+43/‑178) plc → / → USA300HOU_RS00515 1‑phosphatidylinositol phosphodiesterase/putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE109,9640TA 7.2% 815.2 / 38.1 245intergenic (+43/‑178)plc/USA300HOU_RS005151‑phosphatidylinositol phosphodiesterase/putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (139/88);  new base A (9/9);  total (148/97)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

Only 100 of 252 total aligned reads displayed.

CCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTG  >  USA300TCH1516_ALE/109874‑110054
                                                                                           |                                                                                          
ccTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGAt                                                                                            <  1:1019218/91‑1 (MQ=255)
   ggATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCt                                                                                         >  1:889946/1‑91 (MQ=255)
   ggATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCt                                                                                         <  2:903297/91‑1 (MQ=255)
   ggATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCt                                                                                         >  2:3474287/1‑91 (MQ=255)
   ggATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCt                                                                                         <  1:2240531/91‑1 (MQ=255)
    gATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTa                                                                                        >  2:1049376/1‑91 (MQ=255)
     atatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTaa                                                                                       >  1:1272728/1‑91 (MQ=255)
     atatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTaa                                                                                       >  2:3759735/1‑91 (MQ=255)
     atatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTaa                                                                                       <  1:505676/91‑1 (MQ=255)
     atatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTaa                                                                                       >  1:984693/1‑91 (MQ=255)
     atatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTaa                                                                                       <  2:1475743/91‑1 (MQ=255)
      tatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAat                                                                                      >  2:2928537/1‑91 (MQ=255)
       atgatgATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAatt                                                                                     <  2:3135759/91‑1 (MQ=255)
        tgatgaTATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAatta                                                                                    >  2:3125838/1‑91 (MQ=255)
        tgatgaTATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAatta                                                                                    <  2:2479685/91‑1 (MQ=255)
        tgatgaTATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAatta                                                                                    >  1:1476541/1‑91 (MQ=255)
         gatgatATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAattat                                                                                   >  1:2251700/1‑91 (MQ=255)
             atatCGTAAGTGAAATTATAGATAGTAATAAATAAGGCTTCAATAATGATAATAAGACGAGTATGGAAATAGTTAGATTCTAATTATTTTc                                                                               >  1:3611363/1‑91 (MQ=255)
             atatCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTc                                                                               <  1:140665/91‑1 (MQ=255)
             atatCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTc                                                                               >  1:2185604/1‑91 (MQ=255)
             atatCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTc                                                                               >  1:2185703/1‑91 (MQ=255)
             atatCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTAGTTTc                                                                               >  1:3233958/1‑91 (MQ=255)
               atCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTAGTGAAAATAATTAGAATCTAACTATTTTCa                                                                              <  1:563804/91‑1 (MQ=25)
                tCGTCAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCact                                                                            >  1:3060188/1‑91 (MQ=255)
                tCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTAGTGAAAATAATTAGAATCTAACTATTTTCat                                                                             <  2:2019299/91‑2 (MQ=25)
                tCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCact                                                                            >  1:1464661/1‑91 (MQ=255)
                tCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCact                                                                            >  2:2086757/1‑91 (MQ=255)
                tCGTAAGTGAAATTATAGATAGTAATAAATAAGGATGCAATAATGATAGTAAGACGAGTATGAAAATAATTAGAATCTAACTATTTTCata                                                                            <  2:2707098/91‑3 (MQ=37)
                 cGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCacta                                                                           <  2:584030/91‑1 (MQ=255)
                  gTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCactac                                                                          >  1:1918363/1‑91 (MQ=255)
                     aGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCg                                                                       >  1:2430708/1‑91 (MQ=255)
                     aGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCg                                                                       <  1:3026992/91‑1 (MQ=255)
                       tGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTAGTGAAAATAATTAGAATCTAACTATTTTCATACTCGtc                                                                     <  1:2886124/91‑2 (MQ=18)
                       tGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtt                                                                     <  2:354977/91‑1 (MQ=255)
                        gAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTAGTGAAAATAATTAGAATCTAACTATTTTCATACTCGtct                                                                    <  2:3687486/91‑3 (MQ=17)
                        gAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGttt                                                                    <  1:1289615/91‑1 (MQ=255)
                         aaaTTATAGATAGTAATAAATAAGGATTCAATATTGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtttt                                                                   <  1:3442348/91‑1 (MQ=255)
                         aaaTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtttt                                                                   >  2:3439924/1‑91 (MQ=255)
                         aaaTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtttt                                                                   <  2:551104/91‑1 (MQ=255)
                          aaTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGttttt                                                                  <  2:1099650/91‑1 (MQ=255)
                          aaTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGttttt                                                                  >  1:787583/1‑91 (MQ=255)
                           aTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTa                                                                 >  1:3510075/1‑91 (MQ=255)
                              atagatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTAttt                                                              >  1:3687795/1‑91 (MQ=255)
                              atagatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTAttt                                                              >  1:3170837/1‑91 (MQ=255)
                               tagatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTAtttt                                                             <  2:1577444/91‑1 (MQ=255)
                               tagatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTAtttt                                                             >  1:343353/1‑91 (MQ=255)
                                  atagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGaa                                                          >  1:1366190/1‑91 (MQ=255)
                                       aataaataAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATaa                                                     <  2:2268597/91‑1 (MQ=255)
                                       aataaataAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATaa                                                     <  2:2345717/91‑1 (MQ=255)
                                        ataaataaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAg                                                    <  2:2339483/91‑1 (MQ=255)
                                           aataaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGtaa                                                 <  1:2645579/91‑1 (MQ=255)
                                            ataaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGtaat                                                >  1:2732356/1‑91 (MQ=255)
                                             taaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGtaata                                               <  2:1222115/91‑1 (MQ=255)
                                              aaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGtaataa                                              >  2:3069551/1‑91 (MQ=255)
                                               aGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTaataat                                             >  2:2027631/1‑91 (MQ=255)
                                                 gATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTc                                           >  1:1730196/1‑91 (MQ=255)
                                                  aTTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTca                                          >  1:202786/1‑91 (MQ=255)
                                                    tCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTcaac                                        >  2:2934064/1‑91 (MQ=255)
                                                       ataatGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAat                                     <  2:2834566/91‑1 (MQ=255)
                                                       ataatGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAat                                     >  1:2325943/1‑91 (MQ=255)
                                                         aatGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAatat                                   >  1:2088495/1‑91 (MQ=255)
                                                          atGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatt                                  >  2:1292758/1‑91 (MQ=255)
                                                          atGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatt                                  >  2:1976804/1‑91 (MQ=255)
                                                          atGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatt                                  >  2:406091/1‑91 (MQ=255)
                                                          atGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatt                                  >  1:2261654/1‑91 (MQ=255)
                                                           tGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAATAAtatta                                 >  2:3437344/1‑91 (MQ=255)
                                                           tGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatta                                 >  2:2456411/1‑91 (MQ=255)
                                                           tGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatta                                 <  1:10469/91‑1 (MQ=255)
                                                           tGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatta                                 <  1:3470479/91‑1 (MQ=255)
                                                               atTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATaaa                             >  1:3592554/1‑91 (MQ=255)
                                                               atTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATaaa                             >  2:2844137/1‑91 (MQ=255)
                                                               atTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATaaa                             >  1:207422/1‑91 (MQ=255)
                                                               atTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATaaa                             >  2:2749666/1‑91 (MQ=255)
                                                                ttAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAAt                            >  1:3714006/1‑91 (MQ=255)
                                                                ttAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAAt                            >  2:3879592/1‑91 (MQ=255)
                                                                  aaGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTg                          >  2:2134357/1‑91 (MQ=255)
                                                                   aGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGa                         >  1:938853/1‑91 (MQ=255)
                                                                   aGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGa                         >  1:1873257/1‑91 (MQ=255)
                                                                    gACGAGTATGAAAATAGTTAGATTCTAATCATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTACAAATTGaa                        >  2:3379483/1‑91 (MQ=255)
                                                                     aCGAGTAGTGAAAATAATTAGAATCTAACTATTTTCATACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAAc                       >  1:2557722/1‑91 (MQ=18)
                                                                     aCGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAAc                       <  2:3510075/91‑1 (MQ=255)
                                                                      cGAGTAGTGAAAATAATTAGAATCTAACTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAAc                       >  1:649012/1‑91 (MQ=25)
                                                                       gAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAg                     <  2:2189661/91‑1 (MQ=255)
                                                                        aGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGa                    >  2:1381102/1‑91 (MQ=255)
                                                                         gTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTc                                           >  1:957360/1‑67 (MQ=255)
                                                                         gTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAt                   >  1:1753041/1‑91 (MQ=255)
                                                                           agtGAAAATAATTAGAATCTAACTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAtt                  >  1:3344496/3‑91 (MQ=37)
                                                                           agtGAAAATAATTAGAATCTAACTATTTTCATACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttg                 >  2:3789586/3‑91 (MQ=25)
                                                                           aTGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttg                 >  2:319548/1‑91 (MQ=255)
                                                                              gAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttgtt               >  1:3796188/1‑91 (MQ=255)
                                                                               aaaaTAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttgttt              <  2:194533/91‑1 (MQ=255)
                                                                                aaaTAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttgtttg             >  2:1250468/1‑91 (MQ=255)
                                                                                 aaTAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttgtttgt            <  1:2928537/91‑1 (MQ=255)
                                                                                 aaTAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttgtttgt            >  1:28980/1‑91 (MQ=255)
                                                                                      ttAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAAt       >  2:3394178/1‑91 (MQ=255)
                                                                                      ttAGAATCTAACTATTTTCACTACTCGTTTTTATTTTGTAAATAAGTAATAATTCAACAATATTATAAATTGTACAGTTTGTTTGTGAAAt       >  2:2640056/1‑91 (MQ=255)
                                                                                        aGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAAttt     >  1:2396314/1‑91 (MQ=255)
                                                                                        aGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAAGAATTCAACAATATTATAAATTGAACAGATGGTTTGTGAAAttt     >  1:3877018/1‑91 (MQ=255)
                                                                                           ttCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTg  <  1:714191/91‑1 (MQ=255)
                                                                                           ttCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTg  >  1:2150211/1‑91 (MQ=255)
                                                                                           |                                                                                          
CCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTG  >  USA300TCH1516_ALE/109874‑110054

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: