Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,243,499 T→A 7.0% 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,4990TA 7.0% 1224.4 / 63.6 334intergenic (+77/+32)USA300HOU_RS11475/pyrGhypothetical protein/CTP synthase
Reads supporting (aligned to +/- strand):  ref base T (161/147);  new base A (13/10);  total (175/159)
Fisher's exact test for biased strand distribution p-value = 8.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

Only 100 of 351 total aligned reads displayed.

CTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCG  >  USA300TCH1516_ALE/2243410‑2243589
                                                                                         |                                                                                          
  tAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCATTGTTAAAc                                                                                         <  1:591031/91‑1 (MQ=255)
    gAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACaa                                                                                       <  1:4368212/91‑1 (MQ=255)
      aTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATg                                                                                     >  2:884269/1‑91 (MQ=255)
      aTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATg                                                                                     <  1:119708/91‑1 (MQ=255)
      aTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATg                                                                                     <  1:4615441/91‑1 (MQ=255)
      aTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATg                                                                                     <  2:523919/91‑1 (MQ=255)
      aTTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAAATATACGAAAAACCGGTACTTTCAATGTTTAACAATg                                                                                     <  2:4076672/91‑1 (MQ=255)
       ttATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGa                                                                                    <  1:4045614/91‑1 (MQ=255)
        tataAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGaa                                                                                   >  1:1157385/1‑91 (MQ=255)
        tataAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGaa                                                                                   >  1:4561698/1‑91 (MQ=255)
          taAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAg                                                                                 >  2:4896533/1‑91 (MQ=255)
          taAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAg                                                                                 <  2:2583834/91‑1 (MQ=255)
          taAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAg                                                                                 >  2:3604404/1‑91 (MQ=255)
          taAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAg                                                                                 <  2:5390716/91‑1 (MQ=255)
           aaaaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGt                                                                                <  2:2352414/91‑1 (MQ=255)
            aaaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTa                                                                               >  1:5246754/1‑91 (MQ=255)
            aaaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTa                                                                               >  2:4857100/1‑91 (MQ=255)
            aaaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTa                                                                               >  1:4170125/1‑91 (MQ=255)
             aaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAATATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAc                                                                              <  2:5442422/91‑1 (MQ=255)
             aaGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAc                                                                              <  2:755514/91‑1 (MQ=255)
              aGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAcc                                                                             <  1:1109334/91‑1 (MQ=255)
              aGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAcc                                                                             <  2:4613701/91‑1 (MQ=255)
              aGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAcc                                                                             >  2:1747440/1‑91 (MQ=255)
              aGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTAcc                                                                             <  2:4757168/91‑1 (MQ=255)
              aGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATATTTAACAATGACAGTAcc                                                                             >  1:4713915/1‑91 (MQ=255)
               gTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCg                                                                            <  2:4331377/91‑1 (MQ=255)
                tCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCgg                                                                           <  2:3043711/91‑1 (MQ=255)
                   cgaTTGTTATTTCATACCTTGTCACTAACATAATGTTAAAATTATACGAAAAACCGGTACTTTCAATGTTGAACAATGAAAGTACCGGctt                                                                        >  2:4711400/3‑91 (MQ=38)
                   cTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGttt                                                                        <  1:3258522/91‑1 (MQ=255)
                    tATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGtttt                                                                       >  1:4123888/1‑91 (MQ=255)
                    tATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGtttt                                                                       >  1:2338332/1‑91 (MQ=255)
                     aTTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTtta                                                                      >  2:4343290/1‑91 (MQ=255)
                     aTTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTtta                                                                      <  1:2779619/91‑1 (MQ=255)
                      ttGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTttat                                                                     >  2:2250755/1‑91 (MQ=255)
                       tGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTttatt                                                                    <  1:2851739/91‑1 (MQ=255)
                           aTTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCa                                                                >  2:1658404/1‑91 (MQ=255)
                            tttCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCaa                                                               >  1:3118922/1‑91 (MQ=255)
                             ttCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCaaa                                                              >  1:1299437/1‑91 (MQ=255)
                               cATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatt                                                            >  1:3803682/1‑91 (MQ=255)
                                   ccTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttatt                                                        >  2:2564890/1‑91 (MQ=255)
                                     ttGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGTAAGTACCGGTTTTATTAGCAAatttattta                                                      >  2:1318862/1‑91 (MQ=255)
                                     ttGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttattta                                                      >  1:4461318/1‑91 (MQ=255)
                                     ttGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttattta                                                      <  1:1322454/91‑1 (MQ=255)
                                       gTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttatttatt                                                    >  2:3292340/1‑91 (MQ=255)
                                        tCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAatttatttattt                                                   >  2:3250022/1‑91 (MQ=255)
                                           cTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgt                                                >  2:3985102/1‑91 (MQ=255)
                                            tAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgtt                                               <  2:1199582/91‑1 (MQ=255)
                                             aaCATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgttg                                              <  2:4935242/91‑1 (MQ=255)
                                             aaCATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgttg                                              >  2:3653937/1‑91 (MQ=255)
                                             aaCATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTttgttg                                              <  1:3981070/91‑1 (MQ=255)
                                              aCATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGa                                             <  2:3876033/91‑1 (MQ=255)
                                              aCATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGa                                             >  2:562113/1‑91 (MQ=255)
                                              aCATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGa                                             >  2:1391696/1‑91 (MQ=255)
                                                   atatTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATAttt                                        >  2:2028711/1‑91 (MQ=255)
                                                     atTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTa                                      >  2:699453/1‑91 (MQ=255)
                                                     atTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTAGATTGTTGATATTTTa                                      >  1:4333758/1‑91 (MQ=255)
                                                     agtAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATGAGCAAATTTATTTATTTTGTTGATATTTTa                                      <  2:4490053/89‑1 (MQ=255)
                                                      ttAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTaa                                     >  1:4861186/1‑91 (MQ=255)
                                                         aaaTTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGa                                  <  1:4519495/91‑1 (MQ=255)
                                                           aTTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAg                                >  1:5044899/1‑91 (MQ=255)
                                                            ttATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGc                               <  2:1795010/91‑1 (MQ=255)
                                                            ttATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGc                               <  2:1041746/91‑1 (MQ=255)
                                                            ttATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGc                               <  2:3377756/91‑1 (MQ=255)
                                                            ttATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGc                               <  2:5111108/91‑1 (MQ=255)
                                                             tataCGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCt                              <  1:5321508/91‑1 (MQ=255)
                                                                aCGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCa                           >  1:3910811/1‑91 (MQ=255)
                                                                 cGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCaa                          <  2:5334317/91‑1 (MQ=255)
                                                                  gAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAAt                         >  1:3192979/1‑91 (MQ=255)
                                                                  gAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAAt                         <  2:2947068/91‑1 (MQ=255)
                                                                    aaaaCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATaa                       >  2:2774337/1‑91 (MQ=255)
                                                                      aaCCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAAt                     <  1:5152963/91‑1 (MQ=255)
                                                                        ccGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGa                   >  1:3518296/1‑91 (MQ=255)
                                                                        ccGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGa                   <  1:861698/91‑1 (MQ=255)
                                                                        ccGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGa                   <  2:3749273/91‑1 (MQ=255)
                                                                         cGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAt                  >  1:3491451/1‑91 (MQ=255)
                                                                         cGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAt                  <  1:2089361/91‑1 (MQ=255)
                                                                          ggTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAtt                 >  1:3037432/1‑91 (MQ=255)
                                                                          ggTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAtt                 >  2:4435379/1‑91 (MQ=255)
                                                                           gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                <  2:5059532/91‑1 (MQ=255)
                                                                           gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                <  2:2649433/91‑1 (MQ=255)
                                                                           gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                >  1:4818993/1‑91 (MQ=255)
                                                                           gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                >  1:1514759/1‑91 (MQ=255)
                                                                           gTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGAttt                >  1:4878780/1‑91 (MQ=255)
                                                                            tACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTa               >  1:1020438/1‑91 (MQ=255)
                                                                            tACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTa               >  1:3007288/1‑91 (MQ=255)
                                                                            tACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTa               >  1:492810/1‑91 (MQ=255)
                                                                             aCTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGCTATTTTAATGAAGCTTCAATAAATGATTTaa              <  1:5523001/91‑1 (MQ=255)
                                                                              cTTTCATTGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaat             >  2:1810/1‑90 (MQ=37)
                                                                                ttCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTaaaaa           >  1:2026076/1‑91 (MQ=255)
                                                                                ttCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATGTTAATGAAGCTTCAATATATGATTTataaa           >  2:1980862/1‑91 (MQ=255)
                                                                                 tcattGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAAt          >  1:25995/5‑91 (MQ=38)
                                                                                   attGTTAAACATTGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCg        >  1:1809872/3‑91 (MQ=38)
                                                                                   aaTGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTACTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCg        <  1:4249263/91‑1 (MQ=255)
                                                                                    aTGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCgg       >  2:633213/1‑91 (MQ=255)
                                                                                       tttAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTg    <  1:2121489/91‑1 (MQ=255)
                                                                                       tttAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTg    <  1:1643065/91‑1 (MQ=255)
                                                                                         tAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTACAAATCGGGTGCg  >  2:3731535/1‑91 (MQ=255)
                                                                                         tAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCg  <  2:3932781/91‑1 (MQ=255)
                                                                                         tAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCg  <  2:3522171/91‑1 (MQ=255)
                                                                                         tAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCg  >  1:2728193/1‑91 (MQ=255)
                                                                                         |                                                                                          
CTTAGAATTATAAAAGTCACTATTGTTATTTCATACCTTGTCACTAACATAATATTAAAATTATACGAAAAACCGGTACTTTCAATGTTTAACAATGAAAGTACCGGTTTTATTAGCAAATTTATTTATTTTGTTGATATTTTAATGAAGCTTCAATAAATGATTTAAAAATCGGGTGCG  >  USA300TCH1516_ALE/2243410‑2243589

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: