Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 638,614 G→T T14K (ACA→AAA)  pdxK ← Pyridoxine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE638,6140GT99.5% 608.1 / NA 184T14K (ACA→AAA) pdxKPyridoxine kinase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (66/117);  total (67/117)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 184 total aligned reads displayed.

ATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATGAATCATT  >  USA300TCH1516_ALE/638478‑638758
                                                                                                                                        |                                                                                                                                                
aTCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCgg                                                                                                                                      <  2:815856/149‑1 (MQ=255)
  ccATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCa                                                                                                                                    <  2:338377/149‑1 (MQ=255)
     ttGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAAtt                                                                                                                                 >  2:229134/1‑149 (MQ=255)
        gaaaTGGTGTCACATCGTGTGACCATGTCTCTTTATCCATAGTAACGATGGCAGTCAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCCCTTTGCATACCTGCGCCAGCACTTTTGTCAGACCCGGCAATTgct                                                                                                                              <  2:865253/147‑3 (MQ=255)
        gTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATCCCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGtt                                                                                                                              >  1:846501/1‑149 (MQ=255)
        gTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGtt                                                                                                                              <  2:1101549/149‑1 (MQ=255)
           aTGGTGTAACATCGTGTGACTATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTaaa                                                                                                                           <  2:1409551/149‑1 (MQ=255)
           aTGGTGTAACATCGTGTGACTAGGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAAAGTTATAAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTaaa                                                                                                                           <  1:334545/149‑1 (MQ=255)
           aTGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTaaa                                                                                                                           <  1:1268958/149‑1 (MQ=255)
             ggTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAAc                                                                                                                         <  1:1299782/149‑1 (MQ=255)
               tgtAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActt                                                                                                                       <  2:1312084/149‑1 (MQ=255)
                gtAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActtt                                                                                                                      <  1:132864/149‑1 (MQ=255)
                   aCATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActttctt                                                                                                                   <  2:1278651/149‑1 (MQ=255)
                    cATCGTGTGACCATTTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActttcttt                                                                                                                  >  2:954562/1‑149 (MQ=255)
                    cATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAActttcttt                                                                                                                  >  2:898880/1‑149 (MQ=255)
                     aTCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCGTACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTa                                                                                                                 <  1:822464/149‑1 (MQ=255)
                      gcGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTaa                                                                                                                <  1:704151/148‑1 (MQ=255)
                         tgtgACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTaaacc                                                                                                             >  1:1411038/1‑147 (MQ=255)
                         tgtgACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCAGACGTATCTAATTCTTGGAACGTTTTCAAATCGGCTTGCATACCTGCGCCAGCACTTTGGTCAGAGCCGGCAATTGTTAAAACTTTCTTTAAAGc                                                                                                             <  2:625896/149‑1 (MQ=255)
                           tgACCATGTATCTTTATCCCTAGTAACGATGGCAGTTAAAGCGACCATGCGATACGTATCTAATTCTTGGAACGTTTTCAAATCTGATTGCATACCTGCGCCAGCACTTTTGTCAGCACCGGCAATTGTTAAAACTTTCTTTAAAGCCa                                                                                                           >  2:467679/1‑149 (MQ=255)
                           tgACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTATTGGAACGTTTTCAAATCTGCTTGGATACCTGCGCCATCACGTGTATCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCa                                                                                                           >  2:526268/1‑149 (MQ=255)
                                      cTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCAc                                                                                                <  1:1587200/149‑1 (MQ=255)
                                      cTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCAc                                                                                                <  1:619342/149‑1 (MQ=255)
                                      cTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCAc                                                                                                <  2:152501/149‑1 (MQ=255)
                                      cTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTGGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCAc                                                                                                <  1:1177383/149‑1 (MQ=255)
                                      cTTTATCCAGAGTCACGATTCTAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGCACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCAc                                                                                                <  1:646974/149‑1 (MQ=255)
                                       tttATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCACTTGTTAAACCTTTCTTTAAAGCCATTGAGCTCCAcc                                                                                               >  2:1326896/1‑148 (MQ=255)
                                       tttATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACt                                                                                               >  2:38265/1‑149 (MQ=255)
                                       tttATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGAGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACt                                                                                               <  1:1145085/149‑1 (MQ=255)
                                        ttATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAGCTTTCTTTAAAGCCATTGAGCTTCACTc                                                                                              <  1:1046186/149‑1 (MQ=255)
                                        ttATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTc                                                                                              <  1:561263/149‑1 (MQ=255)
                                        taacccATAGGTACGATCGCAGGTAAAGGGACCATGCCTTACGGATCTAATTCTTCGATCGTTTTCAAATCTGCTTGCATACCTGCGCCCGCACTTTTGTCAGAACCGGCCATTGTTAAAACTTTCTTTACAGCCATTGAGCTTCACTc                                                                                              <  2:1496319/145‑1 (MQ=255)
                                          aTCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAATGCGGTTGAGCTTCACTCCt                                                                                            <  2:114900/149‑1 (MQ=255)
                                          aTCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCt                                                                                            <  2:683537/149‑1 (MQ=255)
                                          aTCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCt                                                                                            <  1:744956/149‑1 (MQ=255)
                                          aTCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCt                                                                                            <  1:729738/149‑1 (MQ=255)
                                           tCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTa                                                                                           <  1:1481075/149‑1 (MQ=255)
                                            ccATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCTGCAATTGTTAAAACTTTCTTTAAAACCATTGAGCTTATCTCTTAc                                                                                          >  2:12312/1‑149 (MQ=255)
                                            ccATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTAc                                                                                          <  1:1143709/149‑1 (MQ=255)
                                            ccACAGTAACGACTGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTAc                                                                                          <  2:704324/149‑1 (MQ=255)
                                             cATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCTCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACa                                                                                         <  2:867470/149‑1 (MQ=255)
                                             cATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCTCTCCTACa                                                                                         <  2:872008/149‑1 (MQ=255)
                                                 gTAACGATGGCAGTTAAAACGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACataat                                                                                     >  2:1163061/1‑149 (MQ=255)
                                                  tAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAGTTGTTAAAACTTTCTTTTACGCCATTGCGCTTCACTCCTGCataata                                                                                    >  1:347417/1‑149 (MQ=255)
                                                    aCGATGGCAGGTAAAGCGACCATGCCATACGTATCTAGTTCTTGGGACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACatactaat                                                                                  >  2:1703650/1‑149 (MQ=255)
                                                     cGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCATACataataata                                                                                 >  2:789461/1‑149 (MQ=255)
                                                      gATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAatat                                                                                <  1:1016934/149‑1 (MQ=255)
                                                      gATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAatat                                                                                <  1:102812/149‑1 (MQ=255)
                                                      gATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAatat                                                                                <  1:378992/149‑1 (MQ=255)
                                                       aTGGCAGTTAAAGCGACAATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATAtt                                                                               <  2:643201/149‑1 (MQ=255)
                                                          gCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAACACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTa                                                                            >  1:524099/1‑149 (MQ=255)
                                                                  aGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAACGCCATTGAGCTTCACTCCTACATAATAATATTGTATTcatcat                                                                    >  2:1339712/1‑149 (MQ=255)
                                                                  aGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCGATTGAGCTTCACTCCTACATAATAATATTGTATTcatcat                                                                    >  1:1218583/1‑149 (MQ=255)
                                                                  aGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATGTTGTATTcatcat                                                                    >  2:1003378/1‑149 (MQ=255)
                                                                  aGCGACCATGCCATACGTATCTAATTCATGGAACGTTTTCAAATCTGCTTGCTTACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTcatcat                                                                    >  1:944639/1‑149 (MQ=255)
                                                                  aGCGACCATGCCATACGCATCTCAGTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCCCTCCTCCATAATAATATTGTATTcatcat                                                                    >  2:114517/1‑149 (MQ=255)
                                                                   gCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcata                                                                   <  2:1578283/149‑1 (MQ=255)
                                                                   gCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcata                                                                   >  2:771884/1‑149 (MQ=255)
                                                                   gCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTATCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATATTATTGTATTAatcata                                                                   >  2:522691/1‑149 (MQ=255)
                                                                      aCCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcatatca                                                                >  1:1140645/1‑149 (MQ=255)
                                                                          ttccATTCTTATCTAATTCTTGGTATGTTTTCATATCTGCTTGCATACCTACGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCAtttt                                                            <  2:1173924/147‑1 (MQ=255)
                                                                          tGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCAtttt                                                            >  1:882769/1‑149 (MQ=255)
                                                                          tGCCATACGTATCTAACTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCAtttt                                                            <  1:1667328/149‑1 (MQ=255)
                                                                            ccATACGTATCTAATTCTTGGAACGTTTGCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTa                                                          <  2:36038/149‑1 (MQ=255)
                                                                               tACGTATGTAATTTTTGGAAAGTTTTCAAATTTGTTTGGATATGTGCGTCAGCACTTTTGTAAGAACCGGTAATTGTTAAAACTTTCTTTAAAGGAATTGAGTTTTACTTGTACATAATAATATTGTGTTTATTATATTATTTTTaatt                                                       >  2:1560440/1‑147 (MQ=255)
                                                                                      cTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGa                                                <  1:407640/149‑1 (MQ=255)
                                                                                        aaTTCTTGGAACGTTATCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCAATGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTCTTAACCTAATTGaaa                                              <  1:1075816/149‑1 (MQ=255)
                                                                                             ttGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAgtat                                         >  1:847777/1‑149 (MQ=255)
                                                                                             ttGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAatat                                         >  2:1322110/1‑149 (MQ=255)
                                                                                              tGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTCTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATAtt                                        >  1:1697461/1‑149 (MQ=255)
                                                                                              tGGAACGTTTTCAAAGCTGCTTGCATACCGGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATAtt                                        <  2:1689378/149‑1 (MQ=255)
                                                                                               ggAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTa                                       >  1:878770/1‑149 (MQ=255)
                                                                                               ggAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAAATTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTa                                       >  2:1042116/1‑149 (MQ=255)
                                                                                                gagcGTTTTCAAATCTGCCTGCCTACCTGCGCCGGCACTTTTGTCAGTACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTTCTACATAATAATACTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTaa                                      <  2:991559/146‑1 (MQ=255)
                                                                                                gAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAAGCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAACTTCAATCCTACATAATAATTGTGTATTTTGCATTTCATTTTTAACCTAATTGAAAAATATTaa                                      <  2:340739/149‑1 (MQ=255)
                                                                                                 aaCGTTTTCAAATCTGCTTGCATACTTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAATTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTGATTGAAAAATATTAAg                                     <  2:679190/149‑1 (MQ=255)
                                                                                                 aaCGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAACAAAtttct                                     >  2:424656/1‑146 (MQ=255)
                                                                                                 aaCGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAg                                     >  2:1126069/1‑149 (MQ=255)
                                                                                                  aCGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGc                                    >  1:1606600/1‑149 (MQ=255)
                                                                                                    gTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTCAATGTGACTTTAACGAATTACAAGGTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCAt                                  >  1:827895/1‑149 (MQ=255)
                                                                                                    gTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCAt                                  >  1:788621/1‑149 (MQ=255)
                                                                                                    gTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAACCAt                                  >  2:598914/1‑149 (MQ=255)
                                                                                                     ttttCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCAtt                                 <  1:1387752/149‑1 (MQ=255)
                                                                                                          aaaTCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAata                            <  1:38265/149‑1 (MQ=255)
                                                                                                          aaaTCCGCTTGCACACCTGCGCCAGCACTTTCGTGAGATCCGGCAATTGGTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAata                            <  1:1326896/149‑1 (MQ=255)
                                                                                                              cTGCTTGCACACCTCCGCCAGCACTTTTGTCAGAACCGGCAAGTGTTAAAACTCTCTGTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTtg                        <  1:229134/149‑1 (MQ=255)
                                                                                                                  ttGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTGATCATATCATTTTTAACATCATTGAAAGATATTAAGCATTCCGTATTtgatgt                    >  1:1642030/1‑148 (MQ=255)
                                                                                                                  ttGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTtgatga                    >  2:398046/1‑149 (MQ=255)
                                                                                                                   tGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTgatgat                   >  1:529237/1‑149 (MQ=255)
                                                                                                                   tGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTgatgat                   >  1:484853/1‑149 (MQ=255)
                                                                                                                      aTACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgt                >  1:958034/1‑149 (MQ=255)
                                                                                                                      aTACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCAATCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTACTTGAAAAAAATTAAGCATTCAATATTTGATGTttgt                >  1:580612/1‑149 (MQ=255)
                                                                                                                        aCCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgttg              <  2:1395136/149‑1 (MQ=255)
                                                                                                                          cTGCGCCAGCACTTTTGTCAGACCCGGCAATTGTTAAAACTTACTGTAACGCCATTGAGCTTCTCTCCTACGTAATAAGATTGACTTCAACCTATCATTTTTAACCTAATTGAATCATATTAAGCATGCAAGATTTGATGATTGTTGaa            <  2:347417/149‑1 (MQ=255)
                                                                                                                           tGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGaaa           <  1:247318/149‑1 (MQ=255)
                                                                                                                           tGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGaaa           <  1:1486699/149‑1 (MQ=255)
                                                                                                                           tGCGCCAGCACTTTTGTCAGAACCGGAAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGaaa           <  2:316077/149‑1 (MQ=255)
                                                                                                                            gcgcCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATACTACTATGGTATTCATCAGATCATTTTTAACCTAATTAAAAAAGATTAAGCATTAAATATTTGATGATTGTTGAAAt          >  2:1040310/1‑149 (MQ=255)
                                                                                                                                  gCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATGAAtca    >  1:113843/1‑149 (MQ=255)
                                                                                                                                    aCTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATAGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATGAAtcatt  <  2:882769/149‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                
ATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAATGAATCATT  >  USA300TCH1516_ALE/638478‑638758

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: