Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,390,270 T→C 7.3% E54G (GAA→GGA)  moaA ← GTP 3',8‑cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,390,2700TC 7.3% 185.4 / 3.8 65E54G (GAA→GGA) moaAGTP 3',8‑cyclase
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base C (2/3);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TATCTTCAATACCATCGATTTGATTTAATTTAGCTATAAGTACATCTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTGTCACAGATAACCGTAAGTCACGGATGGGACGTCCTAGT  >  USA300TCH1516_ALE/2390136‑2390407
                                                                                                                                      |                                                                                                                                         
tatCTTCAATACCATCGATTTGATTTAATTTAGCTATAAGTACATCTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATTaa                                                                                                                                     >  2:140805/1‑141 (MQ=255)
 atCTTCAATACCATCGATTTGATTTAATTTAGCTATAAGTACATCTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCaaa                                                                                                                                    >  2:491191/1‑141 (MQ=255)
 atCTTCAATACCATCGATTTGATTTAATTTAGCTATAAGTACATCTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCaaa                                                                                                                                    <  2:441491/141‑1 (MQ=255)
    tttaaTACCTTCGATTTGATTTAATTTAGCTATAAGTACATCTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACGCCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGt                                                                                                                                 <  2:25200/138‑1 (MQ=255)
    ttCAATACCTTCGATTTGATTTAATTTAGCTATAAGTACATCTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGt                                                                                                                                 <  1:277434/141‑1 (MQ=255)
           ccATCGATTTGATTTAATTTAGCTATAAGTACATCTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTCCACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTCCATCACACGTTAAAAGt                                                                                                                          >  2:608207/1‑141 (MQ=255)
           ccATCGATTTGATTTAATTTAGCTATAAGTACATCTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGt                                                                                                                          >  2:241673/1‑141 (MQ=255)
           ccATCGATTTGATTTAATTTAGCTATAAGTACATCTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTCCATCAAACGTTAAAAGt                                                                                                                          >  2:515500/1‑141 (MQ=255)
            cctCGATTTGATTTAATTTAGCTATAAGTACATCTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGtt                                                                                                                         <  1:73012/139‑1 (MQ=255)
                  tttgatttAATTTAGCTATAAGTACATCTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTCAAAGTTCAtttt                                                                                                                   >  1:375398/1‑141 (MQ=255)
                  tttgatttAATTTAGCTATAAGTACATCTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCAtttt                                                                                                                   >  1:535691/1‑141 (MQ=255)
                     gatttAATTTAGCTATAAGTACATCTACATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACTTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCAGTTTTAg                                                                                                                >  1:756822/1‑141 (MQ=255)
                     gatttAATTTAGCTATAAGTACATCTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTGGCGGTTCTAGCCATTTCATCAACCGTTAAAAGTTCTTTTTTAg                                                                                                                >  2:469443/1‑141 (MQ=255)
                       tttaatttaGCTATAAGTACATCTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGt                                                                                                              >  1:581954/1‑141 (MQ=255)
                          aatttaGCTATAAGTACATTTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTaaa                                                                                                           <  1:477572/141‑1 (MQ=255)
                          aatttaGCTATAAGTACATCTAAATCCCGGCGCATCAATGGTTCTCCCCCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTaaa                                                                                                           <  1:438318/141‑1 (MQ=255)
                          aagttaGCTATCAGTACATCTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTaaa                                                                                                           <  1:254228/138‑1 (MQ=255)
                           atttaGCTATAAGTACATCTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTaaaa                                                                                                          >  1:179933/1‑141 (MQ=255)
                            tttaGCTATAAGTACATCTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCTATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTaaaaa                                                                                                         <  2:261878/141‑1 (MQ=255)
                                             cTAAATCCCGTCGCATCAATGTTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCaa                                                                                        <  2:437574/141‑1 (MQ=255)
                                             cTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCaa                                                                                        <  1:685723/141‑1 (MQ=255)
                                                aaTCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATa                                                                                     >  1:37071/1‑141 (MQ=255)
                                                aaTCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTA‑TTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTACAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATAc                                                                                    >  2:167750/1‑141 (MQ=255)
                                                  tCCCGTCGCATGAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAGACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATAcc                                                                                   <  2:367154/141‑1 (MQ=255)
                                                  tCCCGTCGCATCAACTGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACCTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATAcc                                                                                   <  2:496722/141‑1 (MQ=255)
                                                     cGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACctct                                                                                <  1:166642/141‑1 (MQ=255)
                                                     cGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACctct                                                                                >  1:526864/1‑141 (MQ=255)
                                                      gTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCtt                                                                               <  1:449921/141‑1 (MQ=255)
                                                         gCATCAATGGTTCTCCACCTGTGATGAGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATGTCATCAGACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCCAATCCCTCtttgc                                                                            >  1:742876/1‑139 (MQ=255)
                                                            acAATGGTTCTCCACCTGTAATGCGTGTTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCa                                                                         <  1:464262/140‑1 (MQ=255)
                                                                tGGTTCTCCCCCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCa                                                                     <  1:660970/141‑1 (MQ=255)
                                                                 ggTTCTCCACCTGTAATGCGTCTTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCaa                                                                    <  2:398929/141‑1 (MQ=255)
                                                                 ggTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCaa                                                                    >  1:660938/1‑141 (MQ=255)
                                                                   ttCTCCACCTGTAATGCGTATTTTTTTTTTTTTTTATTCTGCATATACCTTAACGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAACCTACTACAATGACCCATTTCGGCATGCaatt                                                                  >  2:766516/1‑140 (MQ=255)
                                                                   ttCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCACTTTAGGCATGCaata                                                                  >  1:29117/1‑141 (MQ=255)
                                                                       ccACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAAAGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCCCTTTAGGCATGCACTAATcc                                                              >  2:273147/1‑140 (MQ=255)
                                                                        cACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATAAACCTTTAGGCATGAAATAATcac                                                             >  2:677136/1‑141 (MQ=255)
                                                                          ccTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCTATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAAAAAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCacac                                                           <  2:276452/141‑1 (MQ=255)
                                                                          ccTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCacac                                                           <  1:459738/141‑1 (MQ=255)
                                                                             gTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTa                                                        <  1:426730/141‑1 (MQ=255)
                                                                                 tGCGTATTTTTTTTACACCTAATTCTGAATATACCTTCGCGATTCTAGCCATTTCCTCAAACGTTAAAAGTTCATTTTTAGGCAAAAATACGAAATCATATCCACATGCCACTTTAGACATGCACTATACACACCTaaact                                                    >  2:750116/1‑139 (MQ=255)
                                                                                  gCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGtt                                                   >  1:49285/1‑141 (MQ=255)
                                                                                              tACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTGt                                       <  2:751806/141‑1 (MQ=255)
                                                                                                     aaTTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTACAGTTACACCGATCTGTCACAGAt                                >  2:607536/1‑141 (MQ=255)
                                                                                                      agtcggCATATACCTTAGCGATTCTAGCCATTCCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTGTCACAGATa                               <  1:608207/136‑1 (MQ=255)
                                                                                                      aTTCTGCATATACCTTAGCGATTCTAGCCATTCCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTGTCACAGATa                               <  1:515500/141‑1 (MQ=255)
                                                                                                          tccATATACCTTAGCGATACTAGCCATTGCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTGTCACAGATAACCg                           <  2:434602/139‑1 (MQ=255)
                                                                                                          tGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTCCACCGATCTGTCACAGATAACCg                           <  2:508474/141‑1 (MQ=255)
                                                                                                           gCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTTTCACAAATAACCga                          >  1:344408/1‑140 (MQ=255)
                                                                                                            cATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTGAAGTTGCACCGACCTGCAACAGATAACTGTa                         >  1:553040/1‑141 (MQ=255)
                                                                                                               ataCCTTCGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTGTCACAGATAACCGGAAGt                      <  1:521713/141‑1 (MQ=255)
                                                                                                               acaCCTTAGCGATTCTAGCCAGTCCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTGTCACAGATAACCGTAAGt                      <  2:410619/139‑1 (MQ=255)
                                                                                                                     tAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTGTCACAGATAACCGTAAGTCACGGa                >  1:134080/1‑141 (MQ=255)
                                                                                                                     tAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTGTCACAGATAACCGTAAGTCACGGa                <  1:491191/141‑1 (MQ=255)
                                                                                                                      aGCGATTCTAGCCATTTCAACAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAAAAATCTCAAAATACCTCTTTAGACAAGCAATAATCACACCTAAAGTAACACCGATCTGTCACAGATAACCGTAAGTCACGAAt               >  2:523145/1‑141 (MQ=255)
                                                                                                                         gATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTGTCACAGATAACCGTAAGTCACGGATggg            <  1:44023/141‑1 (MQ=255)
                                                                                                                         gATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTGTCACAGATAACCGTAAGTCACGGATggg            <  2:746334/141‑1 (MQ=255)
                                                                                                                          aTTCTAGCCATTTCATCAAACGTTAAAATTTCATTTCGAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTGTCACAGAAAACCGTAAGTCACGGATGGGa           <  2:526864/141‑1 (MQ=255)
                                                                                                                            tCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTGTCACAGATAACCGTAAGTCACGGATGGGACg         <  2:581954/141‑1 (MQ=255)
                                                                                                                             cTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTGTCACAGATAACCGTACGTCACGGATGGGACgc        >  2:257989/1‑140 (MQ=255)
                                                                                                                               aGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTGTCACAGATAACCGTAAGTCACGGATGGGACGTcc      <  1:241673/141‑1 (MQ=255)
                                                                                                                                gCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGGCAGAAGCACACCTAAAGTTACACCGATCTGTCACAGATAACCGTAAGTCACGGATGGGACGTCCt     >  1:751441/1‑141 (MQ=255)
                                                                                                                                 ccATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTGTCACAGATAACCGTAAGTCACGGATGGGACGTCCTa    >  1:576305/1‑141 (MQ=255)
                                                                                                                                   aattCATCAAATGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCCCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTGTCACAGATAACCGTAAGTCACGGATGGGACGTCCTAGt  <  2:668881/139‑1 (MQ=255)
                                                                                                                                   aTTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAACGTTACACCGATCTGTCACAGATAACCGTAAGTCACGGATGGGACGTCCTAGt  >  1:299090/1‑141 (MQ=255)
                                                                                                                                   aTTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTGTCACAGATAACCGTAAGTCACGGATGGGACGTCCTAGt  >  2:643771/1‑141 (MQ=255)
                                                                                                                                      |                                                                                                                                         
TATCTTCAATACCATCGATTTGATTTAATTTAGCTATAAGTACATCTAAATCCCGTCGCATCAATGGTTCTCCACCTGTAATGCGTATTTTTTTTACACCTAATTCTGCATATACCTTAGCGATTCTAGCCATTTCATCAAACGTTAAAAGTTCATTTTTAGGTAAAAATACGAAATCATCTCCAAATACCTCTTTAGGCATGCAATAATCACACCTAAAGTTACACCGATCTGTCACAGATAACCGTAAGTCACGGATGGGACGTCCTAGT  >  USA300TCH1516_ALE/2390136‑2390407

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: ‑