Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 109,970 T→C 7.5% intergenic (+49/‑172) 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,9700TC 7.5% 709.6 / 33.9 226intergenic (+49/‑172)plc/USA300HOU_RS005151‑phosphatidylinositol phosphodiesterase/putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (125/84);  new base C (9/8);  total (134/92)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 238 total aligned reads displayed.

TATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTGATAAT  >  USA300TCH1516_ALE/109880‑110059
                                                                                           |                                                                                         
tatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAat                                                                                           >  1:3044767/1‑91 (MQ=255)
tatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAat                                                                                           <  2:1354152/91‑1 (MQ=255)
 atgatgATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAatt                                                                                          <  2:3135759/91‑1 (MQ=255)
  tgatgaTATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAatta                                                                                         >  1:1476541/1‑91 (MQ=255)
   gatgatATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAattat                                                                                        >  1:2252206/1‑91 (MQ=255)
       atatCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTc                                                                                    >  1:2185703/1‑91 (MQ=255)
       atatCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTc                                                                                    >  1:2185604/1‑91 (MQ=255)
       atatCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTc                                                                                    <  1:140665/91‑1 (MQ=255)
       atatCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTAGTTTc                                                                                    >  1:3233958/1‑91 (MQ=255)
        tatCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCa                                                                                   >  2:410923/1‑91 (MQ=255)
         atCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTAGTGAAAATAATTAGAATCTAACTATTTTCa                                                                                   <  1:563921/91‑1 (MQ=25)
          tCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTAGTGAAAATAATTAGAATCTAACTATTTTCat                                                                                  <  2:2019299/91‑2 (MQ=25)
          tCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCact                                                                                 <  2:3323073/91‑1 (MQ=255)
          tCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCact                                                                                 >  1:1464661/1‑91 (MQ=255)
          tCGTAAGTGAAATTATAGATAGTAATAAATAAGGATGCAATAATGATAGTAAGACGAGTATGAAAATAATTAGAATCTAACTATTTTCata                                                                                 <  2:2707098/91‑3 (MQ=37)
           cGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTGAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCacta                                                                                >  2:3417765/1‑91 (MQ=255)
            gTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCactac                                                                               >  1:1918363/1‑91 (MQ=255)
                gTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGt                                                                           <  2:1073188/91‑1 (MQ=255)
                gTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGt                                                                           <  1:2239187/91‑1 (MQ=255)
                 tGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtt                                                                          <  2:354977/91‑1 (MQ=255)
                 tGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtt                                                                          >  1:2827769/1‑91 (MQ=255)
                  gAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTAGTGAAAATAATTAGAATCTAACTATTTTCATACTCGtct                                                                         <  2:3687486/91‑3 (MQ=17)
                  gAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGttt                                                                         <  1:3483792/91‑1 (MQ=255)
                   aaaTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtttt                                                                        <  2:3519628/91‑1 (MQ=255)
                   aaaTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtttt                                                                        >  2:3439882/1‑91 (MQ=255)
                   aaaTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtttt                                                                        <  2:551104/91‑1 (MQ=255)
                   aaaTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtttt                                                                        <  1:3442096/91‑1 (MQ=255)
                    aaTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGttttt                                                                       <  1:1044432/91‑1 (MQ=255)
                     aTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTa                                                                      >  2:400597/1‑91 (MQ=255)
                        atagatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTAttt                                                                   >  1:3687795/1‑91 (MQ=255)
                        atagatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTAttt                                                                   >  1:3170881/1‑91 (MQ=255)
                         tagatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTAtttt                                                                  <  1:3466911/91‑1 (MQ=255)
                         tagatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTAtttt                                                                  <  2:1577444/91‑1 (MQ=255)
                         tagatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTAtttt                                                                  <  1:1143804/91‑1 (MQ=255)
                           gatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGa                                                                >  1:537071/1‑91 (MQ=255)
                            atagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGaa                                                               >  1:1366190/1‑91 (MQ=255)
                              agTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGaaaa                                                             >  1:3370872/1‑91 (MQ=255)
                               gTAATAAATAAGGATTCAATAATGATATTAAGACGGGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAAt                                                            >  2:1275124/1‑91 (MQ=255)
                                 aataaataAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATaa                                                          <  2:2345717/91‑1 (MQ=255)
                                  ataaataaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAg                                                         <  2:3381791/91‑1 (MQ=255)
                                    aaataaGGATTCAATAATGCTATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGta                                                       >  2:1228587/1‑91 (MQ=255)
                                     aataaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGtaa                                                      >  2:3128958/1‑91 (MQ=255)
                                     aataaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGtaa                                                      <  1:2645579/91‑1 (MQ=255)
                                       taaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGtaata                                                    >  1:2471739/1‑91 (MQ=255)
                                        aaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGtaataa                                                   >  2:3069551/1‑91 (MQ=255)
                                        aaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGtaataa                                                   <  2:1464661/91‑1 (MQ=255)
                                         aGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTaataat                                                  >  1:2404548/1‑91 (MQ=255)
                                          ggATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAAtt                                                 <  1:1299028/91‑1 (MQ=255)
                                           gATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTc                                                >  1:1730196/1‑91 (MQ=255)
                                            aTTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTca                                               >  1:202786/1‑91 (MQ=255)
                                              tCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTcaac                                             >  2:2933603/1‑91 (MQ=255)
                                                  taatGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAata                                         >  2:813087/1‑91 (MQ=255)
                                                   aatGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAatat                                        >  1:2088495/1‑91 (MQ=255)
                                                    atGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTTAATTTAAGTAATAATTCAACAAtatt                                       >  2:1978982/1‑91 (MQ=255)
                                                    atGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatt                                       >  1:2261654/1‑91 (MQ=255)
                                                     tGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatta                                      <  1:3593016/91‑1 (MQ=255)
                                                     tGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatta                                      >  2:2456411/1‑91 (MQ=255)
                                                         atTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGATTATAAGTAATAATTCAACAATATTATaaa                                  >  2:3790266/1‑91 (MQ=255)
                                                         atTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATaaa                                  >  1:207422/1‑91 (MQ=255)
                                                         atTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATaaa                                  >  1:2952339/1‑91 (MQ=255)
                                                         atTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATaaa                                  >  2:2844137/1‑91 (MQ=255)
                                                         atTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATaaa                                  >  2:2749666/1‑91 (MQ=255)
                                                          ttAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAAt                                 >  2:2782290/1‑91 (MQ=255)
                                                           tAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAAtt                                >  2:2078086/1‑91 (MQ=255)
                                                             aGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGa                              >  1:938853/1‑91 (MQ=255)
                                                             aGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGa                              >  1:1872787/1‑91 (MQ=255)
                                                             aGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGa                              <  1:2895565/91‑1 (MQ=255)
                                                              gACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGaa                             >  2:2874366/1‑91 (MQ=255)
                                                               aCGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAAc                            <  2:2505083/91‑1 (MQ=255)
                                                               aCGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAACAATTCAACAATATTATAAATTGAAc                            >  1:3617342/1‑91 (MQ=255)
                                                                cGAGTAGTGAAAATAATTAGAATCTAACTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAAc                            >  1:649346/1‑91 (MQ=25)
                                                                cGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACa                           >  1:840579/1‑91 (MQ=255)
                                                                cGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACa                           <  1:3750243/91‑1 (MQ=255)
                                                                 gAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAg                          <  2:2287159/91‑1 (MQ=255)
                                                                  aGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGa                         >  1:2907329/1‑91 (MQ=255)
                                                                   gTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTc                                                >  1:957360/1‑67 (MQ=255)
                                                                   gTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAt                        >  1:1753041/1‑91 (MQ=255)
                                                                     agtGAAAATAATTAGAATCTAACTATTTTCATACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttg                      >  2:3789574/3‑91 (MQ=25)
                                                                     aTGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttg                      >  1:2598789/1‑91 (MQ=255)
                                                                     aTGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttg                      >  2:319548/1‑91 (MQ=255)
                                                                     aTGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTCGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttg                      >  1:3078476/1‑91 (MQ=255)
                                                                        gAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttgtt                    >  1:3796188/1‑91 (MQ=255)
                                                                         aaaaTAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttgttt                   >  1:1168716/1‑91 (MQ=255)
                                                                         aaaaTAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttgttt                   <  2:2179898/91‑1 (MQ=255)
                                                                         aaaaTAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttgttt                   <  2:1976037/91‑1 (MQ=255)
                                                                           aaTAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttgtttgt                 <  2:2088495/91‑1 (MQ=255)
                                                                            aTAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTtgtg                <  2:3447600/91‑1 (MQ=255)
                                                                               gTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGaaa             >  2:3008841/1‑91 (MQ=255)
                                                                                 tAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAAtt           >  1:2344912/1‑91 (MQ=255)
                                                                                  aGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAAttt          >  2:1120660/1‑91 (MQ=255)
                                                                                  aGATTCTAATTATTTTCACTACTCGTCTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAAttt          >  2:2208011/1‑91 (MQ=255)
                                                                                   gATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAAtttt         <  1:3531521/91‑1 (MQ=255)
                                                                                   gATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAAtttt         <  2:2629958/91‑1 (MQ=255)
                                                                                     ttCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTg       <  1:405224/91‑1 (MQ=255)
                                                                                     ttCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTg       >  1:2150211/1‑91 (MQ=255)
                                                                                     ttCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTg       <  1:714307/91‑1 (MQ=255)
                                                                                      tCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTGa      <  2:2347246/91‑1 (MQ=255)
                                                                                       cTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTGat     <  2:1822958/91‑1 (MQ=255)
                                                                                       cTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTGat     >  2:625003/1‑91 (MQ=255)
                                                                                        tAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTGata    <  1:2675289/91‑1 (MQ=255)
                                                                                           |                                                                                         
TATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAATTTTTGATAAT  >  USA300TCH1516_ALE/109880‑110059

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: