Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,448,493 C→T 12.4% intergenic (‑252/+25) yxeP_4 ← / ← hutI putative hydrolase YxeP/Imidazolonepropionase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,448,4930CT12.4% 645.9 / 60.8 232intergenic (‑252/+25)yxeP_4/hutIputative hydrolase YxeP/Imidazolonepropionase
Reads supporting (aligned to +/- strand):  ref base C (112/91);  new base T (17/12);  total (129/103)
Fisher's exact test for biased strand distribution p-value = 8.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 239 total aligned reads displayed.

AACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACTTTTTCTGCGTG  >  USA300TCH1516_ALE/2448403‑2448583
                                                                                          |                                                                                          
aaCTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACACAAAGCTCTAAACATTAGTTTAAAcc                                                                                            >  2:2627494/1‑91 (MQ=255)
 aCTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCa                                                                                           <  1:3231882/91‑1 (MQ=255)
  cTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCaa                                                                                          <  2:857499/91‑1 (MQ=255)
     ttATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGc                                                                                       >  2:2705363/1‑91 (MQ=255)
      tataATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCt                                                                                      >  1:1018053/1‑91 (MQ=255)
      tataATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCt                                                                                      >  1:3742845/1‑91 (MQ=255)
      tataATTGAGACTCTAAACAAAAACGTGTTCTCGAATATACTTGTTTTAAAAAAAAAAAAAAAGCTCTAAACATTAGTTTAAACCAATGCt                                                                                      <  2:1386713/91‑1 (MQ=255)
        taATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCCCAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTa                                                                                    >  2:3112482/1‑91 (MQ=255)
        taATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAGCATTGGTTTAAACTAATGTTTa                                                                                    <  2:1179715/91‑1 (MQ=37)
        taATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTa                                                                                    >  1:2800579/1‑91 (MQ=255)
        taATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAGCTCTAAGCATTGGTTTAAACTAATGTTTag                                                                                   <  1:752815/91‑1 (MQ=25)
         aaTTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAGCATTGGTTTAAACTAATGTTTag                                                                                   <  1:3466787/91‑1 (MQ=37)
           ttGAGACTCTATACAAAAGCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTagag                                                                                 <  2:2742661/91‑1 (MQ=255)
           ttGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAGCATTGGTTTAAACTAATGTTTagag                                                                                 <  2:2536027/91‑1 (MQ=25)
           ttGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAGCATTGGTTTAAACTAATGTTTagag                                                                                 <  1:2698464/91‑1 (MQ=25)
           ttGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTagag                                                                                 >  2:2657426/1‑91 (MQ=255)
            tGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGc                                                                                >  1:2703586/1‑91 (MQ=255)
            tGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGc                                                                                >  1:2178809/1‑91 (MQ=255)
            tGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGc                                                                                <  2:3154728/91‑1 (MQ=255)
               gaCTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCttt                                                                             <  1:744297/91‑1 (MQ=255)
               gaCTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCttt                                                                             >  2:1560298/1‑91 (MQ=255)
               aaCTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCttt                                                                             <  2:669887/90‑1 (MQ=255)
                 ctctATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCt                                                                           >  1:1021814/1‑91 (MQ=255)
                  tctATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTa                                                                          >  1:2579014/1‑91 (MQ=255)
                    tataCAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAat                                                                        <  2:623831/91‑1 (MQ=255)
                    tataCAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAat                                                                        <  2:502952/91‑1 (MQ=255)
                    tataCAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAat                                                                        >  1:1360808/1‑91 (MQ=255)
                     ataCAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAatt                                                                       <  2:907276/91‑1 (MQ=255)
                     ataCAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAatt                                                                       >  2:2089075/1‑91 (MQ=255)
                       aCAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAattat                                                                     <  2:1029000/91‑1 (MQ=255)
                       aCAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAattat                                                                     >  2:2586784/1‑91 (MQ=255)
                        cAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAGATCGCTCTAAACATTAGTTTAAACCAATGCTTAGGGCTTTCTAattatt                                                                    <  2:3290270/91‑1 (MQ=255)
                        cAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAattatt                                                                    >  2:77514/1‑91 (MQ=255)
                        cAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAattatt                                                                    >  1:1907576/1‑91 (MQ=255)
                          aaaaCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTAtttt                                                                  <  1:2551150/91‑1 (MQ=255)
                           aaaCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTa                                                                 <  2:3855738/91‑1 (MQ=255)
                            aaCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTAt                                                                <  2:3357173/91‑1 (MQ=255)
                            aaCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTAt                                                                <  1:133920/91‑1 (MQ=255)
                             aCGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATg                                                               <  1:2886083/91‑1 (MQ=255)
                                 gtTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCttt                                                           >  1:1723568/1‑91 (MQ=255)
                                 gtTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCttt                                                           <  1:2699006/91‑1 (MQ=255)
                                     tcGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTaaaa                                                       <  2:792176/91‑1 (MQ=255)
                                          tataCTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATAc                                                  >  1:3006181/1‑91 (MQ=255)
                                            taCTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACtg                                                <  2:2991907/91‑1 (MQ=255)
                                             aCTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACtgt                                               >  1:332420/1‑91 (MQ=255)
                                              cTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACtgtg                                              >  1:1675690/1‑91 (MQ=255)
                                               ttGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTgtgt                                             >  1:2078910/1‑91 (MQ=255)
                                               ttGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTgtgt                                             >  1:879668/1‑91 (MQ=255)
                                                tGTTTTTACAAACCACAAAAAGCTCTAAAGATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGtt                                            <  1:3167457/91‑1 (MQ=255)
                                                tGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGtt                                            >  1:2827612/1‑91 (MQ=255)
                                                tGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGtt                                            >  1:2555579/1‑91 (MQ=255)
                                                tGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGtt                                            >  1:1418504/1‑91 (MQ=255)
                                                 gTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTa                                           <  2:2823447/91‑1 (MQ=255)
                                                       aCAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACg                                     >  1:2821353/1‑91 (MQ=255)
                                                       aCAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACg                                     >  1:1224084/1‑91 (MQ=255)
                                                       aCAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACg                                     <  1:1825333/91‑1 (MQ=255)
                                                        cAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGa                                    >  1:2269533/1‑91 (MQ=255)
                                                           aCCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGa                                 >  2:3771079/1‑91 (MQ=255)
                                                            ccACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTGTTTTATGCTTTAAAAGAGACTGTGTTATCTACGATGAc                                <  2:193091/91‑1 (MQ=255)
                                                             cacaAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGAcc                               <  2:3264350/91‑1 (MQ=255)
                                                               caAAAAGCTCTAAACATTAGTTTAAACCAATGGTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCtt                             <  2:2928620/91‑1 (MQ=255)
                                                                aaaaaGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTa                            <  1:3775018/91‑1 (MQ=255)
                                                                 aaaaGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTAc                           >  1:1432169/1‑91 (MQ=255)
                                                                  aaaGCTCTAAACATTATTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTAcc                          <  2:3036330/91‑1 (MQ=255)
                                                                  aaaGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTAcc                          <  2:3146779/91‑1 (MQ=255)
                                                                  aaaGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTAcc                          >  2:1710499/1‑91 (MQ=255)
                                                                   aaGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCg                         >  2:2083070/1‑91 (MQ=255)
                                                                    aGCTCTAAACATTAGTTTAAACCAATGCTTCGAGCTTTCTAATTATTCTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGt                        >  2:741840/1‑91 (MQ=255)
                                                                     gCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTc                       <  1:3076809/91‑1 (MQ=255)
                                                                     gCTCTAAAATATAGGCTAAGCCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTATCTACGATGACCTTGCCGTCt                      <  2:1887676/91‑1 (MQ=255)
                                                                        ctAAGCATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCttt                    >  2:3944989/1‑91 (MQ=37)
                                                                         taagcATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTta                   >  1:3532248/5‑91 (MQ=37)
                                                                          aaaCATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaa                  >  1:3860625/1‑91 (MQ=255)
                                                                          aaaCATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaa                  <  2:2487078/91‑1 (MQ=255)
                                                                          aaaCATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaa                  <  2:1453980/91‑1 (MQ=255)
                                                                          aaaCATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTATTtaa                  >  2:3647702/1‑91 (MQ=255)
                                                                           agcATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaat                 >  2:2849826/3‑91 (MQ=38)
                                                                           agcATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaat                 >  2:341766/3‑91 (MQ=37)
                                                                           agcATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaat                 >  2:3318389/3‑91 (MQ=38)
                                                                           agcATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAACATACTGTGTTATCTACGATGACCTTACCGTCTTtaat                 >  2:347884/3‑91 (MQ=37)
                                                                           aaCATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaat                 <  1:136019/91‑1 (MQ=255)
                                                                             cATTGGTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTtaataa               >  1:1505791/1‑91 (MQ=38)
                                                                              attggTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAAc              >  2:466938/5‑91 (MQ=38)
                                                                              aTTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAAc              >  2:1253802/1‑91 (MQ=255)
                                                                               ttggTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACt             >  2:3241482/4‑91 (MQ=38)
                                                                                 aGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACttt           >  2:3495115/1‑91 (MQ=255)
                                                                                 aGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACttt           >  1:826753/1‑91 (MQ=255)
                                                                                 aGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACttt           >  1:827011/1‑91 (MQ=255)
                                                                                  gTTTAAACTAATGTTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACtttt          >  2:2860599/1‑91 (MQ=255)
                                                                                   tttAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACttttt         <  1:631051/91‑1 (MQ=255)
                                                                                   tttAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACttttt         <  1:3399882/91‑1 (MQ=255)
                                                                                    ttAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACTTTTTc        >  1:2588541/1‑91 (MQ=255)
                                                                                    ttAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACTTTTTc        <  2:1553142/91‑1 (MQ=255)
                                                                                     tAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACTTTTTCt       <  1:486525/91‑1 (MQ=255)
                                                                                     tAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACTTTTTCt       <  2:3535557/91‑1 (MQ=255)
                                                                                     tAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACTTTTTCt       >  2:1995330/1‑91 (MQ=255)
                                                                                        aCCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACTTTTTCTGCg    <  1:2673737/91‑1 (MQ=255)
                                                                                        aCCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACTTTTTCTGCg    >  2:2394861/1‑91 (MQ=255)
                                                                                         ccAATGCTTAGAGCTTTCTAATTATTTTATGCTTTATAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACTTTTTCTGCGt   <  2:2407014/91‑1 (MQ=255)
                                                                                          cAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGCTATCTACGATGACCTTACCGTCTTTAATAACTTTTTCTGCGtg  >  1:3908187/1‑91 (MQ=255)
                                                                                          |                                                                                          
AACTTTTATAATTGAGACTCTATACAAAAACGTGTTCTCGAATATACTTGTTTTTACAAACCACAAAAAGCTCTAAACATTAGTTTAAACCAATGCTTAGAGCTTTCTAATTATTTTATGCTTTAAAAGATACTGTGTTATCTACGATGACCTTACCGTCTTTAATAACTTTTTCTGCGTG  >  USA300TCH1516_ALE/2448403‑2448583

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