Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 253,983 A→G 6.7% intergenic (+196/+166) asbF → / ← USA300HOU_RS01125 3‑dehydroshikimate dehydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE253,9830AG 6.7% 399.6 / 13.7 134intergenic (+196/+166)asbF/USA300HOU_RS011253‑dehydroshikimate dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (64/60);  new base G (2/7);  total (67/67)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 141 total aligned reads displayed.

GAATAAATGCATATGTCGCCAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATGGACGTCGATGTAGA  >  USA300TCH1516_ALE/253842‑254121
                                                                                                                                             |                                                                                                                                           
gAATAAATGCATATGTCGCCAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAaga                                                                                                                                             <  2:881639/141‑1 (MQ=255)
 aaTAAATGCATATGTCGCCAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAaga                                                                                                                                             <  1:374534/141‑1 (MQ=255)
   tAAATGCATAAGTCGCCAAGCCATTGCCAAAAATTTCACACCTTAAATAAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTCTATAGCGATGGTCATACCCATTCACAGTAACAATACTCACCATTGACAAGAGt                                                                                                                                           >  2:59247/1‑141 (MQ=255)
      aTGCATATGTCGCCAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtata                                                                                                                                        <  3:268108/141‑1 (MQ=255)
         cATATGTCGCCAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATaa                                                                                                                                     >  3:65614/1‑141 (MQ=255)
         aatatGTCGCCAAGGCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCCATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATaa                                                                                                                                     <  2:1120569/140‑1 (MQ=255)
          atatGTCGCCAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAc                                                                                                                                    >  3:134157/1‑141 (MQ=255)
           tatGTCGCCAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAcc                                                                                                                                   <  1:216703/141‑1 (MQ=255)
           tatGTCGCCAAGCCATGGCCCAAAATTTCACACCTTAAATCCAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGGTtatatactc                                                                                                                                   <  4:501881/141‑4 (MQ=255)
               tCGCCAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtttt                                                                                                                               <  1:585871/141‑1 (MQ=255)
               tCGCCAAGACATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCGCACCATTGAAAAGAGTATATAACCtttt                                                                                                                               <  4:36941/141‑1 (MQ=255)
                 gCCAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCa                                                                                                                             <  4:339309/141‑1 (MQ=255)
                 gCCAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCa                                                                                                                             >  3:482455/1‑141 (MQ=255)
                  ccAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCaa                                                                                                                            <  2:1015257/141‑1 (MQ=39)
                     aGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCAAtgg                                                                                                                         <  4:319542/141‑3 (MQ=38)
                       ccATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTg                                                                                                                       >  2:329166/1‑141 (MQ=255)
                        cATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGa                                                                                                                      <  2:863017/141‑1 (MQ=255)
                        cATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGa                                                                                                                      <  2:995981/141‑1 (MQ=255)
                         aTTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAg                                                                                                                     >  2:88209/1‑141 (MQ=255)
                          ttGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAgg                                                                                                                    <  2:924286/141‑1 (MQ=255)
                             ccAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGtat                                                                                                                 <  1:12874/141‑1 (MQ=255)
                             ccAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGtat                                                                                                                 <  1:881785/141‑1 (MQ=255)
                               aaaaaTTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTatat                                                                                                               >  2:1002931/1‑141 (MQ=255)
                                  aaTTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGat                                                                                                            <  3:116295/141‑1 (MQ=255)
                                    tttCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataa                                                                                                          <  1:114902/141‑1 (MQ=255)
                                       cacaCCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACTCATTCACAGGAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGAtaataa                                                                                                      >  3:400312/1‑141 (MQ=255)
                                       cacaCCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGAtaataa                                                                                                      >  2:418813/1‑141 (MQ=255)
                                        gcacCTTAAATTAAGTCATTGGTTGTAAATAATGTTTAATTTATAAAATTATATAGCGATTTTCATACCCATGCACAGTAACAATCCTCACCATTGAAAAGAGTATATAATCTTTTTTATTTTTAGTTATATATTAATaaa                                                                                                     >  4:427243/2‑141 (MQ=255)
                                        acacCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATATTaaa                                                                                                     >  1:1028267/1‑141 (MQ=255)
                                        acacCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaa                                                                                                     >  1:1023478/1‑141 (MQ=255)
                                              tAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAg                                                                                               <  1:968865/141‑1 (MQ=255)
                                                aaTCAAGACATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGcc                                                                                             <  2:180618/141‑1 (MQ=255)
                                                  tCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtg                                                                                           <  3:29418/141‑1 (MQ=255)
                                                      gTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttgt                                                                                       <  1:622918/141‑1 (MQ=255)
                                                      gTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttgt                                                                                       >  2:859191/1‑141 (MQ=255)
                                                          ttgtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTcaca                                                                                   <  3:502538/141‑1 (MQ=255)
                                                             tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATg                                                                                >  1:454771/1‑141 (MQ=255)
                                                              gagaaaaGAAGGTGTACTTTCTATAAGTATCTAGAGATCGTCATACCCATTCAGAGAAACAAACCTCACCAATGAAAAGAGTATATAACCTTTTCGATAGGAGGGTATATGAAAAAAAAAAAAAGCCTGTTGTCACAATgg                                                                               <  2:1037057/137‑1 (MQ=255)
                                                                gtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCGTTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTTTCACAATGGTc                                                                             >  4:162900/1‑141 (MQ=255)
                                                                 tAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                                                                            <  2:425216/141‑1 (MQ=255)
                                                                  aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAACGAGTCAATCACATTATCAATAGTAAAGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                                                                           >  1:174634/1‑141 (MQ=255)
                                                                        ggTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGacac                                                                     >  3:394504/1‑141 (MQ=255)
                                                                         gtgtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACg                                                                    <  1:212201/141‑1 (MQ=255)
                                                                            tACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACa                                                                 <  1:476149/141‑1 (MQ=255)
                                                                            tACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACa                                                                 <  1:51713/141‑1 (MQ=255)
                                                                                 tatataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACtt                                                            >  3:132100/1‑141 (MQ=255)
                                                                                   tataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTa                                                          <  1:191843/141‑1 (MQ=255)
                                                                                   tataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTa                                                          >  1:435974/1‑141 (MQ=255)
                                                                                       aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCaa                                                                                                                            >  2:534849/1‑72 (MQ=255)
                                                                                       aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCaa                                                                                                                            <  1:534849/72‑1 (MQ=255)
                                                                                       aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAgg                                                      >  1:303321/1‑141 (MQ=255)
                                                                                       aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAgg                                                      <  1:320364/141‑1 (MQ=255)
                                                                                           tataGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAATAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTc                                                  <  4:394563/141‑1 (MQ=255)
                                                                                           tataGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTc                                                  <  4:394485/141‑1 (MQ=255)
                                                                                            ataGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCt                                                 >  1:91273/1‑141 (MQ=255)
                                                                                             taGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTCTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTg                                                >  1:936703/1‑141 (MQ=255)
                                                                                                 gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaata                                            >  4:34958/1‑141 (MQ=255)
                                                                                                 gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaata                                            >  2:860093/1‑141 (MQ=255)
                                                                                                 gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGAATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaata                                            >  3:417978/1‑141 (MQ=255)
                                                                                                 gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGAATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCAGTTGTCACAATGGTCATAGACACTACATACTTTAAAGTTTTCTGatta                                            >  3:417970/1‑141 (MQ=255)
                                                                                                  aTGGTCATACCCATTAACACAAAAAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAAAAAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaata                                            <  4:176160/141‑1 (MQ=255)
                                                                                                    ggTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatataa                                         <  4:235771/141‑1 (MQ=255)
                                                                                                     gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatataat                                        >  1:213399/1‑141 (MQ=255)
                                                                                                     gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaagataat                                        >  1:842695/1‑141 (MQ=255)
                                                                                                      tCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTTTAaatatacta                                       >  1:1004918/1‑141 (MQ=255)
                                                                                                       cATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatataatat                                      <  3:395107/141‑1 (MQ=255)
                                                                                                        aTACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATAtt                                     >  2:686513/1‑141 (MQ=255)
                                                                                                               ttCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAAt                              <  3:298558/141‑1 (MQ=255)
                                                                                                                tCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATg                             <  1:568373/141‑1 (MQ=255)
                                                                                                                acacaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGGATATAATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATg                             <  2:1037141/140‑1 (MQ=255)
                                                                                                                 cacaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGc                            >  4:598/1‑141 (MQ=255)
                                                                                                                 cacaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATATCAGAATGc                            >  2:298228/1‑141 (MQ=255)
                                                                                                                    aGTAACAATCCTCACTATTGAAAAGAGTATATAACCTTTTCAATAGTGACGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACAAGACATACTTTCAAGATTTCTGAATATAATATTTCAGAATGCACt                         >  4:242045/1‑141 (MQ=255)
                                                                                                                     gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATCGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACtt                        >  2:451730/1‑141 (MQ=255)
                                                                                                                     gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACtt                        >  3:273854/1‑141 (MQ=255)
                                                                                                                     gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACtt                        >  4:370000/1‑141 (MQ=255)
                                                                                                                       aaCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCCGTAGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTAAATATAATATTTCAGAATGCACTTTa                      >  1:416059/1‑141 (MQ=255)
                                                                                                                        aCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATgg                                                                               <  2:263325/83‑1 (MQ=255)
                                                                                                                        aCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATCGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGACTGCACTTTaa                     >  2:1113044/1‑141 (MQ=255)
                                                                                                                         aactgcTCGCCATTGAAAAGAGTATATAACCTTTTAAATAGTGAAGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTaaa                    <  4:400312/136‑1 (MQ=255)
                                                                                                                           aTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGa                  >  1:397964/1‑141 (MQ=255)
                                                                                                                           aTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGa                  >  1:874740/1‑141 (MQ=255)
                                                                                                                             ccTCACCATTGAAAAGATTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATg                <  2:1180690/141‑1 (MQ=255)
                                                                                                                                cACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAGAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATGGAc             <  1:860093/141‑1 (MQ=255)
                                                                                                                                 aCCATTGAAAAGAGTATATAACCTTATCAATAGTGAGGCATATGATAATAGAAAAAAGCCCGTTGCCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATACTTCAGAATGCACTTTAAAGATGGACg            <  1:391562/141‑1 (MQ=255)
                                                                                                                                   cATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATGGACGTc          <  1:1153432/141‑1 (MQ=255)
                                                                                                                                    aTTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATGGACGTCg         >  1:89742/1‑141 (MQ=255)
                                                                                                                                    aTTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATGGACGTCg         >  1:14435/1‑141 (MQ=255)
                                                                                                                                     ttGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATGGACGTCGa        <  2:803724/141‑1 (MQ=255)
                                                                                                                                     ttGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATGGACGTCGa        <  2:803779/141‑1 (MQ=255)
                                                                                                                                      tGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATGGACGTCGAt       >  2:1033709/1‑141 (MQ=255)
                                                                                                                                       gAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGa                  >  1:733486/1‑129 (MQ=39)
                                                                                                                                       gAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGa                  <  2:733486/129‑1 (MQ=39)
                                                                                                                                        aaaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATGGACGTCGATGt     >  2:1109983/1‑141 (MQ=255)
                                                                                                                                        aaaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATGGACGTCGATGt     >  4:56990/1‑141 (MQ=255)
                                                                                                                                        aaaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATGGACGTCGATGt     <  3:267377/141‑1 (MQ=255)
                                                                                                                                          aaGAGTATATAACCTTTTCAATAGTGAGGTATATGATACTAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATGGACGTCGATGTAg   >  1:632102/1‑141 (MQ=255)
                                                                                                                                           agagTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATGGACGTCGATGTAGa  <  1:785186/141‑1 (MQ=255)
                                                                                                                                           agagTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATGGACGTCGATGTAGa  >  1:1188280/1‑141 (MQ=255)
                                                                                                                                             agTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTa                                                          <  1:660829/83‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                           
GAATAAATGCATATGTCGCCAAGCCATTGCCAAAAATTTCACACCTTAAATCAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATAATATTTCAGAATGCACTTTAAAGATGGACGTCGATGTAGA  >  USA300TCH1516_ALE/253842‑254121

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