Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,052,886 C→G 5.6% A211A (GCG→GCC)  ubiH ← 2‑octaprenyl‑6‑methoxyphenol hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,052,8860CG 5.6% 304.0 / 7.1 89A211A (GCG→GCC) ubiH2‑octaprenyl‑6‑methoxyphenol hydroxylase
Reads supporting (aligned to +/- strand):  ref base C (41/43);  new base G (3/2);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

CCAGCGTGGGTAATTTTCCCAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATC  >  NC_000913/3052737‑3053035
                                                                                                                                                     |                                                                                                                                                     
ccAGCGTGGGTAATTTTCCCAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATc                                                                                                                                                       >  2:1051626/1‑150 (MQ=255)
        ggTAATTTTCCCAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCgg                                                                                                                                               <  2:1542303/150‑1 (MQ=255)
         gTAATTTTCCCAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCggg                                                                                                                                              >  1:1527207/1‑150 (MQ=255)
              tttCCCAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATg                                                                                                                                         <  2:993939/150‑1 (MQ=255)
                   cAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATGGCCCGCTGGCGATGTTGcc                                                                                                                                    <  1:46361/150‑2 (MQ=255)
                   cAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCg                                                                                                                                    <  1:1035872/150‑1 (MQ=255)
                   cAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCg                                                                                                                                    >  2:102275/1‑150 (MQ=255)
                   cAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCg                                                                                                                                    <  1:698176/150‑1 (MQ=255)
                    aaGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGt                                                                                                                                   <  1:1647164/150‑1 (MQ=255)
                     aGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCGTGTTGCGTa                                                                                                                                  <  2:1755872/150‑1 (MQ=255)
                        cGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAACGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAg                                                                                                                               >  1:1707914/1‑150 (MQ=255)
                         gccagccaAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGc                                                                                                                              >  2:218526/1‑150 (MQ=255)
                         gccagccaAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGc                                                                                                                              >  2:225496/1‑150 (MQ=255)
                             gccaAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTc                                                                                                                          <  1:225496/150‑1 (MQ=255)
                             gccaAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTc                                                                                                                          <  1:218526/150‑1 (MQ=255)
                                aaaGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCCGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCaaa                                                                                                                       >  2:1780340/1‑150 (MQ=255)
                                  aGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAg                                                                                                                     <  2:907379/150‑1 (MQ=255)
                                   ggCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGTCAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAgc                                                                                                                    >  2:602685/1‑150 (MQ=255)
                                    gCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAgcg                                                                                                                   <  1:1672514/150‑1 (MQ=255)
                                           ggAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCtt                                                                                                            >  1:996907/1‑150 (MQ=255)
                                            gAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTc                                                                                                           <  2:1831372/150‑1 (MQ=255)
                                              gTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCAt                                                                                                         >  2:476484/1‑150 (MQ=255)
                                                 cACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATgcg                                                                                                      >  2:187054/1‑150 (MQ=255)
                                                  aCGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATgcgc                                                                                                     <  2:1034350/150‑1 (MQ=255)
                                                   cGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCa                                                                                                    >  1:1099500/1‑150 (MQ=255)
                                                    ggCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCaa                                                                                                   >  1:309746/1‑150 (MQ=255)
                                                      cAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACg                                                                                                 <  2:1262112/150‑1 (MQ=255)
                                                       aaaaCTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACgg                                                                                                <  2:864684/150‑1 (MQ=255)
                                                        aaaCTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGa                                                                                               <  2:1527207/150‑1 (MQ=255)
                                                         aaCTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGat                                                                                              >  2:1662252/1‑149 (MQ=255)
                                                         aaCTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGaa                                                                                              <  1:1391835/150‑1 (MQ=255)
                                                           cTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGt                                                                                            <  2:286444/150‑1 (MQ=255)
                                                              ctcGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCtt                                                                                                            >  2:275916/1‑131 (MQ=255)
                                                              ctcGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCtt                                                                                                            <  1:275916/131‑1 (MQ=255)
                                                                 gTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAAc                                                                                      >  1:1080182/1‑150 (MQ=255)
                                                                 gTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAAc                                                                                      >  1:737865/1‑150 (MQ=255)
                                                                   cACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGt                                                                                    >  1:1717185/1‑150 (MQ=255)
                                                                   cACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGt                                                                                    >  1:1710263/1‑150 (MQ=255)
                                                                   cACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGt                                                                                    >  2:1722785/1‑150 (MQ=255)
                                                                     cTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGTCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTg                                                                                  <  2:936918/150‑1 (MQ=255)
                                                                     cTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTg                                                                                  >  2:287284/1‑150 (MQ=255)
                                                                     cTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCAAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTg                                                                                  >  2:55342/1‑150 (MQ=255)
                                                                       ccACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGc                                                                                <  1:332660/150‑1 (MQ=255)
                                                                       ccACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGc                                                                                >  2:1160790/1‑150 (MQ=255)
                                                                         aCGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCaa                                                                              <  2:1099500/150‑1 (MQ=255)
                                                                            acaacaCCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCa                                                                           >  1:42195/1‑150 (MQ=255)
                                                                            acaacaCCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGAGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCa                                                                           >  1:1839397/1‑150 (MQ=255)
                                                                                cacCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGc                                                                       >  2:83673/1‑150 (MQ=255)
                                                                                cacCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGc                                                                       <  1:1643160/150‑1 (MQ=255)
                                                                                    tctTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAgt                                                                   <  1:1493476/150‑1 (MQ=255)
                                                                                       tCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAgttgt                                                                <  1:1759460/150‑1 (MQ=255)
                                                                                        cgcgCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCa                                                                               <  1:106394/134‑1 (MQ=255)
                                                                                        cgcgCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGAGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCa                                                                               >  2:106394/1‑134 (MQ=255)
                                                                                         gcgcCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTc                                                              <  1:102275/150‑1 (MQ=255)
                                                                                            ccGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTa                                                           <  1:1626119/150‑1 (MQ=255)
                                                                                              gTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAg                                                         <  2:1899197/150‑1 (MQ=255)
                                                                                              gTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAg                                                         <  2:996907/150‑1 (MQ=255)
                                                                                                 ccAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCt                                                      >  2:1161237/1‑150 (MQ=255)
                                                                                                 ccAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCt                                                      >  2:661332/1‑150 (MQ=255)
                                                                                                     tGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCctg                                                  >  1:621454/1‑150 (MQ=255)
                                                                                                      ggATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCctgc                                                 >  1:1490858/1‑150 (MQ=255)
                                                                                                            cacCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAg                                           >  1:288255/1‑150 (MQ=255)
                                                                                                            cacCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAg                                           <  2:1080182/150‑1 (MQ=255)
                                                                                                            cacCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAg                                           <  1:50967/150‑1 (MQ=255)
                                                                                                            cacCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAg                                           >  2:1423800/1‑150 (MQ=255)
                                                                                                            cacCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAg                                           <  2:1490858/150‑1 (MQ=255)
                                                                                                             accagaccagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGt                                          <  2:1379313/150‑1 (MQ=255)
                                                                                                              ccagaccagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTGCTGCTGCCAGTc                                         <  1:476484/150‑1 (MQ=255)
                                                                                                              ccagaccagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTc                                         >  2:1069493/1‑150 (MQ=255)
                                                                                                                 gaccagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAAc                                      >  1:1114648/1‑150 (MQ=255)
                                                                                                                 gaccagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAAc                                      <  2:15755/150‑1 (MQ=255)
                                                                                                                 gaccagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAAc                                      >  2:552630/1‑150 (MQ=255)
                                                                                                                    cagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAAcgcc                                   <  1:1722785/150‑1 (MQ=255)
                                                                                                                        gAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCAc                               <  1:407789/150‑1 (MQ=255)
                                                                                                                         aaCAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCAcg                              >  1:1005775/1‑150 (MQ=255)
                                                                                                                           cAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCAcgcg                            >  1:703300/1‑150 (MQ=255)
                                                                                                                              cgtccgtcAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTg                         <  1:287284/150‑1 (MQ=255)
                                                                                                                               gtccgtcAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTgg                        <  1:921911/150‑1 (MQ=255)
                                                                                                                               gtccgtcAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTgg                        <  1:661332/150‑1 (MQ=255)
                                                                                                                                 ccgtcAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCt                      >  2:429119/1‑150 (MQ=255)
                                                                                                                                    tcAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAAc                   <  2:1710263/150‑1 (MQ=255)
                                                                                                                                    tcAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAAc                   <  2:1717185/150‑1 (MQ=255)
                                                                                                                                      aGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGc                 >  1:659620/1‑150 (MQ=255)
                                                                                                                                       gctttacgCAACATGGCCCGCTGGCGATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCt                <  1:1211155/143‑1 (MQ=255)
                                                                                                                                       gctttacgCAACATGGCCCGCTGGCGATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCt                >  2:1211155/8‑150 (MQ=255)
                                                                                                                                         tttacgCAACATGGCCCGCTGGCGATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGa              >  1:1525799/6‑150 (MQ=255)
                                                                                                                                          aTCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGaa             >  1:1984721/1‑150 (MQ=255)
                                                                                                                                             ggCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATgg          >  2:1077041/1‑150 (MQ=255)
                                                                                                                                              gCAACATGGCCCGCTGGCGATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCCGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATggg         >  2:638180/1‑150 (MQ=255)
                                                                                                                                                aaCATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTg       <  1:1077041/150‑1 (MQ=255)
                                                                                                                                                     cGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATc  <  1:1051626/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAGCGTGGGTAATTTTCCCAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATC  >  NC_000913/3052737‑3053035

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