Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,243,499 T→A 12.4% intergenic (+77/+32) USA300HOU_RS11475 → / ← pyrG hypothetical protein/CTP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,243,4990TA12.4% 730.2 / 71.1 243intergenic (+77/+32)USA300HOU_RS11475/pyrGhypothetical protein/CTP synthase
Reads supporting (aligned to +/- strand):  ref base T (104/106);  new base A (16/14);  total (122/121)
Fisher's exact test for biased strand distribution p-value = 8.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 250 total aligned reads displayed.

ACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCG  >  USA300TCH1516_ALE/2243409‑2243589
                                                                                          |                                                                                          
 cTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTAGTTTCAATGTTTa                                                                                           <  1:2972620/91‑1 (MQ=255)
 cTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTa                                                                                           <  2:1101032/91‑1 (MQ=255)
 cTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTa                                                                                           <  1:1239109/91‑1 (MQ=255)
 cTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTa                                                                                           <  1:264602/91‑1 (MQ=255)
  ttAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTaa                                                                                          >  2:129135/1‑91 (MQ=255)
     gAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACat                                                                                       <  2:63490/91‑2 (MQ=255)
     gAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACat                                                                                       <  2:2680963/91‑2 (MQ=255)
     gAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACat                                                                                       <  1:216998/91‑2 (MQ=255)
     gAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACaa                                                                                       <  1:1027359/91‑1 (MQ=255)
     gAATTATAAAAGTCACTATTGTATTTAATACCTTGTCGCTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACatt                                                                                      <  2:3639693/91‑3 (MQ=37)
       aTTATAAAGGGCACTATTGTTATTTCATACCTTGTCACTAACATAATATTTAAATTATACGAAAAACCGGTACTTTCATTGTTAAACattg                                                                                     <  2:3594438/91‑4 (MQ=37)
        ttATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACATTGa                                                                                    <  1:12082/91‑1 (MQ=38)
          ataAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGaaa                                                                                  >  1:2367667/1‑91 (MQ=255)
           taAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAg                                                                                 >  2:2719352/1‑91 (MQ=255)
           taAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAg                                                                                 <  2:3198939/91‑1 (MQ=255)
            aaaaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGt                                                                                <  2:465454/91‑1 (MQ=255)
             aaaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTa                                                                               >  1:2335825/1‑91 (MQ=255)
              aaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAc                                                                              <  1:641389/91‑1 (MQ=255)
               aGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAcc                                                                             <  2:1986963/91‑1 (MQ=255)
               aGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAcc                                                                             <  1:2835994/91‑1 (MQ=255)
                gTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCg                                                                            <  1:2735119/91‑1 (MQ=255)
                gTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCg                                                                            <  2:3071485/91‑1 (MQ=255)
                 tCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCgg                                                                           <  1:2627975/91‑1 (MQ=255)
                   aCTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGtt                                                                         <  2:3784937/91‑1 (MQ=255)
                    cTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGttt                                                                        <  1:1366801/91‑1 (MQ=255)
                     tATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGtttt                                                                       <  2:3539218/91‑1 (MQ=255)
                     tATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGtttt                                                                       <  1:187813/91‑1 (MQ=255)
                         gTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAACATTGAAAGTACCGGtttttcgt                                                                   <  1:2783867/91‑5 (MQ=21)
                         gTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTttatta                                                                   <  1:1304691/91‑1 (MQ=255)
                          ttATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAg                                                                  <  1:1610172/91‑1 (MQ=255)
                          ttATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAg                                                                  >  1:223274/1‑91 (MQ=255)
                           tATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGc                                                                 <  2:1660045/91‑1 (MQ=255)
                              ttCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCaaa                                                              >  1:1277638/1‑91 (MQ=255)
                                cATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatt                                                            >  2:544314/1‑91 (MQ=255)
                                   aCCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttat                                                         <  1:384433/91‑1 (MQ=255)
                                   aCCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTCCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttat                                                         >  1:3875263/1‑91 (MQ=255)
                                    ccTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttatt                                                        <  1:394731/91‑1 (MQ=255)
                                     cTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttattt                                                       >  2:3176990/1‑91 (MQ=255)
                                       tGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttatttat                                                     <  2:140337/91‑1 (MQ=255)
                                           aCTAACATAATCTTAAAATTATATGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGGACCGGTTTTATTAGCAAATTTATTTATTttg                                                 <  2:2957119/91‑1 (MQ=255)
                                           aCTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttg                                                 >  1:2772794/1‑91 (MQ=255)
                                           aCTAACATAATATTAAAAATATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAGATTTATTTATTttg                                                 <  2:1057477/91‑1 (MQ=255)
                                             tAACATAATATTAAAATTATACTAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAGTTTATTTATTttgtt                                               <  1:2104087/91‑1 (MQ=255)
                                              aaCATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgttg                                              <  1:23124/91‑1 (MQ=255)
                                               aCATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGa                                             <  2:2237076/91‑1 (MQ=255)
                                               aCATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGa                                             >  1:880495/1‑91 (MQ=255)
                                               aCATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGa                                             >  1:2000152/1‑91 (MQ=255)
                                                    atatTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAttt                                        >  2:267503/1‑91 (MQ=255)
                                                    atatTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTACTTATTTTGTTGATAttt                                        >  2:1796292/1‑91 (MQ=255)
                                                     tatTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAtttt                                       >  1:795682/1‑91 (MQ=255)
                                                     tatTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAtttt                                       <  1:2548195/91‑1 (MQ=255)
                                                     tatTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAtttt                                       >  2:60143/1‑91 (MQ=255)
                                                      atTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTa                                      >  1:2210217/1‑91 (MQ=255)
                                                      atTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTa                                      >  2:811701/1‑91 (MQ=255)
                                                        tAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAAt                                    <  1:3498355/91‑1 (MQ=255)
                                                        tAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAAt                                    >  1:2423780/1‑91 (MQ=255)
                                                         aaaaTTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATg                                   >  2:151896/1‑91 (MQ=255)
                                                          aaaTTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGa                                  >  2:2413215/1‑91 (MQ=255)
                                                            aTTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAg                                >  1:614461/1‑91 (MQ=255)
                                                            aTTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAg                                >  2:3220928/1‑91 (MQ=255)
                                                                taCGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTc                            >  2:21576/1‑91 (MQ=255)
                                                                taCGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTc                            >  1:3332155/1‑91 (MQ=255)
                                                                taCGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTc                            <  2:2255515/91‑1 (MQ=255)
                                                                taCGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTc                            >  1:78995/1‑91 (MQ=255)
                                                                  cGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCaa                          >  2:2784125/1‑91 (MQ=255)
                                                                  cGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCaa                          >  1:65896/1‑91 (MQ=255)
                                                                   gAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAAt                         >  1:1789531/1‑91 (MQ=255)
                                                                   gAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAAt                         >  1:2262889/1‑91 (MQ=255)
                                                                   gAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAAt                         >  1:336939/1‑91 (MQ=255)
                                                                    aaaaaCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATa                        >  2:1892727/1‑91 (MQ=255)
                                                                      aaaCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATaaa                      >  1:1700177/1‑91 (MQ=255)
                                                                        aCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATg                    <  1:320394/91‑1 (MQ=255)
                                                                         ccGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGa                   <  1:1580992/91‑1 (MQ=255)
                                                                         ccGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGa                   >  1:1898129/1‑91 (MQ=255)
                                                                          cGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAt                  <  1:2160077/91‑1 (MQ=255)
                                                                           ggTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAtt                 <  1:60143/91‑1 (MQ=255)
                                                                            gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                >  1:2511302/1‑91 (MQ=255)
                                                                            gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                >  1:2913486/1‑91 (MQ=255)
                                                                            gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                <  2:3864660/91‑1 (MQ=255)
                                                                            gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                <  1:668892/91‑1 (MQ=255)
                                                                            gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                <  2:659314/91‑1 (MQ=255)
                                                                               cTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaaa             >  1:1093464/1‑91 (MQ=255)
                                                                               cTTTCAACGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaaa             <  1:1166393/91‑1 (MQ=255)
                                                                                tttCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaaaa            <  1:2257721/91‑1 (MQ=255)
                                                                                 ttCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaaaaa           >  2:3432793/1‑91 (MQ=255)
                                                                                 ttCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaaaaa           >  2:782992/1‑91 (MQ=255)
                                                                                  tCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAAt          <  2:1222272/91‑1 (MQ=255)
                                                                                  tCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAAt          >  2:2051355/1‑91 (MQ=255)
                                                                                   cAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATc         >  1:365788/1‑91 (MQ=255)
                                                                                    attGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCg        >  1:2536147/3‑91 (MQ=38)
                                                                                    attGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCg        >  1:492989/3‑91 (MQ=38)
                                                                                     aTGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCgg       >  2:2521426/1‑91 (MQ=255)
                                                                                     aTGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCgg       >  2:2314418/1‑91 (MQ=255)
                                                                                      tGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCggg      >  2:2531310/1‑91 (MQ=255)
                                                                                       gTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAACATCGGGt     >  1:1936794/1‑91 (MQ=255)
                                                                                       gTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGt     >  2:3674305/1‑91 (MQ=255)
                                                                                        tttAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTg    <  2:3277819/91‑1 (MQ=255)
                                                                                         ttAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGc   >  1:1817040/1‑91 (MQ=255)
                                                                                          tAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCg  >  1:3709947/1‑91 (MQ=255)
                                                                                          tAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCg  >  1:959065/1‑91 (MQ=255)
                                                                                          |                                                                                          
ACTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCG  >  USA300TCH1516_ALE/2243409‑2243589

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: