Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,830,450 T→C 7.2% T224A (ACT→GCT)  ezrA ← Septation ring formation regulator EzrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,830,4500TC 7.2% 177.8 / 3.4 67T224A (ACT→GCT) ezrASeptation ring formation regulator EzrA
Reads supporting (aligned to +/- strand):  ref base T (30/32);  new base C (3/2);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

TGGTTCAACGAAACTAAGCTCTGTTTTTAAGCTTTGTAATGTACTGTCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGTTGTGCGCTTGTACATAATTACCATCAGCTTTTAGTACTTCATATTGCTCTAACCTTGGCTC  >  USA300TCH1516_ALE/1830310‑1830588
                                                                                                                                            |                                                                                                                                          
tGGTTCAACGAAACTAAGCTCTGTTTTTAAGCTTTGTAATGTACTGTCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGt                                                                                                                                            <  1:74914/141‑1 (MQ=255)
     cAACGAAACTAAGCTCTGTTTTTAAGCTTTGTAATGTACTGTCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGCTtctc                                                                                                                                       >  1:709292/1‑141 (MQ=255)
     cAACGAAACTAAGCTCTGTTTTTAAGCTTTGTAATGTACTGTCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGCTtctc                                                                                                                                       >  1:395807/1‑141 (MQ=255)
           aaCTAAGCTCTGTTTTGAAGCTTTGTAATGTACTGTCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCtaatta                                                                                                                                 <  2:627053/141‑1 (MQ=255)
                aGCTCTGTTTTTAAGCTTTGTAATGTACTGTCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAGTTCTTTTTGTGTTTCTCTaatttattgt                                                                                                                            >  2:888313/1‑139 (MQ=255)
                aGCTCTGTTTTTAAGCTTTGTAATGTACTGTCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTTAGTTTCTCTAATTAATTCt                                                                                                                            >  2:756260/1‑141 (MQ=255)
                  ctctTTTTTTAAGCTTGTTAATGTACTGTCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTgg                                                                                                                          <  2:160697/141‑1 (MQ=255)
                  ctctGTTTTTAAGCTTTGTAATGTACTGTCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTgg                                                                                                                          <  1:826401/141‑1 (MQ=255)
                  ctctGTTTTTAAGCTTTGTAATGTACTGTCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTgg                                                                                                                          <  2:444195/141‑1 (MQ=255)
                  ctctGTTTTTAAGCTTTGTAATGTACTGTCTACTTTCACGTGATCCAGACACGACACGTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTgg                                                                                                                          <  1:26342/141‑1 (MQ=255)
                   tctGTTTTTAAGCTTTGTAATGTACTGTCTACTTTCACGTGATCCAGATCATACCCTTCAAGTTTAAGATCACGGCAACCATATTTAAAATCTTGGAATTGACCAGGTAATTTTTTTTGAGTTTCTCTAATTAATTCTGGt                                                                                                                         <  1:588415/141‑1 (MQ=255)
                    ctGTTTTTAAGCTTTGCAATGTACTGTCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTa                                                                                                                        <  2:786132/141‑1 (MQ=255)
                         tttAAGCTTTGTAATGTACTGTCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTCATTCTGGTATttct                                                                                                                   >  2:473494/1‑141 (MQ=255)
                         tttAAGCTTTGTAATGTACTGTCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTTTttct                                                                                                                   >  2:126396/1‑141 (MQ=255)
                         tttAAGCTTTGTAATGTACTGTCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATttct                                                                                                                   >  2:473489/1‑141 (MQ=255)
                            aaGCTTTGTAATGTACTGTCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTcc                                                                                                                >  2:386563/1‑141 (MQ=255)
                             aGCTTTGTAATGTACTGTCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCa                                                                                                               <  1:259119/141‑1 (MQ=255)
                                        gTACTGTCTACTTTCACGTGATCAAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGGGTTTCTCTAGTTAATTCTGGTATTTCTTCCATATAAGATCtt                                                                                                    >  1:33819/1‑141 (MQ=255)
                                            tGTCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCt                                                                                                >  2:302560/1‑141 (MQ=255)
                                            tGTCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCt                                                                                                <  1:366001/141‑1 (MQ=255)
                                              tCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGt                                                                                              <  1:909902/141‑1 (MQ=255)
                                                tACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCACCCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTGATTATTTCTGGTATTTCTTTCATTTTAGTTCTTAGCTGttt                                                                                            >  2:19635/1‑141 (MQ=255)
                                                tACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTCTTTCTTCCATCTAAGATCTAAGCTCttt                                                                                            >  2:353370/1‑141 (MQ=255)
                                                tACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGttt                                                                                            >  1:516528/1‑141 (MQ=255)
                                                tACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTAAGCTTCTATAATTAATTCTGGTATTTCTTCCGTCTAAGATATTAGCTGttt                                                                                            >  2:450035/1‑141 (MQ=255)
                                                   tttCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAACATCTTAGCTGTTTCAt                                                                                         <  1:465697/141‑1 (MQ=255)
                                                   tttCACGTGACCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCAt                                                                                         <  2:694082/141‑1 (MQ=255)
                                                      cacaTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTTAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTg                                                                                      <  2:841333/141‑1 (MQ=255)
                                                         gTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGttc                                                                                   <  2:348551/141‑1 (MQ=255)
                                                         gTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCCCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGttc                                                                                   >  2:416084/1‑141 (MQ=255)
                                                          tGATCCAGATCATAGCCTTCCACTTTAACATTACGGCAACCACATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTTAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGttca                                                                                  <  2:400129/141‑1 (MQ=255)
                                                            aTCCAGATCATACCTTTCAACTTTAAGAGCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGttcatt                                                                                <  2:472402/141‑1 (MQ=255)
                                                            aTCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAATATCTTAGCTGTTTCATTTGttcatt                                                                                >  2:386777/1‑141 (MQ=255)
                                                            aTCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTTttcatt                                                                                >  1:630332/1‑141 (MQ=255)
                                                            aTCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGttcatt                                                                                >  1:3517/1‑141 (MQ=255)
                                                            aTCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGttaaaa                                                                                >  2:224364/1‑137 (MQ=255)
                                                             tCCATATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGttcattc                                                                               <  2:719708/141‑1 (MQ=255)
                                                             tCCAGATCATACCCTTCAACTTTAAGATCGCGGGACACAGATTTTAAATCTTGGAGTTGACCAGGTAATTCTTTTTGAGTTTCTCTAATGAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTGTCATTTgttcgttc                                                                               >  1:215041/1‑141 (MQ=255)
                                                             tCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGttcattc                                                                               <  1:712259/141‑1 (MQ=255)
                                                             tCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGttcattc                                                                               <  1:712231/141‑1 (MQ=255)
                                                             tCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGttcattc                                                                               <  1:648506/141‑1 (MQ=255)
                                                              ccAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGttcattca                                                                              >  1:215332/1‑141 (MQ=255)
                                                                  aTCATACCCTTCAACATTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGc                                                                          >  2:885885/1‑141 (MQ=255)
                                                                      tACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGAAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTTTTCCATATAATATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCt                                                                      >  1:626165/1‑141 (MQ=255)
                                                                       aCCCTTCATCTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCta                                                                     <  1:236253/141‑1 (MQ=255)
                                                                          cTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCGAGGCAGCTatat                                                                  >  2:820248/1‑141 (MQ=255)
                                                                          cTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTatat                                                                  >  2:725693/1‑141 (MQ=255)
                                                                             cAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGt                                                               <  1:773850/141‑1 (MQ=255)
                                                                                cTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGTtgt                                                            >  2:887242/1‑141 (MQ=255)
                                                                                 tttAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGTtgtg                                                           >  1:317491/1‑141 (MQ=255)
                                                                                  ttgataccACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGTTGTgc                                                          <  2:223758/134‑1 (MQ=255)
                                                                                   tAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGTTGTgcg                                                         <  1:456785/141‑1 (MQ=255)
                                                                                    aaGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGTTGTgcgc                                                        >  1:420300/1‑141 (MQ=255)
                                                                                        tCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTTTTTTTTTTTTATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCCTTCAAGGCCGCTATATGGTTGTGCGCTTGt                                                    >  2:576732/1‑141 (MQ=255)
                                                                                                 ccATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGTTGTGCGCTTGTACATAATTa                                           <  2:849085/141‑1 (MQ=255)
                                                                                                   agatTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGTTGTGCGCTTGTACATAATTAcc                                         <  2:173862/139‑1 (MQ=255)
                                                                                                      ttttAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGCTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGTTGTGCGCTTGTACATAATTACCATc                                      <  2:709292/141‑1 (MQ=255)
                                                                                                      ttttAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGCTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGTTGTGCGCTTGTACATAATTACCATc                                      <  2:395807/141‑1 (MQ=255)
                                                                                                             tCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGTTGTGCGCTTGTACATAATTACCATCAGCtttt                               <  2:58177/141‑1 (MQ=255)
                                                                                                                 ggAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGTTGTGCGCTTGTACATGATTACCATCCACTGTTAGtc                           >  2:716316/1‑140 (MQ=255)
                                                                                                                  gAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGTTGTGCGCTTGTACATAATTACCATCAGCTTTTAGTAc                          <  1:841551/141‑1 (MQ=255)
                                                                                                                  gAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGTTGTGCGCTTGTACATAATTACCATCAGCTTTTAGTAc                          <  1:237559/141‑1 (MQ=255)
                                                                                                                     ttGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGTTGTGCGCTTGTACATAATTACCATCAGCTTTTAGTACTTc                       <  2:808668/141‑1 (MQ=255)
                                                                                                                       gACCAGGTAATTCTTTTTGAGTTTCTCTAATTAGTTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCACCTATATGGTTGTGCGCTTGTACATAATTACCATCAGCTTTTAGTACTTCat                     >  2:602500/1‑141 (MQ=255)
                                                                                                                          cAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATCAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGTTGTGCGCTTGTACATAATTACCATCAGCTTTTAGTACTTCATAtt                  <  1:171425/141‑1 (MQ=255)
                                                                                                                          cAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGTTGTGCGCTTGTACATAATTACCATCAGCTTTTAGTACTTCATAtt                  >  1:317889/1‑141 (MQ=255)
                                                                                                                                aTTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGTTGTGCGCTTGTACATAATTACCATCAGCTTTTAGTACTTCATATTGCTCTa            >  2:563215/1‑141 (MQ=255)
                                                                                                                                          aGTTTCTCTAATTAATTATCGGACTTCGTCCATATAATAGCGTAGCGGTGTCATTGCTTCATTCAAGGCAGCTATATGGTTGTGCGCTTGTACATAATTACCATCAGCTTTTAGTACTTCATATTGCTCTAACCTTGGCTc  <  2:626165/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                          
TGGTTCAACGAAACTAAGCTCTGTTTTTAAGCTTTGTAATGTACTGTCTACTTTCACGTGATCCAGATCATACCCTTCAACTTTAAGATCACGGCAACCATATTTTAAATCTTGGAATTGACCAGGTAATTCTTTTTGAGTTTCTCTAATTAATTCTGGTATTTCTTCCATATAAGATCTTAGCTGTTTCATTTGTTCATTCAAGGCAGCTATATGGTTGTGCGCTTGTACATAATTACCATCAGCTTTTAGTACTTCATATTGCTCTAACCTTGGCTC  >  USA300TCH1516_ALE/1830310‑1830588

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