Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,107,417 G→A 9.1% intergenic (+52/‑116) pdhD → / → USA300HOU_RS05475 Dihydrolipoyl dehydrogenase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,107,4170GA 9.1% 661.7 / 40.8 219intergenic (+52/‑116)pdhD/USA300HOU_RS05475Dihydrolipoyl dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (103/95);  new base A (15/5);  total (118/101)
Fisher's exact test for biased strand distribution p-value = 6.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

Only 100 of 230 total aligned reads displayed.

GCAGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCGATGCTAACAA  >  USA300TCH1516_ALE/1107327‑1107507
                                                                                           |                                                                                          
gcagAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACg                                                                                            >  1:3564615/1‑91 (MQ=255)
 cagAAAAAGCTATCGGATACCCGATCCATACAATGTAATAACTGATTATGTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACgg                                                                                           <  2:2949526/91‑1 (MQ=255)
 cagAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACgg                                                                                           >  1:84653/1‑91 (MQ=255)
 cagAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACgg                                                                                           <  1:1553468/91‑1 (MQ=255)
 cagAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACgg                                                                                           <  2:3179561/91‑1 (MQ=255)
  agAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGc                                                                                          <  1:1529764/91‑1 (MQ=255)
  agAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGc                                                                                          <  1:2432252/91‑1 (MQ=255)
  agAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGc                                                                                          <  2:2958269/91‑1 (MQ=255)
   gAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAATGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCt                                                                                         <  1:3266716/91‑1 (MQ=255)
   gAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCt                                                                                         <  1:585266/91‑1 (MQ=255)
   gAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCt                                                                                         <  1:1331437/91‑1 (MQ=255)
   gAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCt                                                                                         <  2:2187790/91‑1 (MQ=255)
    aaaaaGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCtt                                                                                        <  1:639534/91‑1 (MQ=255)
    aaaaaGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCCGTAGCATATGCTACAGCtt                                                                                        <  2:18394/91‑1 (MQ=255)
     aaaaGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCttt                                                                                       >  2:1652580/1‑91 (MQ=255)
       aaGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCttttt                                                                                     <  1:39112/91‑1 (MQ=255)
        aGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCtttttt                                                                                    <  2:2301468/91‑1 (MQ=255)
         gCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTg                                                                                   >  2:1265967/1‑91 (MQ=255)
         gCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTg                                                                                   >  1:1983797/1‑91 (MQ=255)
          cTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGt                                                                                  >  1:3461702/1‑91 (MQ=255)
          cTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGt                                                                                  >  1:3459295/1‑91 (MQ=255)
            aTCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGttt                                                                                >  1:3806301/1‑91 (MQ=255)
            aTCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGttt                                                                                >  2:544870/1‑91 (MQ=255)
             tCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGtttt                                                                               >  2:1052034/1‑91 (MQ=255)
             tCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGtttt                                                                               >  1:1211153/1‑91 (MQ=255)
             tCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGtttt                                                                               >  2:2138085/1‑91 (MQ=255)
             tCGGATACCCAAGCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGtttt                                                                               <  2:2649176/91‑1 (MQ=255)
               ggATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAg                                                                             >  1:3256150/1‑91 (MQ=255)
                  tACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTa                                                                          <  1:3673102/91‑1 (MQ=255)
                    cccAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTaaa                                                                        >  2:2465623/1‑91 (MQ=255)
                     ccAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAg                                                                       >  2:2922007/1‑91 (MQ=255)
                     ccAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTCTTTGTTTTAGGTAAAg                                                                       >  2:3715922/1‑91 (MQ=255)
                        aTCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTaa                                                                    >  2:1527267/1‑91 (MQ=255)
                         tCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAAt                                                                   >  2:3595812/1‑91 (MQ=255)
                          ccATACAATGTAATAACTGATTATCTATAAAGATTCAGTTATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATg                                                                  <  2:2415594/91‑1 (MQ=255)
                          ccATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATg                                                                  >  1:689306/1‑91 (MQ=255)
                          ccATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATg                                                                  >  1:1067276/1‑91 (MQ=255)
                             tACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTaa                                                               >  2:2847618/1‑91 (MQ=255)
                              aCAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAGTGTaa                                                               >  1:2164662/1‑91 (MQ=255)
                              aCAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAg                                                              >  1:2099276/1‑91 (MQ=255)
                              aCAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAg                                                              >  1:1811822/1‑91 (MQ=255)
                              aCAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAg                                                              >  2:3106865/1‑91 (MQ=255)
                               cAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAgg                                                             <  1:3396721/91‑1 (MQ=255)
                                aaTGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGa                                                            <  2:3481977/91‑1 (MQ=255)
                                   gTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAAt                                                         <  1:1612828/91‑1 (MQ=255)
                                     aataaCTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTg                                                       <  2:3571610/91‑1 (MQ=255)
                                        aaCTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGAtt                                                    <  2:927600/91‑1 (MQ=255)
                                         aCTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGAttt                                                   >  2:2807007/1‑91 (MQ=255)
                                         aCTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGAttt                                                   <  1:985642/91‑1 (MQ=255)
                                            gATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTgag                                                >  2:2043024/1‑91 (MQ=255)
                                               tATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGata                                             <  1:1197447/91‑1 (MQ=255)
                                               tATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGata                                             >  1:2617836/1‑91 (MQ=255)
                                                aTCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGatat                                            >  1:2494250/1‑91 (MQ=255)
                                                 tCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTTGGTAAAGTAATGTAAGGAAATTGATTTGAGATATc                                           >  2:2871070/1‑91 (MQ=255)
                                                 tCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATc                                           <  1:3772984/91‑1 (MQ=255)
                                                    ataAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGa                                                            <  2:1514722/71‑1 (MQ=255)
                                                    ataAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGtt                                        >  2:2096177/1‑91 (MQ=255)
                                                    ataAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGtt                                        <  2:1771462/91‑1 (MQ=255)
                                                     taAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTa                                       >  1:296448/1‑91 (MQ=255)
                                                      aaaGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTaa                                      >  2:638469/1‑91 (MQ=255)
                                                      aaaGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTaa                                      <  2:3256150/91‑1 (MQ=255)
                                                      aaaGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTaa                                      >  1:802557/1‑91 (MQ=255)
                                                      aaaGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTCGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTaa                                      >  1:1131591/1‑91 (MQ=255)
                                                         gATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACAt                                   >  2:1651692/1‑91 (MQ=255)
                                                          aTTCAGTCATTAAAAGCTGTAGCATATGCTACTGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACAtg                                  <  1:3299711/91‑1 (MQ=255)
                                                          aTTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACAtg                                  <  2:3103344/91‑1 (MQ=255)
                                                          aTTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACAtg                                  <  2:1106420/91‑1 (MQ=255)
                                                           ttCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGaaa                                                          <  1:2147666/66‑1 (MQ=255)
                                                           ttCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACAtgt                                 <  2:2758778/91‑1 (MQ=255)
                                                            tCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACAtgtg                                >  1:481934/1‑91 (MQ=255)
                                                            tCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACAtgtg                                >  2:1084991/1‑91 (MQ=255)
                                                             cAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGa                               >  1:2683576/1‑91 (MQ=255)
                                                               gTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGAca                             <  1:467018/91‑1 (MQ=255)
                                                                      aaaaGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGAcatgcatg                       <  1:128028/91‑1 (MQ=255)
                                                                      aaaaGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGAcatgcatg                       >  1:3630561/1‑91 (MQ=255)
                                                                      aaaaGCCGTAGCATATGCTACAGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGAcatgcatg                       >  1:1416801/1‑91 (MQ=255)
                                                                       aaaGCCGTAGCATATGCTACAGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGt                      >  1:507107/1‑91 (MQ=255)
                                                                         aGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTta                    <  1:222487/91‑1 (MQ=255)
                                                                          gCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTtat                   <  1:2082275/91‑1 (MQ=255)
                                                                           cTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTtata                  >  1:2980759/1‑91 (MQ=255)
                                                                           cTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTtata                  >  2:3211495/1‑91 (MQ=255)
                                                                           cTGTAGCATATGCTACGGCTTTTTTGTTTTAAGTAAAGTAACGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATACATGTtata                  >  2:3184650/1‑91 (MQ=255)
                                                                            tGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGTCATGCATGTTATAc                 <  2:3292389/91‑1 (MQ=255)
                                                                             gTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGGTCTCGTTAACATGTGACATGCATGTTATACt                >  1:3105761/1‑91 (MQ=255)
                                                                             gTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACt                >  1:1872397/1‑91 (MQ=255)
                                                                             gTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACt                >  1:225874/1‑91 (MQ=255)
                                                                             gTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACt                <  1:2847618/91‑1 (MQ=255)
                                                                             gTAGCATATGCTACAGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACt                >  2:710507/1‑91 (MQ=255)
                                                                             gTAGCATATGCTACAGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACt                >  2:2227637/1‑91 (MQ=255)
                                                                             gTAGCATATGCTACAGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACt                >  1:858250/1‑91 (MQ=255)
                                                                              tAGCATATGCTACAGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTa               >  1:37965/1‑91 (MQ=255)
                                                                               aGCATATGCTACAGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAg              >  1:3351746/1‑91 (MQ=255)
                                                                                gCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTGACATGTGACAGGCATGTAATACTAGc             >  2:3370292/1‑91 (MQ=255)
                                                                                gCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGc             <  1:1916046/91‑1 (MQ=255)
                                                                                 cATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCg            >  1:2030364/1‑91 (MQ=255)
                                                                                    atGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCGATg         <  2:699757/91‑1 (MQ=255)
                                                                                     tGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCGATGc        <  2:377764/91‑1 (MQ=255)
                                                                                       cTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCGATGCTa      >  1:2239396/1‑91 (MQ=255)
                                                                                         aCGGCTTTTTTGTTTTAGGTAAAGTAATGTAATGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCGATGCTAAc    <  2:351427/91‑1 (MQ=255)
                                                                                           ggCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCGATGCTAACaa  <  2:1422810/91‑1 (MQ=255)
                                                                                           |                                                                                          
GCAGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCGATGCTAACAA  >  USA300TCH1516_ALE/1107327‑1107507

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: