Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,921,750 G→A 6.0% G353G (GGC→GGT menE ← 2‑succinylbenzoate‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,921,7500GA 6.0% 251.5 / 4.3 84G353G (GGC→GGTmenE2‑succinylbenzoate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base G (34/44);  new base A (3/2);  total (37/47)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

ATAATTAAATCTTTACGTCGGTCATAAATCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTGCTGTTAAAAATTGCGAACATGTCT  >  USA300TCH1516_ALE/1921610‑1921889
                                                                                                                                            |                                                                                                                                           
ataatTAAATCTTTACGTCGGTCATATATCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGgcg                                                                                                                                             <  2:858813/141‑1 (MQ=255)
ataatTAAATCTTTACGTCGGTCATAAATCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGgcg                                                                                                                                             >  2:255890/1‑141 (MQ=255)
  acttAAGTCTTTACGTCGGTCATAAAACATGACATAACCTTCATGATCTATTGCAGCAATGTCACCCGTATTAAAATAACCATTCTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCGTCACATTGGCGcc                                                                                                                                           <  4:112214/139‑1 (MQ=255)
  aatTAAATCTTTACGTCGGTCATAAATCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGcc                                                                                                                                           <  2:156295/141‑1 (MQ=255)
  aatTAAATCTTTACGTCGGTCATAAATCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGcc                                                                                                                                           <  2:1007826/141‑1 (MQ=255)
   atTAAATCTTTACGTCGGTCATAAATCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCt                                                                                                                                          >  2:579305/1‑141 (MQ=255)
    ttAAATCTTTACGTCGGTCATAAATCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCCGTTGGACACAAATATACATTCATCACATTGTCGCCtt                                                                                                                                         >  1:313643/1‑141 (MQ=255)
      aaaTCTTTACGTCGGTCATAAATCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCATTCGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTa                                                                                                                                       >  1:237374/1‑141 (MQ=255)
      aaaTCTTTACGTCGGTCATAAATCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTa                                                                                                                                       >  2:142000/1‑141 (MQ=255)
        aTCTTTACGTCGGTCATAAATCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAAt                                                                                                                                     <  2:936557/141‑1 (MQ=255)
           tttACGTCTGTCATAAATCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCAt                                                                                                                                  <  1:227909/141‑1 (MQ=255)
           tttACGTCGGTCATAAATCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCAt                                                                                                                                  <  4:307442/141‑1 (MQ=255)
             tACGTCGGTCATAAATCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTa                                                                                                                                >  1:380047/1‑141 (MQ=255)
              aCGTCGGTCATAAATCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTaa                                                                                                                               >  4:145254/1‑141 (MQ=255)
                 tCGGTCATAAATCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATtc                                                                                                                            >  1:943825/1‑141 (MQ=255)
                      cATAAATCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCAt                                                                                                                       <  2:1004235/141‑1 (MQ=255)
                           aTCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACct                                                                                                                  <  3:209369/141‑1 (MQ=255)
                           aTCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACct                                                                                                                  <  2:104867/141‑1 (MQ=255)
                           aTCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACct                                                                                                                  <  2:386550/141‑1 (MQ=255)
                             cATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACcttc                                                                                                                <  1:995716/141‑1 (MQ=255)
                              aTGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACcttct                                                                                                               <  1:243682/141‑1 (MQ=255)
                               tGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACcttctt                                                                                                              >  1:832742/1‑141 (MQ=255)
                                   aTAACCTTCATGATCTATTTCAGCAATGTCACCTGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATCCATTCATCACATTGGCACCTTTAATCATTAATTCTCCATGACCTTCTttatt                                                                                                          <  3:11156/141‑1 (MQ=255)
                                   aTAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTttatt                                                                                                          <  1:648349/141‑1 (MQ=255)
                                         ttCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGatt                                                                                                    <  1:676342/141‑1 (MQ=255)
                                          tCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGattt                                                                                                   <  1:771897/141‑1 (MQ=255)
                                          tCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGattt                                                                                                   >  3:275323/1‑141 (MQ=255)
                                            aTGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTTGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGattttt                                                                                                 >  3:263904/1‑141 (MQ=255)
                                            aTGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGattttt                                                                                                 <  3:101746/141‑1 (MQ=255)
                                                 cTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAGATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGatttttaattt                                                                                            >  1:843210/1‑141 (MQ=255)
                                                 cTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGatttttaattt                                                                                            >  1:895032/1‑141 (MQ=255)
                                                     ttCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTAc                                                                                        <  1:395474/141‑1 (MQ=255)
                                                           aaTGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAAc                                                                                  >  1:925610/1‑141 (MQ=255)
                                                            aTGTCATCCGTATGAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACa                                                                                 <  2:39846/141‑1 (MQ=255)
                                                            aTGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACa                                                                                 <  4:16660/141‑1 (MQ=255)
                                                             tGTCACCTGCATTAAATTAACCATTTTCAAATGTACTTGTTTTTTTTGTTGGATACAAATATCCATTCATCACATTGGCACCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTAACGTCAACAt                                                                                >  4:169491/1‑141 (MQ=255)
                                                             tGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTTAACAt                                                                                >  1:579891/1‑141 (MQ=255)
                                                             tGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACAt                                                                                <  2:354996/141‑1 (MQ=255)
                                                                 aCCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGc                                                                            <  1:694234/141‑1 (MQ=255)
                                                                     gTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTGTATAAGGATTTTTAATTTGTACGTCAACCtttagcgaa                                                                        <  2:661916/141‑8 (MQ=255)
                                                                       aTTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTgg                                                                      <  2:422221/141‑1 (MQ=255)
                                                                        ttAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGc                                                                     <  2:260385/141‑1 (MQ=255)
                                                                         tAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCa                                                                    >  3:62258/1‑141 (MQ=255)
                                                                          aaaaTAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCAt                                                                   >  4:336091/1‑141 (MQ=255)
                                                                           aaaTAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATc                                                                  >  4:101599/1‑141 (MQ=255)
                                                                            aaTAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACATCAACAATGGCACTTCGCATcc                                                                 >  3:47101/1‑141 (MQ=255)
                                                                              tagccATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCt                                                               <  2:959360/138‑1 (MQ=255)
                                                                                aCCATTTTCAAACGTGCCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTAc                                                             <  1:899807/141‑1 (MQ=255)
                                                                                  cATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAg                                                           <  2:752082/141‑1 (MQ=255)
                                                                                    ttttCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAgtt                                                         >  2:1053776/1‑140 (MQ=255)
                                                                                    ttttCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAgtg                                                         >  2:510703/1‑141 (MQ=255)
                                                                                    ttttCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAgtg                                                         >  1:441272/1‑141 (MQ=255)
                                                                                         aaaCGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAgg                                                    <  2:307141/141‑1 (MQ=255)
                                                                                          aaCGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGa                                                   >  3:156113/1‑141 (MQ=255)
                                                                                          aaCGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGa                                                   >  3:48863/1‑141 (MQ=255)
                                                                                            cGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACg                                                 >  1:788167/1‑141 (MQ=255)
                                                                                                  cGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGtgcatg                                           <  1:278168/141‑1 (MQ=255)
                                                                                                     tAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCaa                                        <  1:77173/141‑1 (MQ=255)
                                                                                                       aaTCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACa                                      <  4:62258/141‑1 (MQ=255)
                                                                                                        aTCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACAt                                     <  2:515681/141‑1 (MQ=255)
                                                                                                          cTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACAttt                                   >  3:138009/1‑141 (MQ=255)
                                                                                                          cTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACAttt                                   >  2:539486/1‑141 (MQ=255)
                                                                                                          cTGTTGGATACAAATATACATTCAGCACATTGTCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCTCTACATTGTCAGGACGTGCATGCAACtttt                                   >  2:651621/1‑141 (MQ=255)
                                                                                                           tGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTc                                  >  2:122803/1‑141 (MQ=255)
                                                                                                            gTTGGATACAAATATACATTCATCTCATTGGCTCCTTTGATCATTATTTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTcc                                 <  3:89131/141‑1 (MQ=255)
                                                                                                             ttGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCg                                <  1:331085/141‑1 (MQ=255)
                                                                                                                gATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGgtg                             <  3:354809/141‑1 (MQ=255)
                                                                                                                 aTACAAATATCCATTCATCACATTGGCACCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTAACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGgtgt                            >  2:940438/1‑141 (MQ=255)
                                                                                                                  tACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGtt                           <  1:591216/141‑1 (MQ=255)
                                                                                                                   aCAAATATACATTCATCACGTTGGCGCCTTTGATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTg                          <  2:199848/141‑1 (MQ=255)
                                                                                                                    cAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGATTTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTGc                         >  2:606961/1‑141 (MQ=255)
                                                                                                                        tatccATTCATCACATTGGCACCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTAACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTGCTGtt                     <  3:169491/137‑1 (MQ=255)
                                                                                                                        tataCATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTGCTGtt                     >  3:238143/1‑141 (MQ=255)
                                                                                                                        tataCATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTACTAGGATTTTTAACTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTGCTGtt                     <  1:487347/141‑1 (MQ=255)
                                                                                                                         ataCATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTGCTGTTa                    <  2:198006/141‑1 (MQ=255)
                                                                                                                           actttCATCACATTCGCGCCTTTAATCATTAATTCTCCATTACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTGCTGTTaaa                  <  4:371013/138‑1 (MQ=255)
                                                                                                                           aCATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTGCTGTTaaa                  <  2:460488/141‑1 (MQ=255)
                                                                                                                           aCATTCATCACATTGGCGCATTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTGCTGTTaaa                  <  2:562294/141‑1 (MQ=255)
                                                                                                                              ttCATCACATTGGCACCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTAACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTGCTGTTAAAAAt               >  1:250838/1‑141 (MQ=255)
                                                                                                                                    acaTTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTGCTGTTAAAAATTGCGaa         <  2:216661/141‑1 (MQ=255)
                                                                                                                                    acaTTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTGCTGTTAAAAATTGCGaa         <  3:328134/141‑1 (MQ=255)
                                                                                                                                      aTTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTGCTGTTAAAAATTGCGAACa       >  2:572058/1‑141 (MQ=255)
                                                                                                                                      aTTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTGCTGTTAAAAATTGCGAACa       >  1:5283/1‑141 (MQ=255)
                                                                                                                                      aTTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTGCTGTTAAAAATTGCGAACa       >  4:55584/1‑141 (MQ=255)
                                                                                                                                       ttGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTGCTGTTAAAAATTGCGAACAt      >  1:829724/1‑141 (MQ=255)
                                                                                                                                         ggCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTGCTGTTAAAAATTGCGAACATGt    <  4:139908/141‑1 (MQ=255)
                                                                                                                                           cgcCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTGCTGTTAAAAATTGCGAACATGtct  >  3:89951/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                           
ATAATTAAATCTTTACGTCGGTCATAAATCATGACATAACCTTCATGATCTATTTCAGCAATGTCACCCGTATTAAAATAACCATTTTCAAACGTACCCGTTAAATCTGTTGGATACAAATATACATTCATCACATTGGCGCCTTTAATCATTAATTCTCCATGACCTTCTTTATTAGGATTTTTAATTTTTACGTCAACATTGGCACTTGGCATCCCTACAGTGTCAGGACGTGCATGCAACATTTCCGGTGTTGCTGTTAAAAATTGCGAACATGTCT  >  USA300TCH1516_ALE/1921610‑1921889

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