Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,840,841 T→A 5.4% intergenic (‑162/+143) hisZ ← / ← USA300HOU_RS14555 ATP phosphoribosyltransferase regulatory subunit/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,840,8410TA 5.4% 751.3 / 19.2 225intergenic (‑162/+143)hisZ/USA300HOU_RS14555ATP phosphoribosyltransferase regulatory subunit/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (122/90);  new base A (7/6);  total (129/96)
Fisher's exact test for biased strand distribution p-value = 7.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 244 total aligned reads displayed.

CCACGCTAAAGTGTAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGTACTTCTTGCACATATTGT  >  USA300TCH1516_ALE/2840751‑2840931
                                                                                          |                                                                                          
ccACGCTAAAGTGTAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAAt                                                                                            <  1:1706514/91‑1 (MQ=255)
ccACGCTAAAGTGTAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAAt                                                                                            <  2:1634009/91‑1 (MQ=255)
 cACGCTAAAGTGTAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAAtt                                                                                           <  1:3166914/91‑1 (MQ=255)
 cACGCTAAAGTGTAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAAtt                                                                                           <  2:616828/91‑1 (MQ=255)
 cACGCTAAAGTGTAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAAtt                                                                                           <  2:616501/91‑1 (MQ=255)
  aCGCTAAAGTGTAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAAttt                                                                                          >  2:137076/1‑91 (MQ=255)
  aCGCTAAAGTGTAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAAttt                                                                                          >  2:1455782/1‑91 (MQ=255)
  aCGCTAAAGTGTAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAAttt                                                                                          <  1:2582533/91‑1 (MQ=255)
  aCGCTAAAGTGTAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAAttt                                                                                          >  2:2020509/1‑91 (MQ=255)
    gCTAAAGTGTAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAAttttt                                                                                        >  2:2087936/1‑91 (MQ=255)
      tAAAGTGTAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTta                                                                                      >  1:3764579/1‑91 (MQ=255)
       aaaGTGTAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTtaa                                                                                     >  2:3161826/1‑91 (MQ=255)
       aaaGTGTAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTtaa                                                                                     <  2:2301970/91‑1 (MQ=255)
        aaGTGTAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTtaaa                                                                                    >  2:2840473/1‑91 (MQ=255)
         aGTGTAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTtaaaa                                                                                   >  1:1048264/1‑91 (MQ=255)
         aGTGTAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTtaaaa                                                                                   >  1:3057909/1‑91 (MQ=255)
         aGTGTAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAAAATTTtaaaa                                                                                   <  2:2197572/91‑1 (MQ=255)
          gtgtAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTtaaaat                                                                                  <  1:2325742/91‑1 (MQ=255)
           tgtAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTtaaaata                                                                                 <  2:415954/91‑1 (MQ=255)
            gtAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTtaaaataa                                                                                <  2:3002166/91‑1 (MQ=255)
            gtAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTtaaaataa                                                                                <  2:124907/91‑1 (MQ=255)
                 atGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAg                                                                           >  2:1437064/1‑91 (MQ=255)
                 atGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAAATCGGTGTAGGGATTGTAATAATTGTTGAAGTAATAAAg                                                                           <  2:721836/91‑1 (MQ=255)
                         taaCTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACa                                                                   <  1:3813375/91‑1 (MQ=255)
                         taaCTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACa                                                                   >  1:1575477/1‑91 (MQ=255)
                            cTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACAttt                                                                >  2:1776567/1‑91 (MQ=255)
                             ttATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTtg                                                               >  2:778518/1‑91 (MQ=255)
                              tATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTtgt                                                              >  1:2932139/1‑91 (MQ=255)
                                 cATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGcc                                                           <  1:381686/91‑1 (MQ=255)
                                 cATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGcc                                                           >  1:2107644/1‑91 (MQ=255)
                                  aTTTTTAGCAATACATTACCACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCt                                                          >  2:3396833/1‑91 (MQ=255)
                                  aTTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCt                                                          <  2:391696/91‑1 (MQ=255)
                                  aTTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCt                                                          >  2:303407/1‑91 (MQ=255)
                                   tttttAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCta                                                         <  1:2955733/91‑1 (MQ=255)
                                     tttAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCtata                                                       >  1:141704/1‑91 (MQ=255)
                                     tttAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTAAAAATAAAAAAGGCAAGACATTTGTGCCtata                                                       <  2:1209734/91‑1 (MQ=255)
                                     tttAGCAATACATTACAACCTTTTTCAGAAAATTCGGAGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCtata                                                       >  1:3273722/1‑91 (MQ=255)
                                       tAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATaaa                                                     <  1:1797176/91‑1 (MQ=255)
                                         gCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATaaaaa                                                   >  1:38713/1‑91 (MQ=255)
                                            aTACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAATATGCAAGACATTTGTGCCTATAAAAATGc                                                >  2:1685968/1‑91 (MQ=255)
                                            aTACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGc                                                >  2:2436548/1‑91 (MQ=255)
                                            aTACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGc                                                >  2:2574314/1‑91 (MQ=255)
                                              aCATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCtt                                              >  2:1400096/1‑91 (MQ=255)
                                               cATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTa                                             <  2:695252/91‑1 (MQ=255)
                                                aTTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTaa                                            <  2:2863284/91‑1 (MQ=255)
                                                  tACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAAcc                                          <  1:365243/91‑1 (MQ=255)
                                                  tACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAAcc                                          <  2:3189894/91‑1 (MQ=255)
                                                  tACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAAcc                                          >  1:3633099/1‑91 (MQ=255)
                                                   aCAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAAc                                           >  2:3574853/1‑89 (MQ=255)
                                                   aCAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAAc                                           <  1:3574853/89‑1 (MQ=255)
                                                    cAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCaa                                        >  2:156830/1‑91 (MQ=255)
                                                    cAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCaa                                        >  2:1666848/1‑91 (MQ=255)
                                                      aCCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAAcc                                          <  1:3491181/87‑1 (MQ=255)
                                                      aCCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAAcc                                          >  2:3491181/1‑87 (MQ=255)
                                                        cTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGAtt                                    <  1:2247009/91‑1 (MQ=255)
                                                         tttttCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGAttt                                   >  1:794755/1‑91 (MQ=255)
                                                          ttttCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGAtttt                                  <  1:2524881/91‑1 (MQ=255)
                                                           tttCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGAttttt                                 >  1:808748/1‑91 (MQ=255)
                                                            ttCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTta                                >  1:3240972/1‑91 (MQ=255)
                                                             tCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTtat                               <  1:3381220/91‑1 (MQ=255)
                                                             tCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTtat                               >  2:2953938/1‑91 (MQ=255)
                                                              cAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTtata                              <  2:2494306/91‑1 (MQ=255)
                                                               aGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTatat                             >  1:260363/1‑91 (MQ=255)
                                                               aGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTatat                             >  2:1093833/1‑91 (MQ=255)
                                                               aGAAAATTCGGTGTATTGATTTTAAAAATATATAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTatat                             <  2:3737300/91‑1 (MQ=255)
                                                                gAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATAtt                            <  1:2751135/91‑1 (MQ=255)
                                                                 aaaaTTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTg                           <  2:3773943/91‑1 (MQ=255)
                                                                 aaaaTTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTg                           >  1:2362825/1‑91 (MQ=255)
                                                                 aaaaTTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTg                           >  1:962242/1‑91 (MQ=255)
                                                                 aaaaTTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTg                           >  1:1383800/1‑91 (MQ=255)
                                                                  aaaTTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGa                          <  2:956629/91‑1 (MQ=255)
                                                                  aaaTTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGa                          >  2:2270593/1‑91 (MQ=255)
                                                                     ttCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGt                       >  1:1257516/1‑91 (MQ=255)
                                                                     ttCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGt                       >  1:1815228/1‑91 (MQ=255)
                                                                      tCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGtt                      <  1:2901543/91‑1 (MQ=255)
                                                                      tCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGtt                      <  2:2787587/91‑1 (MQ=255)
                                                                      tCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGtt                      <  2:2130255/91‑1 (MQ=255)
                                                                       tggTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAATAATGCTTAACCAAGATTTTTATATTGAAGTTg                     <  2:3808629/90‑1 (MQ=255)
                                                                       cGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTg                     >  1:2115890/1‑91 (MQ=255)
                                                                       cGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTg                     >  1:1833798/1‑91 (MQ=255)
                                                                        ggTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGt                    >  2:2026724/1‑91 (MQ=255)
                                                                        ggTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGt                    >  2:2425941/1‑91 (MQ=255)
                                                                         gtgtATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGTa                   >  1:1325800/1‑91 (MQ=255)
                                                                         gtgtATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGTa                   <  2:1169963/91‑1 (MQ=255)
                                                                         gtgtATTGACTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCTTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTCTa                   >  2:2622924/1‑91 (MQ=255)
                                                                           gtATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGTAct                 <  1:1964665/91‑1 (MQ=255)
                                                                            tATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGTActt                >  2:2732961/1‑91 (MQ=255)
                                                                             aTTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGTActtc               <  2:2380985/91‑1 (MQ=255)
                                                                               tGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGTActtctt             >  1:2415924/1‑91 (MQ=255)
                                                                                gATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGTACTTCTTg            >  1:249063/1‑91 (MQ=255)
                                                                                 aTTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGTACTTCTTGc           >  2:1198484/1‑91 (MQ=255)
                                                                                  ttttAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGTACTTCTTGca          <  1:2974991/91‑1 (MQ=255)
                                                                                   tttAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGTACTTCTTGcac         <  2:308127/91‑1 (MQ=255)
                                                                                   tttAAAAAATTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGTACTTCTTGcac         >  1:3754718/1‑91 (MQ=255)
                                                                                   tttAAAAAATTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGTACTTCTTGcac         >  1:654247/1‑91 (MQ=255)
                                                                                     tAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGTACTTCTTGCACat       >  2:1315109/1‑91 (MQ=255)
                                                                                     tAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGTACTTCTTGCACat       >  2:2349749/1‑91 (MQ=255)
                                                                                      aaaaTTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGTACTTCTTGCACata      >  1:1517120/1‑91 (MQ=255)
                                                                                         aTTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGTACTTCTTGCACATATTg   >  1:2222168/1‑91 (MQ=255)
                                                                                          ttttttAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGTACTTCTTGCACATATTGt  >  1:543401/1‑91 (MQ=255)
                                                                                          |                                                                                          
CCACGCTAAAGTGTAATATGATTAATAACTTATCATTTTTAGCAATACATTACAACCTTTTTCAGAAAATTCGGTGTATTGATTTTAAAATTTTTTAAAATAAAAAAGGCAAGACATTTGTGCCTATAAAAATGCTTAACCAAGATTTTTATATTGAAGTTGTACTTCTTGCACATATTGT  >  USA300TCH1516_ALE/2840751‑2840931

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: