Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,646,545 G→A 6.4% T36T (ACC→ACT ypdF ← Aminopeptidase YpdF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,646,5450GA 6.4% 179.3 / 3.5 68T36T (ACC→ACTypdFAminopeptidase YpdF
Reads supporting (aligned to +/- strand):  ref base G (28/34);  new base A (2/2);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-01

ATTCTTTAATCTCACCAATAATAGTAGATTTACGATTAATAATTTCATAATTTGGCGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCTCATCTTATTACCTCCTAGTGAAACGTTTGTGCATTT  >  USA300TCH1516_ALE/1646409‑1646685
                                                                                                                                        |                                                                                                                                            
aTTCTTTAATCTCACCAATAATAGTAGATTTACGATTAATAATTTCATAATTTGGCGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTAc                                                                                                                                          >  1:914246/1‑141 (MQ=255)
   cTTTAATTTCACCAATAATAGTAGATTTACGATTAATTATTTCATAATTTGGCGTCTGCTTAGTTGCTTGATCAATATATCTAGAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAg                                                                                                                                       <  2:135591/141‑1 (MQ=255)
          ctcACCAATAATAGTCGATCTACGATTAATAATTTCATAATTTGCCGCCTGCTTAGTTGCTTGACCAATATATCTAAAGCCTGTTATCAAATATTGTTTATCTTTAGCTATGCTACGTGCTCCACTCGTACCAGCAAAAcc                                                                                                                                >  2:687835/1‑141 (MQ=255)
          ctcACCAATAATAGTAGATTTACGATTAATAATTTCATAATTTGGCGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAAcc                                                                                                                                >  2:793447/1‑141 (MQ=255)
            cACCAATAATAGTAGATTTACGATTAATAATTTCATAATTTGGCGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTg                                                                                                                              <  1:423809/141‑1 (MQ=255)
            cACCAATAATAGTAGATTTACGATTAATAATTTCATAATTTGGCGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTg                                                                                                                              <  1:440186/141‑1 (MQ=255)
             aCCAATAATAGTAGATTTACGATTAATAATTTCATAATTTGGCGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGa                                                                                                                             <  1:939870/141‑1 (MQ=255)
             aCCAATAATAGTAGATTTACGATTAATAATTTCATAATTTGGCGCCTGCTAAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGGTCCACTGGTACCAGTAAAAccggg                                                                                                                             >  1:31872/1‑138 (MQ=255)
             aCCAATAATAGTAGATTTAAACTTAAAATGGACATGATTTGGCGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGa                                                                                                                             >  1:254096/1‑141 (MQ=255)
                aataatAGTAGATTTACGATTAATAATTTCATAATTTGGCGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATaa                                                                                                                          <  1:887265/141‑1 (MQ=255)
                aataatAGTAGATTTACGATTAATAATTTCATAATTTGGCGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATaa                                                                                                                          <  2:246087/141‑1 (MQ=255)
                      agtagATTTACGATTAATAATTTCATAATTTGGCGCCTGCTTAGTTGATTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATAtct                                                                                                                    >  2:680860/1‑141 (MQ=255)
                        tagATTTACGATTAATAATTTCATAATTTGGCGCCTGCTTAGTTGCTTGATCAATATAGCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATAtcttc                                                                                                                  >  1:640002/1‑141 (MQ=255)
                          gATTTACTATTAATAATTTCTTAATTTGGCGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTa                                                                                                                <  2:18114/141‑1 (MQ=255)
                          gATTTACGATTAATAATTTCATAATTTGGCGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGAGCGCCACTGGTACCAGTAAAACCTGATAAATATCTTCTa                                                                                                                >  1:818929/1‑141 (MQ=255)
                           aTTTACGATTAATAATTTCATAATTTGGCGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTAt                                                                                                               <  2:572244/141‑1 (MQ=255)
                                 gATTAATAATTTCATAATTTGGCGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAAt                                                                                                         <  1:540306/141‑1 (MQ=255)
                                 caTTAATAATTTCATAATTTGACGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAAt                                                                                                         <  2:626495/140‑1 (MQ=255)
                                  tctaataaTTTCATAATTGGGCGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCCTTTGTAATc                                                                                                        <  2:917381/139‑1 (MQ=255)
                                       taatTTCATAATTTGGCGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGaaa                                                                                                   >  2:247085/1‑141 (MQ=255)
                                        aatTTCATAATTTGGCGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAACTATCTTCTATTGTAATCCGAAAg                                                                                                  >  2:463015/1‑141 (MQ=255)
                                             cATAATTTGGCGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGa                                                                                             <  1:970285/141‑1 (MQ=255)
                                              aTAATTTGGCGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATGGTTTATCTTTAGATATGATAAGTGCTCCACTTGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGAt                                                                                            <  1:610395/141‑1 (MQ=255)
                                                         ccTGCTTAGTTGCTTGATCAATTTATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTa                                                                                 <  1:87334/141‑1 (MQ=255)
                                                          cTGCTTAGTTTCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTaa                                                                                <  2:456706/141‑1 (MQ=255)
                                                            gCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAAt                                                                              <  1:693355/141‑1 (MQ=255)
                                                             cTTAGTTGCTTGATCAATTTAACTAGAGTCGTTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATg                                                                             <  2:1011639/141‑1 (MQ=255)
                                                                 gTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGAAATGCCAAGTGCACCCCTGGTACCAGTAAAACCTGATAAATATCTCCTATTGTAATCCGAAAGAAACATAATCACATCTAAATGtttt                                                                         >  2:294977/1‑141 (MQ=255)
                                                                   tGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTg                                                                       >  1:46023/1‑141 (MQ=255)
                                                                     cTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATTTTTTTGtt                                                                     >  2:301870/1‑141 (MQ=255)
                                                                      ttGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGCTTTTCTGGCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTc                                                                    <  1:244799/141‑1 (MQ=255)
                                                                      ttGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACTAGTAAAACCTTATTAATATCTTCTTTTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTc                                                                    <  2:410789/141‑1 (MQ=255)
                                                                      ttGATCAATATATCCCACGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTc                                                                    <  1:929386/141‑1 (MQ=255)
                                                                         aTCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGTTACTAGTAAACCTTGTTAAATATTTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTaa                                                                 <  1:506323/141‑1 (MQ=255)
                                                                          tCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTAGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTaaa                                                                >  2:92387/1‑141 (MQ=255)
                                                                          tCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTAGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTaaa                                                                >  2:506972/1‑141 (MQ=255)
                                                                           aaacatatCCAAAGCCTGTTATCAAATAATGCTCATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTaaaa                                                               <  1:87762/137‑1 (MQ=255)
                                                                                        gTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTtg                                                  <  2:833151/141‑1 (MQ=255)
                                                                                           tGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTgtgt                                               >  2:323929/1‑141 (MQ=255)
                                                                                              tATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATGTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTACAATACGATGCAGTTGTGTTAt                                            >  2:8501/1‑141 (MQ=255)
                                                                                               ttCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTAtt                                           <  1:844277/140‑1 (MQ=255)
                                                                                               aTCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTAtt                                           <  2:914246/141‑1 (MQ=255)
                                                                                                    ctatCGTTTATCTTTTGATATGATAAGTGCTCCACTGGTAACAGTAAAACCTGCTAAATATCTTCTATTGTAATCAGAAAGAATGCTATTCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCt                                      >  2:399995/2‑141 (MQ=255)
                                                                                                    atatTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCt                                      >  1:905086/1‑141 (MQ=255)
                                                                                                    atatTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCt                                      <  1:834610/141‑1 (MQ=255)
                                                                                                    atatTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCt                                      <  2:299238/141‑1 (MQ=255)
                                                                                                        tGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCGTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCTCATc                                  >  1:84075/1‑141 (MQ=255)
                                                                                                          tttATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTTTATAATACGATGCACTTGTGTTATTCTGCTCATCtt                                >  2:707998/1‑141 (MQ=255)
                                                                                                             aTCTTTAGATATGATAAGTGCTCCACTAGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCTCATCttatt                             <  1:92387/141‑1 (MQ=255)
                                                                                                             aTCTTTAGATATGATAAGTGCTCCACTAGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCTCATCttatt                             <  1:506972/141‑1 (MQ=255)
                                                                                                                  tAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTTTGTTATTCTGCTCATCTTATTAcctc                        >  1:389596/1‑141 (MQ=255)
                                                                                                                   aGATATGATAAGTGCTCCACCGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTGTAAAATACGATGCACTTGTGTTATTCTGCTCATCTTATTAcctca                       >  1:745177/1‑140 (MQ=255)
                                                                                                                       ttGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCTCATCTTATTACCTCCTAGt                   <  1:463750/140‑1 (MQ=255)
                                                                                                                       atGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCTCATCTTATTACTTCCTAGt                   >  2:860351/1‑141 (MQ=255)
                                                                                                                       atGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGAAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCTCATCTTATTACCTCCTAGt                   <  1:1040731/141‑1 (MQ=255)
                                                                                                                        tGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCTCATCTTATTACCTCCTAGTg                  <  2:46023/141‑1 (MQ=255)
                                                                                                                        tGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCTCATCTTATTACCTCCTAGTg                  >  1:8518/1‑141 (MQ=255)
                                                                                                                         gATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTACTAAAATACGATGCACTTGTGTTTTTCTGCTCATCTTATTACCTCCTAGCGa                 >  2:859158/1‑141 (MQ=255)
                                                                                                                             aGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCTCATCTTATTACCTCCTAGTGAAACg             >  1:414148/1‑141 (MQ=255)
                                                                                                                                 ctcctcTGGTACCAGTAAAACATGATAAATATCTTCTATTGTAATCCGAAAGAATGATAAGCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCTCATCTTACTACCTCCTAGTGAAACGTTtg         <  1:547332/141‑1 (MQ=255)
                                                                                                                                  tCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTATGCTCATCTTATTACCTCCTAGTGAAACGTTtgt        >  1:435711/1‑141 (MQ=255)
                                                                                                                                   ccACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCTCATCTTATTACCTCCTAGTGAAACGTTtgtg       <  2:56358/141‑1 (MQ=255)
                                                                                                                                   ccACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCTCATCTTATTACCTCCTAGTGAAACGTTtgtg       >  1:128/1‑141 (MQ=255)
                                                                                                                                   ccACTGGTACCAGTAAAACCTCATAAATATCTTCTTTTGTAATCCGAAAGAATGATAATCGCATCTAAATTTTTTTGTTCTAAAATTCGATGCACTTGTATTATTCTGCTTATCTTATTACCTCCTAGTGACACGTTtgtg       <  2:153688/141‑1 (MQ=255)
                                                                                                                                    cACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCTCATCTTATTACCTCCTAGTGAAACGTTTGTGc      <  1:247085/141‑1 (MQ=255)
                                                                                                                                     aCTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCTCATCTTATTACCTCCTAGTGAAACGTTTGTGCa     >  1:1054671/1‑141 (MQ=255)
                                                                                                                                       tGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCTCATCTTATTACCTCCTAGTGCAACGTTTGTGCAtt   <  1:556075/141‑1 (MQ=255)
                                                                                                                                       tGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCTCATCTTATTACCTCCTAGTGAAAGGTTTGTGCAtt   >  1:231900/1‑141 (MQ=255)
                                                                                                                                        ggTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCTCATCTTATTACCTCCTAGTGAAACGTTTGTGCAttt  >  1:701733/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                            
ATTCTTTAATCTCACCAATAATAGTAGATTTACGATTAATAATTTCATAATTTGGCGCCTGCTTAGTTGCTTGATCAATATATCTAAAGTCTGTTATCAAATATTGTTTATCTTTAGATATGATAAGTGCTCCACTGGTACCAGTAAAACCTGATAAATATCTTCTATTGTAATCCGAAAGAATGATAATCGCATCTAAATGTTTTTGTTCTAAAATACGATGCACTTGTGTTATTCTGCTCATCTTATTACCTCCTAGTGAAACGTTTGTGCATTT  >  USA300TCH1516_ALE/1646409‑1646685

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: