Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 973,213 C→A 6.6% E175* (GAG→TAG)  yitU_1 ← 5‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase YitU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE973,2130CA 6.6% 784.0 / 32.6 226E175* (GAG→TAG) yitU_15‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase YitU
Reads supporting (aligned to +/- strand):  ref base C (101/110);  new base A (9/6);  total (110/116)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 230 total aligned reads displayed.

ATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTC  >  USA300TCH1516_ALE/973123‑973303
                                                                                          |                                                                                          
aTTGTGTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTc                                                                                            <  2:3137646/91‑1 (MQ=255)
aTTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTc                                                                                            <  2:342426/91‑1 (MQ=255)
 ttGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCa                                                                                           <  2:1868468/91‑1 (MQ=255)
   gTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAat                                                                                         >  1:1317690/1‑91 (MQ=255)
   gTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAat                                                                                         >  1:440775/1‑91 (MQ=255)
    tCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAata                                                                                        <  2:3639677/91‑1 (MQ=255)
           tGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCgg                                                                                 <  2:1212806/91‑1 (MQ=255)
           tGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCgg                                                                                 <  2:160436/91‑1 (MQ=255)
           tGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCgg                                                                                 <  2:3072816/91‑1 (MQ=255)
            gCTCAATGCCTCTTGCTTTATTAATACCAAGTTCTACCATTTCAATCACAGGGAATGGTCCACCCCAGCGTCGATGCTCAATATGAtagat                                                                                >  1:3101603/1‑87 (MQ=255)
              tCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCAt                                                                              <  2:92008/91‑1 (MQ=255)
               cAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATa                                                                             >  2:1818197/1‑91 (MQ=255)
               cAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATa                                                                             >  1:198721/1‑91 (MQ=255)
               cAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATa                                                                             >  2:1709241/1‑91 (MQ=255)
                   gCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATaaaaa                                                                         >  1:3693178/1‑91 (MQ=255)
                    ccTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAAt                                                                        <  1:1844805/91‑1 (MQ=255)
                    ccTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATAACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAAt                                                                        <  2:10981/91‑1 (MQ=255)
                      tctTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGa                                                                      >  2:2861407/1‑91 (MQ=255)
                         tGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTa                                                                   >  2:3327594/1‑91 (MQ=255)
                         tGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTa                                                                   >  2:3250596/1‑91 (MQ=255)
                          gCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTaa                                                                  <  2:125225/91‑1 (MQ=255)
                          gCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTaa                                                                  >  1:1164317/1‑91 (MQ=255)
                           cTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAg                                                                 <  2:806071/91‑1 (MQ=255)
                              tattaATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCat                                                              >  2:818904/1‑91 (MQ=255)
                                 taATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATAtt                                                           <  1:3430622/91‑1 (MQ=255)
                                 taATACCAAGTTTTACAATGTCAATGACAGGGAATGGTGCGCCCCAGCGTCGACGCTCAATATGATCGGCATAAAAATGAGGAAGCATAtt                                                           >  1:2220138/1‑91 (MQ=255)
                                  aaTACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATAttt                                                          <  1:2861407/91‑1 (MQ=255)
                                     aCCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTg                                                       <  1:183621/91‑1 (MQ=255)
                                      cctaGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGa                                                      <  1:1890441/88‑1 (MQ=255)
                                       cAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGAt                                                     >  2:2557704/1‑91 (MQ=255)
                                       cAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGAt                                                     <  1:2968468/91‑1 (MQ=255)
                                        aaGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGAtt                                                    >  1:205051/1‑91 (MQ=255)
                                         aGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGAttt                                                   >  1:631063/1‑91 (MQ=255)
                                         aGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGAttt                                                   >  1:2128588/1‑91 (MQ=255)
                                         aGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGAttt                                                   <  1:2669253/91‑1 (MQ=255)
                                           ttttACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCa                                                 >  1:3387280/1‑91 (MQ=255)
                                             ttACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAgg                                               <  1:1818197/91‑1 (MQ=255)
                                              tACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGt                                              <  2:203217/91‑1 (MQ=255)
                                              tACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGt                                              >  1:3492058/1‑91 (MQ=255)
                                               aCAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTa                                             >  2:2398867/1‑91 (MQ=255)
                                               aCAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTa                                             >  1:1064227/1‑91 (MQ=255)
                                                  aTTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTAttt                                          <  2:1884905/91‑1 (MQ=255)
                                                  aTTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTAttt                                          <  2:475614/91‑1 (MQ=255)
                                                  aTTTCAATGACAGGGAATGGTGCGCCCCAGCGACGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTAttt                                          <  1:3299161/91‑1 (MQ=255)
                                                    ttCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTa                                        >  2:2909713/1‑91 (MQ=255)
                                                     tCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTAc                                       <  1:2573648/91‑1 (MQ=255)
                                                        aTGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACttt                                    <  2:1188151/91‑1 (MQ=255)
                                                        aTGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACttt                                    <  2:2633904/91‑1 (MQ=255)
                                                        aTGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACttt                                    <  1:3250596/91‑1 (MQ=255)
                                                         tGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTttc                                   >  2:513405/1‑91 (MQ=255)
                                                          gACAGGGAATGGTGCGCCCCAGCGTCGTTGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTttct                                  <  2:70262/91‑1 (MQ=255)
                                                          gACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTttct                                  <  1:445715/91‑1 (MQ=255)
                                                          gACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTttct                                  >  2:2078082/1‑91 (MQ=255)
                                                           aCAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTttctt                                 <  2:1352547/91‑1 (MQ=255)
                                                            cAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTttcttc                                <  1:142794/91‑1 (MQ=255)
                                                               ggAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCggg                             >  2:3499682/1‑90 (MQ=255)
                                                               ggAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGc                             <  1:931243/91‑1 (MQ=255)
                                                                gAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCt                            >  1:2370646/1‑91 (MQ=255)
                                                                gAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCt                            >  2:3025704/1‑91 (MQ=255)
                                                                gAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCt                            >  2:3823330/1‑91 (MQ=255)
                                                                gAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCt                            >  2:1492794/1‑91 (MQ=255)
                                                                  aTGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTc                          >  1:210274/1‑91 (MQ=255)
                                                                   tGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCa                         >  2:653808/1‑91 (MQ=255)
                                                                   tGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCa                         <  2:2615716/91‑1 (MQ=255)
                                                                   tGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCa                         <  2:2003025/91‑1 (MQ=255)
                                                                     gTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAt                       >  1:939346/1‑91 (MQ=255)
                                                                      tGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAtt                      <  2:940681/91‑1 (MQ=255)
                                                                      tGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAtt                      >  1:963805/1‑91 (MQ=255)
                                                                      tGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAtt                      <  1:1615323/91‑1 (MQ=255)
                                                                      tGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAtt                      <  1:479874/91‑1 (MQ=255)
                                                                      tGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAtt                      <  1:690349/91‑1 (MQ=255)
                                                                       gcgcCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTa                     <  1:594731/91‑1 (MQ=255)
                                                                         gcCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTaaa                   <  1:1984423/91‑1 (MQ=255)
                                                                         gcCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTaaa                   >  1:3552366/1‑91 (MQ=255)
                                                                          ccccAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGc                                                                >  1:319985/1‑45 (MQ=255)
                                                                          ccccAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTaaaa                  >  2:2421860/1‑91 (MQ=255)
                                                                          ccccAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTaaaa                  >  1:2482266/1‑91 (MQ=255)
                                                                           cccAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAAt                 >  1:610461/1‑91 (MQ=255)
                                                                            ccAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAAtt                <  2:667574/91‑1 (MQ=255)
                                                                             cAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTg               <  2:2504248/91‑1 (MQ=255)
                                                                             cAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTg               >  1:2307201/1‑91 (MQ=255)
                                                                              aGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGa              <  2:1152848/91‑1 (MQ=255)
                                                                                  tCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTa          >  1:2614470/1‑91 (MQ=255)
                                                                                   cGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAg         <  1:192442/91‑1 (MQ=255)
                                                                                   cGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAg         >  1:3011880/1‑91 (MQ=255)
                                                                                    gATGCTCAATATGATCGGCATAAAAATGAGTAATCATATTTTTGATTTCAGGTATTTTAGTTTCTTCGGCTTCAATTAAAATTGAGGTAgg        <  1:3441715/91‑1 (MQ=255)
                                                                                    gATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAgg        >  2:2275151/1‑91 (MQ=255)
                                                                                    gATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAgg        <  2:3558887/91‑1 (MQ=255)
                                                                                    gATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAgg        <  1:352662/91‑1 (MQ=255)
                                                                                    gATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAgg        <  2:2371943/91‑1 (MQ=255)
                                                                                    gATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCCACGGCTTCAACTAAAATTGAGGTAgg        >  2:3362177/1‑91 (MQ=255)
                                                                                    gATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTCCAGGTATTTTACTTTCTTCGGCTTCAATTCAAATTGAGGTAgg        >  1:3873851/1‑91 (MQ=255)
                                                                                     aTGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAggg       >  1:1009053/1‑91 (MQ=255)
                                                                                     aTGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAggg       <  2:1033180/91‑1 (MQ=255)
                                                                                     aTGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAggg       <  2:95453/91‑1 (MQ=255)
                                                                                     aTGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAggg       <  1:1415796/91‑1 (MQ=255)
                                                                                       gCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGa     >  2:3642212/1‑91 (MQ=255)
                                                                                       gCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGa     <  1:3686155/91‑1 (MQ=255)
                                                                                         tCAATATGATCGGCATAAAAATGTCTAAGCATATTTTTGATTTTAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTCGGGGAtt   <  2:3615517/91‑1 (MQ=255)
                                                                                          cAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTc  <  1:378369/91‑1 (MQ=255)
                                                                                          |                                                                                          
ATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTC  >  USA300TCH1516_ALE/973123‑973303

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: