Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,516,802 A→T 9.0% intergenic (‑8/+6) gpsB ← / ← USA300HOU_RS07400 Cell cycle protein GpsB/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,516,8020AT 9.0% 178.0 / 2.2 71intergenic (‑8/+6)gpsB/USA300HOU_RS07400Cell cycle protein GpsB/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/37);  new base T (5/2);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTTTGATAATCAGCAATAATGTCATCTAAAAATGCATCTACTTCTTCTCTTCTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACAATATCACAAGTATAGTTTGTTTTATCCATAAAATCAACTAACATCTGCTTGAAGAACTTAGGACTCGC  >  USA300TCH1516_ALE/1516662‑1516942
                                                                                                                                            |                                                                                                                                            
ttttttGATAATCAGCAATAATGTCATCTAAAAATGCATCTACTTCTTCTCTTCTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATCAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCa                                                                                                                                              >  1:480005/1‑141 (MQ=255)
ttttttGATAATCAGCAATAATGTCATCTAAAAATGCATCTACTTCTTCTCTTCTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCa                                                                                                                                              >  1:790212/1‑141 (MQ=255)
     tGATAATCAGCAATAATGTCATCTAAAAATGCATCTACTTCTTCTCTTCTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCa                                                                                                                                         >  1:239762/1‑141 (MQ=255)
       aTAATCAGCAATAATGTCATCTAAAAATGCATCTACTTCTTCTCTTCTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCAtt                                                                                                                                       <  2:375956/141‑1 (MQ=255)
          aTCAGCAATAATGTCATCTAAAAATGCATCTACTTCTTCTCTTCTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATGTCAATGAAACATCTGACATTTTTTCCACCTCATTAGa                                                                                                                                    <  2:395011/141‑1 (MQ=255)
             aGCAATAATGTCATCTAAAAATGCATCTACTTCTTCTCTTCTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAAc                                                                                                                                 <  1:711035/141‑1 (MQ=255)
             aGCAATAATGTCATCTAAAAATGCATCTACTTCTTCTCTTCTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAAc                                                                                                                                 >  1:603149/1‑141 (MQ=255)
                      tttatttAAAAATGCATCTACTTCTTCTCTTCTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCtt                                                                                                                        <  1:269369/135‑1 (MQ=255)
                      gTCATCTAAAAATGCATCTACTTCTTCTCTTCTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCtt                                                                                                                        <  1:143775/141‑1 (MQ=255)
                         aTCTAAAAATGCATCTACTTCTTCTCTTCTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAg                                                                                                                     >  1:429455/1‑141 (MQ=255)
                             aaaaaTGCATCTACTTCTTCTCTTCTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCa                                                                                                                 <  2:180108/141‑1 (MQ=255)
                              acaaTGCAGCTACTTCTTCTCTTCTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATCTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCAc                                                                                                                <  2:210138/139‑1 (MQ=255)
                               aaaTGCATCTACTTCTTCTCTTCTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACt                                                                                                               >  1:128300/1‑141 (MQ=255)
                                 aTGCATCTACTTCTTCTCTTCTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCCGACCACTGt                                                                                                             >  2:64950/1‑141 (MQ=255)
                                   gCATCTACTTCTTCTCTTCTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTaa                                                                                                           <  2:753355/141‑1 (MQ=255)
                                     aTCTACTTCTTCTCTTCTATAGCCACGAGCCATCGTTTTTTTTTTTTTTTTTTCATAAAAAGCTGTTGCTGATAATTTTAATGAAACATCTGACATTTTTTCGTGTTCTTTAGAAACTTTGAAATTCAAACCACTGttatt                                                                                                         >  2:165216/1‑137 (MQ=255)
                                       ttACTGCTTCTCTTCTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCa                                                                                                       <  2:142385/140‑1 (MQ=255)
                                       cTACTTCTTCTCTTCTATAGTCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCa                                                                                                       <  2:428388/141‑1 (MQ=255)
                                           ttcttcTCTTCTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATtga                                                                                                   <  1:442790/141‑1 (MQ=255)
                                            tcttctcttcTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATtgat                                                                                                  <  2:654044/141‑1 (MQ=255)
                                                 tcttctATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGCAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTgatggatg                                                                                             >  1:68756/1‑141 (MQ=255)
                                                 tcttctATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCCCCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATg                                                                                             >  1:91487/1‑141 (MQ=255)
                                                    tctATAGCCACGAGCCATCGTTTTTTGAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTCATGAATGCTg                                                                                          <  1:578045/141‑1 (MQ=255)
                                                         aGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAAt                                                                                     <  1:40768/141‑1 (MQ=255)
                                                         aGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAAt                                                                                     >  1:164743/1‑141 (MQ=255)
                                                         aGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTAATTTTATGCTTGGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAAt                                                                                     <  1:742392/141‑1 (MQ=255)
                                                          gCCACGAGCCATCGTTTTTTCAATATCTTTTTCATAAATATCTTCTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATCAGAAACTTTGATCTTCACACCACTGTAACTCATTTATGAATGCTGTTAAtt                                                                                    >  2:122049/1‑141 (MQ=255)
                                                              cGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATc                                                                                <  1:385971/141‑1 (MQ=255)
                                                              cGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTCTGATCTTCAGACCACTGTAAGTCATTGATGAATTCTGTTAATTCATc                                                                                >  1:594621/1‑141 (MQ=255)
                                                               gAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTCTAAGTCATTGATGAATGCTGTTAATTCCTCg                                                                               >  1:16457/1‑141 (MQ=255)
                                                                ggCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGGTAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAAGTTGGGTGGTGGGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGa                                                                              <  2:2632/140‑1 (MQ=255)
                                                                aGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGTTGAATGCTGTTAATTCATCGa                                                                              >  1:460808/1‑141 (MQ=255)
                                                                aGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGGATGCTGTTAATTCATCGa                                                                              >  1:35806/1‑141 (MQ=255)
                                                                 gCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGaa                                                                             >  2:253304/1‑141 (MQ=255)
                                                                 gCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGaa                                                                             >  1:753588/1‑141 (MQ=255)
                                                                  ccATCGTTTTTTTTTTTTTTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTCGAAACTTTGATCTTCAGACCACTTTTATTCTTTTTTGTTTGCTGTTAATTCATCGAAc                                                                            >  2:430674/1‑141 (MQ=255)
                                                                   gagcGTGTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACg                                                                           <  1:176132/138‑1 (MQ=255)
                                                                   cATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACg                                                                           <  1:109217/141‑1 (MQ=255)
                                                                    aTCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGt                                                                          >  1:308208/1‑141 (MQ=255)
                                                                    aTCGTTTTTTCAAAATCTTTTTCATAAATAACTTTTGCTGATAATTTCAATGTAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGt                                                                          <  2:495035/141‑1 (MQ=255)
                                                                      cGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTca                                                                        <  2:207171/141‑1 (MQ=255)
                                                                        ttttttCAAATTCTGTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTcaca                                                                      <  2:632068/141‑1 (MQ=255)
                                                                        ttttttCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTcaca                                                                      <  2:73777/141‑1 (MQ=255)
                                                                         tttttCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACaa                                                                     <  1:9027/141‑1 (MQ=255)
                                                                           tttCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACAata                                                                   <  1:475594/141‑1 (MQ=255)
                                                                              cAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCGGACCACTGTAAGTCATTGCTGAATGCTGTTAGTTCATCGAACGTCAACATATcc                                                                >  1:67607/1‑140 (MQ=255)
                                                                                    cTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACAATATCACAAGta                                                          <  2:430204/141‑1 (MQ=255)
                                                                                      ttttCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTTTTTTTTTTTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTGATCGAACCTCACATTATCACAAGtttt                                                        >  2:600017/1‑138 (MQ=255)
                                                                                       tttCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCTAACGTCACAATATCACAAGTATAg                                                       >  1:749414/1‑141 (MQ=255)
                                                                                        ttCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTTTTTTTTTTTTTAAACTTTGATCTTCAAACCACTTTAAGTCATTGATGAATGCTGTTAATTCATCGATCGTCCCCCTATCACAAGtttact                                                      >  2:745551/1‑139 (MQ=255)
                                                                                        ttCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACTTCACAATACCACAAGTATAgt                                                      >  1:656787/1‑141 (MQ=255)
                                                                                        ttCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACAATATCACAAGTATAgt                                                      <  2:621978/141‑1 (MQ=255)
                                                                                           aTAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACAATATCACAAGTATAgtttt                                                   >  1:443315/1‑140 (MQ=255)
                                                                                             aaaTATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACAATATCACAAGTATAgtttgtt                                                 <  1:386978/141‑1 (MQ=255)
                                                                                                  tCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCGTATAGGACACGTTGGTCGTCATAGCACTATAAGAAATTGATGAATGTTGTTTATTCATTGAACGTCACAATATCACAGGGATAGTTTGTTTTATc                                            >  2:44006/1‑141 (MQ=255)
                                                                                                   tttttGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACAATATCACAAGTATAGTTTGTTTTATcc                                           <  2:271227/140‑1 (MQ=255)
                                                                                                   cTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTTTTTTTTTTTTGAAAATTTGATCTTCAGACCACTGTAAGGCATTGATGAGTGTTGTTAATTGTTCGCACCTCACAATATCAACAGTATAGTTTGTTTTaact                                           >  2:177189/1‑138 (MQ=255)
                                                                                                   cTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAACCTTTGATTTTCAGACCACTGTAAGTCATTGATGAATGCAGTTAATTCATCGAACGTCACAATATAACAAGTCAAGTTTGGTTTCTcc                                           >  2:635639/1‑141 (MQ=255)
                                                                                                         cTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACAATATCACAAGTATAGTTTGTTTTATCCATaaaa                                     >  1:76325/1‑141 (MQ=255)
                                                                                                          tGATAATTTCAATGAAACATCTGACATTTTTTTTTTTTTTTTTTAAACTTTGATCTTCAGACCACTTTAATTGATTGATGGATGTTGTTACTTTTTCGCTTGTCACATTTTCTCTAGTATATTTTGTTTTATCCTTAAAAt                                    >  2:19631/1‑141 (MQ=255)
                                                                                                               aTTTCAATGAAACATCTGACATTTTTTCCTCCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACAATATCACAAGTATAGTTTGTTTTATCCATAAAATCAACt                               <  2:472065/141‑1 (MQ=255)
                                                                                                               aTTTCAATGAAACATCTGACATTTTTTCCAGCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACAATATCACAAGTATAGTTTGTTTTATCCATAAAATCAACt                               <  2:500387/141‑1 (MQ=255)
                                                                                                                 ttCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACAATATCACAAGTATAGTTTGTTTTATCCATAAAATCAACTaa                             >  2:66238/1‑141 (MQ=255)
                                                                                                                     aTGAAACATCTGACATTTTTTCCAGCTGATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACAATATCAAAAGTATAGTTTGTTTTATCAATACAATCAACTAACATc                         >  2:559958/1‑141 (MQ=255)
                                                                                                                       gAAACATCTTACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACAATATCACAAGTATAGTTTGTTTTATCCATAAAATCAACTAACATCTg                       <  2:270132/141‑1 (MQ=255)
                                                                                                                        aaaCATCTGACATTTTTTCCTCCTCATTAGATACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGACCGTCTCAATATCACAAGTATAGTTTGTTTTATCCATAAAATCAACTAACATCTGc                      <  2:86951/141‑1 (MQ=255)
                                                                                                                         aaCATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACAATATCACAAGTATAGTTTGTTTTATCCATAAAATCAACTAACATCTGCt                     >  1:463772/1‑141 (MQ=255)
                                                                                                                         aaCATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACAATATCACAAGTATAGTTTGTTTTATCCATAAAATCAACTAACATCTGCt                     <  1:138666/141‑1 (MQ=255)
                                                                                                                         aaCATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACAGTAAGTCATTGACGAATGCTGTTAATTCATCGAACGTCACAATATCACAAGTATAGTTTGTTTTATCCATAAAATTAACTAACATCTGCt                     <  2:539080/141‑1 (MQ=255)
                                                                                                                           cATCTGACCTTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACAATATCACAAGTATAGTTTGTTTTATCCATAAAATCATCTAACATCTGCTTg                   <  1:436539/141‑1 (MQ=255)
                                                                                                                                   cttctttCCACCTCATTAGCAACTTTGATCTTCAGACCACTGTAAGTCATTGTTGAATGCTGTTAATTCATCGAACGTCACAATATCACAAGTATAGTTTGTTTTATCCATAAAATCAACCAACATCTGCTTGAAGAACtt           <  2:451938/137‑1 (MQ=255)
                                                                                                                                      ttttCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACAATATCACAAGTATAGTTTGTTTTATCCATAAAATCAACTAACATCTGCTTGAAGAACTTAgg        <  1:125662/141‑1 (MQ=255)
                                                                                                                                       tctCCACCTCATTAGAAACTTTTATCTTCAGCCAACTGTAGATCTTTAATTTCTCCTGCTAATTCTCCGCATTTCTCATTATTACAATTATAGTTTGTTTTATCCGTAAAATCAACTAACATTTGTTTTAAAAACTTAGGa       <  2:86192/139‑1 (MQ=255)
                                                                                                                                            aCCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACAATATCACAAGTATAGTTTGTTTTATCCATAAAATCAACTAACATCTGCTTGAAGAACTTAGGACTCGc  <  1:430674/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                            
TTTTTTGATAATCAGCAATAATGTCATCTAAAAATGCATCTACTTCTTCTCTTCTATAGCCACGAGCCATCGTTTTTTCAAAATCTTTTTCATAAATATCTTTTGCTGATAATTTCAATGAAACATCTGACATTTTTTCCACCTCATTAGAAACTTTGATCTTCAGACCACTGTAAGTCATTGATGAATGCTGTTAATTCATCGAACGTCACAATATCACAAGTATAGTTTGTTTTATCCATAAAATCAACTAACATCTGCTTGAAGAACTTAGGACTCGC  >  USA300TCH1516_ALE/1516662‑1516942

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