Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 370,891 C→A 5.9% P205P (CCC→CCA mhpB → 2,3‑dihydroxyphenylpropionate 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913370,8910CA 5.9% 264.0 / 5.0 84P205P (CCC→CCAmhpB2,3‑dihydroxyphenylpropionate 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (36/43);  new base A (3/2);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

GGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGG  >  NC_000913/370751‑371036
                                                                                                                                            |                                                                                                                                                 
ggACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGt                                                                                                                                             >  8:27395/1‑147 (MQ=255)
 gACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAg                                                                                                                                          >  8:302294/1‑149 (MQ=255)
 gACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAg                                                                                                                                          >  3:312811/1‑149 (MQ=255)
 gACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCAGCCAGTGAg                                                                                                                                          >  6:72342/1‑149 (MQ=255)
        cACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGaa                                                                                                                                    <  8:33595/148‑1 (MQ=255)
         aCCAGCACTGTCAATTAACGCGGGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGTCCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAAtt                                                                                                                                  <  2:78462/149‑1 (MQ=255)
           cAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGc                                                                                                                                <  8:236546/149‑1 (MQ=255)
           cAGCACTCTCAATAAACGCGTGCTGTTCCTGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCg                                                                                                                               >  3:20534/1‑149 (MQ=255)
               aCTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCa                                                                                                                            >  4:277533/1‑149 (MQ=255)
                ctctCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGt                                                                                                                              >  2:355388/1‑146 (MQ=255)
                ctctCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCa                                                                                                                            <  1:355388/148‑1 (MQ=255)
                 tctcAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAAtt                                                                                                                                  <  8:184466/141‑1 (MQ=255)
                 tctcAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAAtt                                                                                                                                  >  7:184466/1‑141 (MQ=255)
                   tcAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCaa                                                                                                                        >  2:135436/1‑149 (MQ=255)
                     aaTAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTACCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACg                                                                                                                      >  3:136466/1‑149 (MQ=255)
                       tAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACggg                                                                                                                    <  8:71143/149‑1 (MQ=255)
                         aaCGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGt                                                                                                                   <  3:277533/148‑1 (MQ=255)
                           cgcgTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGTGTGAt                                                                                                                <  4:312811/149‑1 (MQ=255)
                           cgcgTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGAt                                                                                                                >  8:169979/1‑149 (MQ=255)
                            gcgTGCTGTTCCTGGGTTCCGGTGGGCTTTGCCATCAGCCGCCGGTGCCCGAACTTGCTAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGTAAGATTTACCCTCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGAt                                                                                                                <  3:127320/148‑1 (MQ=255)
                             cgTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTa                                                                                                              <  7:169979/149‑1 (MQ=255)
                             cgTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCAGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTa                                                                                                              <  5:72342/149‑1 (MQ=255)
                              gTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAg                                                                                                             >  3:275297/1‑149 (MQ=255)
                               tGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAgc                                                                                                            >  5:210491/1‑149 (MQ=255)
                                    ttCCTGGGTTCCGGTGGGCTTTCCCATCAGACGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGcgccgc                                                                                                       >  6:109064/1‑149 (MQ=255)
                                          ggTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGaa                                                                                                 <  5:55747/149‑1 (MQ=255)
                                            ttCCGGTGGGCTTTCCCATCAGCCGCCGGTACCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGt                                                                                               <  4:136466/149‑1 (MQ=255)
                                               cGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGttt                                                                                             <  5:87093/148‑1 (MQ=255)
                                                ggTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttgt                                                                                           <  4:157254/149‑1 (MQ=255)
                                                    ggCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAg                                                                                       <  5:38108/149‑1 (MQ=255)
                                                    ggCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAg                                                                                       >  2:61544/1‑149 (MQ=255)
                                                           ccATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCaga                                                                                <  7:27395/149‑1 (MQ=255)
                                                            cATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCa                                                                                  >  6:364001/1‑146 (MQ=255)
                                                            cATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCCCCAGTCAGCGCGAATTCCTAAAGAAACGGGGCATTAGCGCCGCTGAGAAGTTTGTTGAGGATCaga                                                                                >  7:128189/1‑148 (MQ=255)
                                                             aTCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGTATCagag                                                                               >  4:53755/1‑148 (MQ=255)
                                                                aGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAAc                                                                            <  5:109064/148‑1 (MQ=255)
                                                                aGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACg                                                                           >  1:138926/1‑149 (MQ=255)
                                                                aGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACg                                                                           >  6:187827/1‑149 (MQ=255)
                                                                   cgccgGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCAGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGttt                                                                                             >  6:74610/1‑128 (MQ=255)
                                                                   cgccgGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCAGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGttt                                                                                             <  5:74610/128‑1 (MQ=255)
                                                                    gccgGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGaa                                                                                                 <  1:239107/123‑1 (MQ=255)
                                                                    gccgGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGaa                                                                                                 >  2:239107/1‑123 (MQ=255)
                                                                            ccGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCt                                                                >  4:339240/1‑148 (MQ=255)
                                                                            ccGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCt                                                                >  4:344521/1‑148 (MQ=255)
                                                                            ccGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGGTTGTTGAGGATCAGAGAACGCTGc                                                                       <  6:54016/141‑1 (MQ=255)
                                                                            ccGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGGTTGTTGAGGATCAGAGAACGCTGc                                                                       >  5:54016/1‑141 (MQ=255)
                                                                               aaCTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAATGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAAc                                                            <  4:20534/149‑1 (MQ=255)
                                                                                  tGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAAccc                                                          >  1:138770/1‑148 (MQ=255)
                                                                                      gAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACggg                                                                                                                    <  8:252366/86‑1 (MQ=255)
                                                                                      gAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACggg                                                                                                                    >  7:252366/1‑86 (MQ=255)
                                                                                      gAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTAGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGAttt                                                     >  2:164827/1‑149 (MQ=255)
                                                                                        aaGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCg                                                                  >  2:308912/1‑134 (MQ=255)
                                                                                        aaGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCg                                                                  <  1:308912/134‑1 (MQ=255)
                                                                                        aaGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTgg                                                   <  1:135436/149‑1 (MQ=255)
                                                                                         aGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTggg                                                  <  1:61544/149‑1 (MQ=255)
                                                                                         aGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTggg                                                  <  1:95771/149‑1 (MQ=255)
                                                                                           ccGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGAt                                                <  6:210491/149‑1 (MQ=255)
                                                                                             gATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATaa                                              >  5:45563/1‑149 (MQ=255)
                                                                                             gATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATaa                                              <  3:344521/149‑1 (MQ=255)
                                                                                             gATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCGGCATCCGCTCAACCCGATTTGGGATaa                                              <  3:339240/149‑1 (MQ=255)
                                                                                              aTGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAAc                                             <  7:302294/149‑1 (MQ=255)
                                                                                                  ccATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttgtt                                                                                          <  5:370957/100‑1 (MQ=255)
                                                                                                  ccATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttgtt                                                                                          >  6:370957/1‑100 (MQ=255)
                                                                                                        gcgcGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATg                                    <  1:44986/148‑1 (MQ=255)
                                                                                                        gcgcGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATg                                    >  2:44986/1‑148 (MQ=255)
                                                                                                         cgcgACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGa                                   <  5:364001/148‑1 (MQ=255)
                                                                                                            taCCGTCTGTTGGGGAGCGGGAATGATTTACCCGCCAGTGAGCGCTATTTGCGTCAGCAACGGGTGAGTAGCGCGGCTGAGAAGTTTGTTGAGGAGCAGAGAACGCTGCATCGGCTCAACCCGATTTGGGATAACCAGTTCATGACttt                               <  3:53755/148‑1 (MQ=255)
                                                                                                              ccGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTAGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGc                             <  1:164827/149‑1 (MQ=255)
                                                                                                                gTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTg                           <  6:45563/149‑1 (MQ=255)
                                                                                                                  cTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGt                                         <  2:76129/133‑1 (MQ=255)
                                                                                                                   tGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCag                                                                                 <  8:209409/92‑1 (MQ=255)
                                                                                                                    gTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAg                        >  7:189726/1‑148 (MQ=255)
                                                                                                                     ttGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGGTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGa                         <  5:187827/146‑1 (MQ=255)
                                                                                                                       ggggAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAg                     <  3:265167/148‑1 (MQ=255)
                                                                                                                        gggAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAggg                   <  8:189726/149‑1 (MQ=255)
                                                                                                                               ggtaataTTTACCCGCGAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATa            <  8:128189/143‑1 (MQ=255)
                                                                                                                               ggAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACTGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCAt             >  5:244664/1‑148 (MQ=255)
                                                                                                                               ggAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttgtt                                                                                          >  2:59429/1‑71 (MQ=255)
                                                                                                                               ggAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttgtt                                                                                          <  1:59429/71‑1 (MQ=255)
                                                                                                                                gAAAGATTTACCCGCCTGTGAGATCGATTTGCGTCAGCAACGGGTTATTAGCGCCGCTGAGTAGTTTGTTGAGGATCAGAGAACGCTGc                                                                       <  4:160521/89‑1 (MQ=255)
                                                                                                                                gAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTAGTTGAGGATCAGAGAACGCTGc                                                                       >  3:160521/1‑89 (MQ=255)
                                                                                                                                gAAAGATTTACCAGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATa            >  7:35388/1‑148 (MQ=255)
                                                                                                                                 aaaGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATa            >  5:306087/1‑147 (MQ=255)
                                                                                                                                      tttACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAAc     >  3:41634/1‑149 (MQ=255)
                                                                                                                                          cccGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTgg  <  6:306087/148‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
GGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGGCTTTCCCATCAGCCGCCGGTGCCCGAACTGGCGAAAGCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAAGATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGG  >  NC_000913/370751‑371036

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