Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,447,767 T→C R83R (CGT→CGC feaB → phenylacetaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,447,7670TC98.5% 238.4 / ‑4.1 67R83R (CGT→CGCfeaBphenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (27/39);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-01

AAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCC  >  NC_000913/1447620‑1447914
                                                                                                                                                   |                                                                                                                                                   
aaaaCGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGa                                                                                                                                                   >  1:838344/1‑150 (MQ=255)
 aaaCGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAAc                                                                                                                                                 <  2:721763/151‑1 (MQ=255)
 aaaCGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAAc                                                                                                                                                 <  2:1956/151‑1 (MQ=255)
     ggTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTAt                                                                                                                                             <  2:245451/151‑1 (MQ=255)
      gTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTAtt                                                                                                                                            >  1:1035613/1‑151 (MQ=255)
        tGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTc                                                                                                                                           <  2:756518/150‑1 (MQ=255)
                 ttGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGttt                                                                                                                                 >  2:408436/1‑151 (MQ=255)
                  tGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGtt                                                                                                                                  >  2:916180/1‑149 (MQ=255)
                  tGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGtt                                                                                                                                  <  1:916165/149‑1 (MQ=255)
                     tCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCt                                                                                                                             >  2:187150/1‑151 (MQ=255)
                      ccGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCt                                                                                                                             <  2:775055/150‑1 (MQ=255)
                      ccGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTg                                                                                                                            >  1:281991/1‑151 (MQ=255)
                        ggCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGAt                                                                                                                          <  1:1025337/151‑1 (MQ=255)
                         gCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATc                                                                                                                         <  1:1060509/151‑1 (MQ=255)
                             ccGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCtggt                                                                                                                     >  2:333654/1‑151 (MQ=255)
                              cGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCtggtg                                                                                                                    >  2:1117764/1‑151 (MQ=255)
                               gggCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCtggtgg                                                                                                                   <  2:526925/151‑1 (MQ=255)
                                 gCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGag                                                                                                                 <  2:223959/151‑1 (MQ=255)
                                  cAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGCTCTGGTGGagc                                                                                                                >  1:530706/1‑151 (MQ=255)
                                    aGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGagcag                                                                                                              <  2:106858/151‑1 (MQ=255)
                                    aGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGagc                                                                                                                >  1:844323/1‑149 (MQ=255)
                                      aaaTTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGagcagc                                                                                                             >  1:751270/1‑150 (MQ=255)
                                           gCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGt                                                                                                        <  2:240036/150‑1 (MQ=255)
                                             gTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTgag                                                                                                     <  1:397069/151‑1 (MQ=255)
                                                       gATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTCCCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGgcgc                                                                                           <  2:243613/151‑1 (MQ=255)
                                                       gATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGgcg                                                                                            <  1:1003455/150‑1 (MQ=255)
                                                            cAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGagc                                                                                                                <  2:670037/125‑1 (MQ=255)
                                                            cAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGagc                                                                                                                >  1:670027/1‑125 (MQ=255)
                                                            cAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTg                                                                                      <  1:1196125/151‑1 (MQ=255)
                                                               cGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGaa                                                                                   <  2:281994/151‑1 (MQ=255)
                                                                gAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGaa                                                                                   >  1:508065/1‑150 (MQ=255)
                                                                   gCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAAccc                                                                               >  1:594705/1‑151 (MQ=255)
                                                                   gCGGATGTAGATAACGCAGTCATGTCTGCCTGGAGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAAccc                                                                               >  1:230172/1‑151 (MQ=255)
                                                                    cGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCt                                                                              <  1:1133339/151‑1 (MQ=255)
                                                                     ggATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTg                                                                             <  1:627948/151‑1 (MQ=255)
                                                                      gaggTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACGGGAAACCCttg                                                                            <  1:1005198/148‑3 (MQ=255)
                                                                      gATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTgg                                                                            >  2:1114558/1‑151 (MQ=255)
                                                                       aTGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTgg                                                                            >  2:267822/1‑150 (MQ=255)
                                                                               aaCGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAgcaaggc                                                                    >  2:1180743/1‑150 (MQ=255)
                                                                                 cGCAGTCATGTCTGCCTGGAGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAgcaaggcaag                                                                 <  1:555005/151‑1 (MQ=255)
                                                                                     gTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCaa                                                             <  2:595401/151‑1 (MQ=255)
                                                                                      tCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAAt                                                            >  1:334542/1‑151 (MQ=255)
                                                                                       cATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAAtt                                                           <  1:200499/151‑1 (MQ=255)
                                                                                            cTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCAt                                                      >  2:186503/1‑151 (MQ=255)
                                                                                            cTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCAt                                                      >  2:185916/1‑151 (MQ=255)
                                                                                            cTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCAt                                                      >  2:529878/1‑151 (MQ=255)
                                                                                              gCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCa                                                       <  1:464546/148‑1 (MQ=255)
                                                                                              gCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCa                                                       <  1:1123630/148‑1 (MQ=255)
                                                                                              gCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCa                                                       <  1:698683/148‑1 (MQ=255)
                                                                                                  ggCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTcc                                                  <  2:329880/149‑1 (MQ=255)
                                                                                                      ggCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCt                                            <  2:410038/151‑1 (MQ=255)
                                                                                                             gTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAg                                     <  2:356207/151‑1 (MQ=255)
                                                                                                                 cgcgTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTg                                   <  2:806721/149‑1 (MQ=255)
                                                                                                                   cgTCGCTGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGCt                               >  2:848199/1‑151 (MQ=255)
                                                                                                                       gCTGAGCCGGGCGATTACCCGCAGAGCGCGAACGAATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGCTGTAc                           >  2:545887/1‑151 (MQ=255)
                                                                                                                         tGGGCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGCTGTACGc                         >  2:422006/1‑151 (MQ=255)
                                                                                                                          gggCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGCTGTACGCt                        <  2:319751/151‑1 (MQ=255)
                                                                                                                           ggCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGCTGTACGCTg                       <  1:369292/151‑1 (MQ=255)
                                                                                                                           ggCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGCTGTACGCTg                       <  2:58510/151‑1 (MQ=255)
                                                                                                                           ggCCGGGCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGCTGTACGCTg                       >  2:771663/1‑151 (MQ=255)
                                                                                                                                 gCGATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTgg                 <  1:728324/151‑1 (MQ=255)
                                                                                                                                   gATTACCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGAt               <  2:684991/151‑1 (MQ=255)
                                                                                                                                       aCCCGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGt           <  2:400959/151‑1 (MQ=255)
                                                                                                                                        cccGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGtt          >  1:172833/1‑151 (MQ=255)
                                                                                                                                         ccGCAGAGCGCGAACGTATTCTGCTACGTTATGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTta         <  1:1089411/151‑1 (MQ=255)
                                                                                                                                          cGCAGAGCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTtat        <  1:344468/151‑1 (MQ=255)
                                                                                                                                             agagCGCGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATAcc     <  2:1025017/151‑1 (MQ=255)
                                                                                                                                                gcgcGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATAccgcc  <  2:1080528/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGGAGCAGCACAGTGAGGAGCTGGCGCAACTGGAAACCCTGGAGCAAGGCAAGTCAATTGCCATTTCCCGTGCTTTTGAAGTGGGCTGTACGCTGAACTGGATGCGTTATACCGCC  >  NC_000913/1447620‑1447914

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: