Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,213,105 G→A 5.6% F105F (TTC→TTT thiE ← Thiamine‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,213,1050GA 5.6% 204.2 / 3.5 71F105F (TTC→TTTthiEThiamine‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (39/28);  new base A (2/2);  total (41/30)
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.60e-01

CTACTGGTATATGCGCGTCATGCTTCGACGGCGTAGGATAGATTGGTCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATTGATGACAGAGATGCTGTAATTCT  >  USA300TCH1516_ALE/2212972‑2213243
                                                                                                                                     |                                                                                                                                          
cTACTGGTATATGCGCGTCATGCTTGGACGGCGTAGGATAGATTGGTCCCACACCTATACAAGCTACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTACTACTCAGTCCAATAATTTTCTCCCTGACATAtt                                                                                                                                     >  2:138684/1‑141 (MQ=255)
 tACTGGTATATGCGCGTCATGCTTCGACGGCGTAGGATAGATTGGTCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATtg                                                                                                                                    <  1:138073/141‑1 (MQ=255)
 gattGGTATTTGCGAGTCATGATTCTACGGCTTAGTATAGATTGGTCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAATTCCAATAATTTTATCAGTGAAATATtg                                                                                                                                    <  2:560456/138‑1 (MQ=255)
          atGCGCGTCATGCTTCGACGGCGTAGGATAGATTGGTCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTAtctc                                                                                                                           >  2:98794/1‑141 (MQ=255)
           tGCGCGTCATGCTTCGACGGCGTAGGATAGATTGGTCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGCAATATTGTGCTATctct                                                                                                                          >  1:28587/1‑141 (MQ=255)
           tGCGCGTCATGCTTCGACGGCGTAGGATAGATTGGTCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATctct                                                                                                                          >  1:189880/1‑141 (MQ=255)
             cgcgTCATGCTTCGACGGCGTAGGATAGATTGGTCCAACACCTATATAATCAACATGAGTTAACTCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTCTCAGTGAAATATTGTGCTATCTCttt                                                                                                                        >  1:100766/1‑141 (MQ=255)
                     tcTTCGACGGCGTAGGATAGATTGGTCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATTAGTGAAATATTGTGCTATCTCTTTGACTTTCg                                                                                                                <  2:91508/140‑1 (MQ=255)
                     gcgtCGACGGCGTAGGATAGATTGGTCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCg                                                                                                                <  1:310603/138‑1 (MQ=255)
                     gCTTCGACGGCGTAGGATAGATTGGTCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTATGCTATCTCTTTGACTTTct                                                                                                                >  2:631539/1‑140 (MQ=255)
                      cTTCGACGGCGTAGGATAGATTGGTCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGc                                                                                                               >  1:209540/1‑141 (MQ=255)
                            cgtcGTAGGATAGATTGGTCCAACACCTATATCATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCatcatc                                                                                                         <  2:469977/138‑1 (MQ=255)
                            cGGCGTAGGATAGATTGGTCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCGAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCatcatc                                                                                                         >  1:174183/1‑141 (MQ=255)
                            cGGCGTAGGATAGATTGGTCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCatcatc                                                                                                         >  1:100846/1‑141 (MQ=255)
                                gTAGGATAGATTGGTCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGa                                                                                                     >  2:463792/1‑141 (MQ=255)
                                gTAGGATAGATTGGTCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGa                                                                                                     >  1:559520/1‑141 (MQ=255)
                                gTAGGATAGATTGGTCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGa                                                                                                     <  1:342073/141‑1 (MQ=255)
                                   ggATAGATTGGTCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCg                                                                                                  <  1:137281/141‑1 (MQ=255)
                                    gATAGATTGGTCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGa                                                                                                 >  1:294622/1‑141 (MQ=255)
                                              tCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATAcc                                                                                       <  1:564900/141‑1 (MQ=255)
                                              tCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATAcc                                                                                       <  1:557384/141‑1 (MQ=255)
                                              accAACACCTATATAATCGACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATAcc                                                                                       <  2:99704/140‑1 (MQ=255)
                                                      cTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGCCATGTATACCACCCGCGt                                                                               >  2:392873/1‑141 (MQ=255)
                                                       tatataATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGtt                                                                              <  1:385664/141‑1 (MQ=255)
                                                             aTCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCAGCGTTAATTTc                                                                        >  1:13041/1‑141 (MQ=255)
                                                               cAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTTACCGACATGTATACCATCCGCGTTAATTTCtt                                                                      >  1:15794/1‑141 (MQ=255)
                                                               cAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTAAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCtt                                                                      >  2:401752/1‑141 (MQ=255)
                                                               cAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTAAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCtt                                                                      >  2:218445/1‑141 (MQ=255)
                                                                 aCATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCtttt                                                                    <  1:98794/141‑1 (MQ=255)
                                                                  cATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTg                                                                   >  1:154945/1‑141 (MQ=255)
                                                                  cATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTg                                                                   >  1:454537/1‑141 (MQ=255)
                                                                   aTGAGTTAAATCAGATTTCGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGc                                                                  <  1:248594/141‑1 (MQ=255)
                                                                     gggatAAACCAGATATTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTATATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCa                                                                <  1:462744/137‑1 (MQ=255)
                                                                      aGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCaa                                                               >  1:63894/1‑141 (MQ=255)
                                                                      aGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCaa                                                               >  1:325627/1‑141 (MQ=255)
                                                                       gTTAAATCAGATTTTGCATACTCATCCAAATGACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTTACCGACATGTATACCATCAGCGTTAATTTCTTTTGCCAAg                                                              >  2:232547/1‑141 (MQ=255)
                                                                           aaTCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATa                                                          >  1:447902/1‑141 (MQ=255)
                                                                              cAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATAcat                                                       <  1:380603/141‑1 (MQ=255)
                                                                              cAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATAcat                                                       >  1:509315/1‑141 (MQ=255)
                                                                               aGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATAcatc                                                      <  1:455819/141‑1 (MQ=255)
                                                                                gATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGAACGACAGATATACCATCCACGTTCATTTCTTTTGCCAAGGATAcctca                                                     >  2:453532/1‑141 (MQ=255)
                                                                                     tGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCAtt                                                >  1:512279/1‑141 (MQ=255)
                                                                                           cTCACCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCTTCATCTTGTCTGATATGTATACCATCTGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGAt                                          <  1:594097/141‑1 (MQ=255)
                                                                                             cATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGCATACATCATCATTAACGATaa                                        >  2:221183/1‑141 (MQ=255)
                                                                                               tCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATataa                                      >  1:518736/1‑141 (MQ=255)
                                                                                                     tCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATTAAAGGcaca                                >  2:447366/1‑141 (MQ=255)
                                                                                                     tCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGcaca                                <  2:497966/141‑1 (MQ=255)
                                                                                                             aCTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCACGGATACATCATCATTAACGATACAACGCACATCATATtg                        >  2:255000/1‑141 (MQ=255)
                                                                                                              cTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATCAAAGGCACATCATATtga                       >  1:210665/1‑141 (MQ=255)
                                                                                                              cTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATtga                       >  1:511848/1‑141 (MQ=255)
                                                                                                              cTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATtga                       <  2:277912/141‑1 (MQ=255)
                                                                                                               tAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTCACGATAAAAGGCACATCCTATtgat                      >  1:576867/1‑141 (MQ=255)
                                                                                                               tAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATtgat                      <  1:463792/141‑1 (MQ=255)
                                                                                                                 aGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATtgatga                    >  1:364966/1‑141 (MQ=255)
                                                                                                                 aGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCGTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATtgatga                    >  1:324571/1‑141 (MQ=255)
                                                                                                                      aataatTTTATCAGTAAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATTGATGACagag               <  1:401752/141‑1 (MQ=255)
                                                                                                                      aataatTTTATCAGTAAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATTGATGACagag               <  1:218445/141‑1 (MQ=255)
                                                                                                                       ataatTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATTGATGACAgaga              >  2:419912/1‑141 (MQ=255)
                                                                                                                       ataatTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATTGATGACAgaga              >  1:421765/1‑141 (MQ=255)
                                                                                                                         aatTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAAGGATAAAAGGCACATGATATTGATGACAgagagg            >  1:197751/1‑139 (MQ=255)
                                                                                                                         aatTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATTTATGACagcaatg            >  1:614145/1‑136 (MQ=255)
                                                                                                                            tttGTCACTGAAATAGTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATTGATGACAGAGATGCTg         <  2:189880/141‑1 (MQ=255)
                                                                                                                            tttATCAGTGAAAGATTGTGCTATCGCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATTGATGACAGAGATGCTg         <  2:28587/141‑1 (MQ=255)
                                                                                                                             ttctgcgTGAAATGTTGTGCTATCTCTGTGATCGTCGGATCATCTTGACCGACATGTATACCATCCGCGTTAAGTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATTGATGACAGAGATGCTGt        <  2:100766/135‑1 (MQ=255)
                                                                                                                             ttATCAGTGAAATAGTGCGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATTGATGACAGAGATGCTGt        <  1:48643/141‑1 (MQ=255)
                                                                                                                              tATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATTGATGACAGAGATGCTGTa       >  2:347370/1‑141 (MQ=255)
                                                                                                                               aTCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGAGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATTGATGACAGAGATGCTGTaa      <  2:246725/141‑1 (MQ=255)
                                                                                                                                tCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATTGATGACAGAGATGCTGTAAt     <  1:414678/141‑1 (MQ=255)
                                                                                                                                tCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATTGATGACAGAGATGCTGTAAt     <  2:509315/141‑1 (MQ=255)
                                                                                                                                  aGTGAGATATTGTGCTATCGCTTTGAGTTTCGCATCTTCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATTGATGACAGAGATGCTGTAATTc   <  2:282127/141‑1 (MQ=255)
                                                                                                                                   gTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATTGATGACAGAGATGCTGTAATTCt  >  1:28305/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                                          
CTACTGGTATATGCGCGTCATGCTTCGACGGCGTAGGATAGATTGGTCCAACACCTATATAATCAACATGAGTTAAATCAGATTTTGCATACTCATCTAAATCACTAATACTAAGTCCAATAATTTTATCAGTGAAATATTGTGCTATCTCTTTGACTTTCGCATCATCTTGACCGACATGTATACCATCCGCGTTAATTTCTTTTGCCAAGGATACATCATCATTAACGATAAAAGGCACATCATATTGATGACAGAGATGCTGTAATTCT  >  USA300TCH1516_ALE/2212972‑2213243

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