Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,797,772 G→C R7G (CGC→GGC)  phoR ← Alkaline phosphatase synthesis sensor protein PhoR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,797,7720GC95.3% 294.4 / 6.0 86R7G (CGC→GGC) phoRAlkaline phosphatase synthesis sensor protein PhoR
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base C (42/40);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

AAATATATGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACATCAACAATTCT  >  USA300TCH1516_ALE/1797637‑1797909
                                                                                                                                       |                                                                                                                                         
aaataTATGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCcgtggt                                                                                                                                      <  2:26584/141‑1 (MQ=255)
 aataTATGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCgtggtg                                                                                                                                     >  2:125218/1‑141 (MQ=255)
  atatatGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGa                                                                                                                                    >  4:112698/1‑141 (MQ=255)
     tatGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAAc                                                                                                                                 <  1:424773/141‑1 (MQ=255)
     tatGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAAc                                                                                                                                 <  2:403966/141‑1 (MQ=255)
      atGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACt                                                                                                                                >  2:168262/1‑141 (MQ=255)
       tGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCCGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACtt                                                                                                                               >  1:307263/1‑141 (MQ=255)
           cTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTtcatc                                                                                                                           >  2:416848/1‑141 (MQ=255)
              gCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCAtt                                                                                                                        <  1:293962/141‑1 (MQ=255)
               cATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAAGATTAAGCCGTGGTGAAACTTCATCATTg                                                                                                                       >  3:269911/1‑141 (MQ=255)
                aTCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATAGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGTAACTTCATCATTGt                                                                                                                      <  4:28712/141‑1 (MQ=255)
                 tCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGtt                                                                                                                     >  1:541237/1‑141 (MQ=255)
                   tGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCt                                                                                                                   >  2:403348/1‑141 (MQ=255)
                    gTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAATCTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCtt                                                                                                                  <  3:158026/141‑1 (MQ=255)
                    gTTGCTAGAGCTTGTTCGGGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCtt                                                                                                                  <  2:592528/141‑1 (MQ=255)
                     ttGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCttt                                                                                                                 <  2:178117/141‑1 (MQ=255)
                     ttGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCttt                                                                                                                 >  2:50066/1‑141 (MQ=255)
                     ttGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCttt                                                                                                                 >  1:50500/1‑141 (MQ=255)
                      tGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTACTTTa                                                                                                                <  2:529338/141‑1 (MQ=255)
                       gCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAg                                                                                                               >  3:118032/1‑141 (MQ=255)
                           gagCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCt                                                                                                           >  1:544955/1‑141 (MQ=255)
                             gCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCttt                                                                                                         >  1:415035/1‑141 (MQ=255)
                             gCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCttt                                                                                                         >  4:268627/1‑141 (MQ=255)
                              cTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTc                                                                                                        >  2:269601/1‑141 (MQ=255)
                              cTTGATCGTGTTTTTCTGATAAAACATATGTGATGGAAGCCTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCTATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTc                                                                                                        <  2:701278/141‑1 (MQ=255)
                                tGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCa                                                                                                      <  1:20728/141‑1 (MQ=255)
                                tGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCa                                                                                                      <  1:272218/141‑1 (MQ=255)
                                 gTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCaa                                                                                                     <  3:143476/141‑1 (MQ=255)
                                  ttCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAAt                                                                                                    >  2:490046/1‑141 (MQ=255)
                                   tCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAAtt                                                                                                   <  4:100079/141‑1 (MQ=255)
                                    cGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAAttt                                                                                                  <  4:21584/141‑1 (MQ=255)
                                             cTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTaa                                                                                         >  1:409360/1‑141 (MQ=255)
                                              tGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAg                                                                                        >  1:632504/1‑141 (MQ=255)
                                                  aaaaCATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCAc                                                                                    <  3:168021/141‑1 (MQ=255)
                                                    aaCATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGc                                                                                  >  1:349633/1‑141 (MQ=255)
                                                     aCATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGca                                                                                 <  3:49314/141‑1 (MQ=255)
                                                      cATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGcac                                                                                <  3:261302/141‑1 (MQ=255)
                                                        tatGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAg                                                                              <  1:30266/141‑1 (MQ=255)
                                                        tatGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAg                                                                              <  1:204053/141‑1 (MQ=255)
                                                        tatGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACCTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAg                                                                              <  1:318677/141‑1 (MQ=255)
                                                         atGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGt                                                                             <  2:575033/141‑1 (MQ=255)
                                                          tgtgATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGtt                                                                            >  2:79349/1‑141 (MQ=255)
                                                          tgtgATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGtt                                                                            >  3:56908/1‑141 (MQ=255)
                                                           gtgATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTGCTCTATTGCTTATCAATAACATTAAGCCGTGGTGAAAATTCATAATTGTTATTTAGGTCTTTTAAATTTATAGCATAAGCAACGAACAGttt                                                                           >  1:339413/1‑141 (MQ=255)
                                                               tGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTtaat                                                                       <  3:112698/141‑1 (MQ=255)
                                                                ggAAGCGTGAATAATAATGCCTAATGTAACAAAACTTATAATTAATATACTTCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTtaata                                                                      <  2:339413/141‑1 (MQ=255)
                                                                   aGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGt                                                                   >  2:50922/1‑141 (MQ=255)
                                                                     cGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTtg                                                                 <  4:302984/141‑1 (MQ=255)
                                                                      gTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTtgt                                                                <  1:167369/141‑1 (MQ=255)
                                                                         aataataatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGc                                                             >  3:121293/1‑141 (MQ=255)
                                                                           taataatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGctt                                                           <  1:279110/141‑1 (MQ=255)
                                                                             ataatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGcttct                                                         <  1:623112/141‑1 (MQ=255)
                                                                               aatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTa                                                       >  2:147616/1‑141 (MQ=255)
                                                                               aatGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTa                                                       >  1:656974/1‑141 (MQ=255)
                                                                                atGCCTAATGTAAAAAAGCGGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGCAACTTCATCATTGTTCTTTACGTCTTTCCACTTTATAGCCTAAGCCACGAACAGTTTTAATAAGTTGTGGCTTCTTAg                                                      <  4:12686/141‑1 (MQ=255)
                                                                                 tGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAAAATTAAGCCGTGGTGAGACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAgg                                                     >  2:477023/1‑141 (MQ=255)
                                                                                    cTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGAtt                                                  <  2:94829/141‑1 (MQ=255)
                                                                                     tAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTa                                                 <  1:514462/141‑1 (MQ=255)
                                                                                       aTGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTAtc                                               >  1:696653/1‑141 (MQ=255)
                                                                                       aTGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTAtc                                               >  1:696620/1‑141 (MQ=255)
                                                                                         gTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTAtctt                                             <  1:256364/141‑1 (MQ=255)
                                                                                         gTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTAtctt                                             <  2:544955/141‑1 (MQ=255)
                                                                                                   tGATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTTtc                                   >  4:282512/1‑141 (MQ=255)
                                                                                                    gATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCttt                                                                                                         >  2:349324/1‑70 (MQ=255)
                                                                                                    gATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCttt                                                                                                         <  1:349324/70‑1 (MQ=255)
                                                                                                    gATAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTAtct                                  >  3:209391/1‑141 (MQ=255)
                                                                                                     aTAATTAATATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTCGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTGAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTAtctc                                 <  1:608411/141‑1 (MQ=255)
                                                                                                      taattaatATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATctct                                <  1:269601/141‑1 (MQ=255)
                                                                                                        attaatATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTa                              <  1:79349/141‑1 (MQ=255)
                                                                                                          taatATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTaaa                            >  2:458080/1‑141 (MQ=255)
                                                                                                          taatATACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTATGTCTTTCCAATTTATTGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTaaa                            <  3:110296/141‑1 (MQ=255)
                                                                                                            atatACTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATg                          >  4:135309/1‑141 (MQ=255)
                                                                                                                 cTGCTTATCAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGa                     >  2:413848/1‑141 (MQ=255)
                                                                                                                   gCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGata                   <  1:458080/141‑1 (MQ=255)
                                                                                                                         cAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACa             >  1:22493/1‑141 (MQ=255)
                                                                                                                         cAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACa             >  2:485517/1‑141 (MQ=255)
                                                                                                                         cAATAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATATTCTAATTTATCTCTTAAATGACTGATATGTACa             >  1:554990/1‑141 (MQ=255)
                                                                                                                           aTAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATctct                                >  2:402217/1‑120 (MQ=255)
                                                                                                                           aTAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATctct                                <  1:402217/120‑1 (MQ=255)
                                                                                                                           aTAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACATc           >  3:168125/1‑141 (MQ=255)
                                                                                                                           aTAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACATc           >  3:127229/1‑141 (MQ=255)
                                                                                                                           aTAACATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGAATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACATc           >  2:672103/1‑141 (MQ=255)
                                                                                                                              aCATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACATcaac        <  1:403348/141‑1 (MQ=255)
                                                                                                                              aCATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACATcaac        <  1:469289/141‑1 (MQ=255)
                                                                                                                               cATTAAGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACATcaaca       <  1:605553/141‑1 (MQ=255)
                                                                                                                                   aaGCCGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACATCAACAATTc   >  1:222207/1‑141 (MQ=255)
                                                                                                                                    aGCCGTGGTGAAACTTCATCATTGCTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACGGTTTTAATAAGTTGTGGCTTCTTAGGTTTATCTTCTAATTTATCTCTTAAATGCCTGATATGAACATCAACAATTCt  >  1:518739/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                         
AAATATATGATCTTGCATCTGTTGCTAGAGCTTGTTCGTGTTTTTCTGATAAAACATATGTGATGGAAGCGTGAATAATAATGCCTAATGTAACAAAACTGATAATTAATATACTGCTTATCAATAACATTAAGCGGTGGTGAAACTTCATCATTGTTCTTTAGGTCTTTCCAATTTATAGCCTAAGCCACGCACAGTTTTAATAAGTTGTGGCTTCTTAGGATTATCTTCTAATTTATCTCTTAAATGACTGATATGTACATCAACAATTCT  >  USA300TCH1516_ALE/1797637‑1797909

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: