Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,608,829 C→T 7.5% intergenic (+110/+150) USA300TCH1516_02485 → / ← USA300HOU_RS13450 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,608,8290CT 7.5% 577.1 / 22.4 186intergenic (+110/+150)USA300TCH1516_02485/USA300HOU_RS13450hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (77/94);  new base T (11/3);  total (89/97)
Fisher's exact test for biased strand distribution p-value = 2.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 186 total aligned reads displayed.

AACAAGTTGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGGTAGTAATCATTG  >  USA300TCH1516_ALE/2608694‑2608968
                                                                                                                                       |                                                                                                                                           
aaCAAGTTGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGc                                                                                                                                        >  2:236031/1‑141 (MQ=255)
  aaaGTTGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACGTCCTGTATTGGATGTGGCtt                                                                                                                                      <  4:284407/140‑1 (MQ=255)
    aGTTGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCtttt                                                                                                                                    >  1:936893/1‑141 (MQ=255)
    aGTTGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCtttt                                                                                                                                    <  3:380386/141‑1 (MQ=255)
      ttGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAg                                                                                                                                  <  2:1643850/141‑1 (MQ=255)
      ttGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAg                                                                                                                                  <  2:865352/141‑1 (MQ=255)
            aTGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGt                                                                                                                            >  3:618727/1‑141 (MQ=255)
             tGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGtt                                                                                                                           <  2:1520539/141‑1 (MQ=255)
                aGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTGAAAAATAATATAAGTATCACTAAAGCAACATCCAATATAGGACGTCGCTTGTAGTATCGATCCa                                                                                                                        >  2:788869/1‑141 (MQ=255)
                 ggATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAg                                                                                                                       <  2:533487/141‑1 (MQ=255)
                 ggATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAg                                                                                                                       <  4:140282/141‑1 (MQ=255)
                   aTTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGAt                                                                                                                     <  3:447710/141‑1 (MQ=255)
                       tAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTtgt                                                                                                                 <  2:1289733/141‑1 (MQ=255)
                        aaTATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTtgtg                                                                                                                <  4:588508/141‑1 (MQ=255)
                        aaTATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTtgtg                                                                                                                >  3:11140/1‑141 (MQ=255)
                          tatCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTtgtaaa                                                                                                              >  1:1350497/1‑138 (MQ=255)
                           atCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                                                                             >  2:1095363/1‑141 (MQ=255)
                             cATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAg                                                                                                           >  1:644104/1‑141 (MQ=255)
                              aTTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGt                                                                                                          <  2:1587020/141‑1 (MQ=255)
                                  cGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAgg                                                                                                      >  1:73037/1‑141 (MQ=255)
                                  cGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAgg                                                                                                      <  4:346854/141‑1 (MQ=255)
                                   gTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAggg                                                                                                     <  3:647469/141‑1 (MQ=255)
                                   gTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAggg                                                                                                     <  1:720805/141‑1 (MQ=255)
                                   gTTTATATATAGCAAACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAggg                                                                                                     <  1:117865/141‑1 (MQ=255)
                                         atataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACt                                                                                               <  3:426750/141‑1 (MQ=255)
                                          tataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACtt                                                                                              <  2:463448/141‑1 (MQ=255)
                                            taGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaa                                                                                            >  1:329488/1‑141 (MQ=255)
                                             aGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaat                                                                                           <  1:1099544/141‑1 (MQ=255)
                                                aGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATGGTTCCAGATTTGTGGGGAGTTAGGGTTTACttaattca                                                                                        >  2:645822/1‑139 (MQ=255)
                                                  aCATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAg                                                                                      >  1:195129/1‑141 (MQ=255)
                                                       aTAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTc                                                                                 >  3:641913/1‑141 (MQ=255)
                                                         aGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTg                                                                               <  1:80921/141‑1 (MQ=255)
                                                         aGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTg                                                                               >  1:1315967/1‑141 (MQ=255)
                                                          gAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGc                                                                              <  2:58064/141‑1 (MQ=255)
                                                            aTTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGccc                                                                            <  3:72445/141‑1 (MQ=255)
                                                             ttttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGtt                                                                                                                           <  2:1372230/93‑1 (MQ=255)
                                                             ttttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGtt                                                                                                                           >  1:1372230/1‑93 (MQ=255)
                                                             ttttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCt                                                                           <  2:1537970/141‑1 (MQ=255)
                                                             ttttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCt                                                                           <  1:511769/141‑1 (MQ=255)
                                                              tttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTc                                                                          <  1:519918/141‑1 (MQ=255)
                                                               ttATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCa                                                                         <  1:1469393/141‑1 (MQ=255)
                                                                   atGTAAATCTTGTAGGTAAACGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg                                                                     <  1:89796/141‑1 (MQ=255)
                                                                    tGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGc                                                                    >  1:964893/1‑141 (MQ=255)
                                                                      tAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTc                                                                  <  4:417279/141‑1 (MQ=255)
                                                                      tAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCAAATCGAATTTAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTc                                                                  >  1:380420/1‑141 (MQ=255)
                                                                          tCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCa                                                              >  2:322353/1‑141 (MQ=255)
                                                                              gTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGa                                                          >  4:21999/1‑141 (MQ=255)
                                                                              gTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGa                                                          <  2:1229744/141‑1 (MQ=255)
                                                                              gTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGa                                                          >  1:486249/1‑141 (MQ=255)
                                                                              gTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGATTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGa                                                          >  1:852390/1‑141 (MQ=255)
                                                                                  gTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCGATTGCTCACCAACGAGGGc                                                      <  1:714452/141‑1 (MQ=255)
                                                                                     aTCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAg                                                   >  4:242661/1‑141 (MQ=255)
                                                                                     aTCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAg                                                   >  1:234491/1‑141 (MQ=255)
                                                                                       cGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTa                                                                                                        >  1:59983/1‑86 (MQ=255)
                                                                                        gTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTa                                                >  1:1043343/1‑141 (MQ=255)
                                                                                           ttAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACat                                             >  3:579751/1‑141 (MQ=255)
                                                                                             aaaaaTAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTGAAGATTAGTGTGGATGTAAGGTTTTCTTAATTATAGGAGTCTTCCCTCAATTGCTCGCCAACGAGGGCAAGTTACgtca                                           >  1:1653313/1‑141 (MQ=255)
                                                                                             aaaaaTAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACatca                                           >  3:72660/1‑141 (MQ=255)
                                                                                              aaaaTAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACatcaa                                          <  2:83828/141‑1 (MQ=255)
                                                                                                aataatATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTAACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACatcaatc                                        >  1:864936/1‑141 (MQ=255)
                                                                                                     tataAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACt                                   >  1:1006168/1‑141 (MQ=255)
                                                                                                     tataAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACt                                   <  3:486071/141‑1 (MQ=255)
                                                                                                     tataAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACt                                   <  2:763880/141‑1 (MQ=255)
                                                                                                     tataAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACt                                   <  3:31401/141‑1 (MQ=255)
                                                                                                     tataAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACt                                   <  3:267437/141‑1 (MQ=255)
                                                                                                         aGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCaa                               >  3:544448/1‑141 (MQ=255)
                                                                                                         aGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCaa                               >  3:142598/1‑141 (MQ=255)
                                                                                                          gTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAAt                              <  1:1054515/141‑1 (MQ=255)
                                                                                                             tGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTgc                           <  4:582894/141‑1 (MQ=255)
                                                                                                               aCTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGt                                                  <  1:1145074/116‑1 (MQ=255)
                                                                                                               aCTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGt                                                  >  2:1145074/1‑116 (MQ=255)
                                                                                                               aCTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTgccg                         <  1:1268519/141‑1 (MQ=255)
                                                                                                               aCTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTgccg                         <  1:987661/141‑1 (MQ=255)
                                                                                                                 tAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTgccgcc                       >  3:340331/1‑141 (MQ=255)
                                                                                                                  aaaGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGccgccg                      >  1:402717/1‑141 (MQ=255)
                                                                                                                  aaaGCCACATCCAATATAGGACGTGGCTTGTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCGATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGccgccg                      <  4:347964/141‑1 (MQ=255)
                                                                                                                   aaGCCACATCCAATATAGGGCGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGa                     <  2:73037/141‑1 (MQ=255)
                                                                                                                    aGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGAc                    >  2:1139160/1‑141 (MQ=255)
                                                                                                                      ccACGTCCTATATTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCaa                                                             <  4:547813/98‑1 (MQ=38)
                                                                                                                         cgtCCTATATTGGATGGGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAggg               >  2:289541/3‑141 (MQ=255)
                                                                                                                         caTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAggg               >  1:1511258/1‑141 (MQ=255)
                                                                                                                         caTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAggg               >  2:1588546/1‑141 (MQ=255)
                                                                                                                          aTCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAgggg              >  1:1318533/1‑141 (MQ=255)
                                                                                                                           tcctatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCa                                                                 >  2:982579/5‑89 (MQ=255)
                                                                                                                           tcctatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGAGgt             >  1:262862/5‑141 (MQ=255)
                                                                                                                           tCCTATATTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCa                                                                 <  1:982579/89‑1 (MQ=255)
                                                                                                                            cctatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGgta            >  1:323735/4‑141 (MQ=255)
                                                                                                                            ccAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGgta            <  4:693130/141‑1 (MQ=255)
                                                                                                                             ctatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGgtag           >  1:187109/3‑141 (MQ=255)
                                                                                                                             ctatatTGGATGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGgtag           >  2:951189/3‑141 (MQ=255)
                                                                                                                             cAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGgtag           >  2:1027545/1‑141 (MQ=255)
                                                                                                                                 ataGGACGTGGCTTTTAGTATCGTTCCCGATTTGTGGGGCGTTCGGGTTTCCTTAATTAAAGGGGTCTTCCCTCAATTGCTCTCCGCCGCGGGCAAGTTACATCTATCGCACTTATATTGCCGCCGACTAGGGgtagttag       <  2:600960/141‑4 (MQ=255)
                                                                                                                                 ataGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCTATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGGTAGTAAt       <  1:976245/141‑1 (MQ=255)
                                                                                                                                    ggACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGGTAGTAATCAt    >  3:225241/1‑141 (MQ=255)
                                                                                                                                    ggACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGGTAGTAATCAt    >  4:310419/1‑141 (MQ=255)
                                                                                                                                    ggACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGGTAGTAATCAt    >  3:695611/1‑141 (MQ=255)
                                                                                                                                    ggACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGGTAGTAATCAt    <  2:969043/141‑1 (MQ=255)
                                                                                                                                      aCGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGGTAGTAATCATTg  <  4:641913/141‑1 (MQ=255)
                                                                                                                                      aCGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGGTAGTAATCATTg  >  1:291119/1‑141 (MQ=255)
                                                                                                                                      aCGTGGCTTTTAGTATCGTTCCAGAATTGTGGGGAGATAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCCACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGGTAGTAATCATTg  >  2:117990/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                           
AACAAGTTGAAAATGAAGGATTCTAATATCATTTCGTTTATATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACATCAATCACACTTCAATTGCCGCCGACTAGGGGTAGTAATCATTG  >  USA300TCH1516_ALE/2608694‑2608968

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