Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,052,893 G→T 5.5% P209Q (CCG→CAG)  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,8930GT 5.5% 307.3 / 5.6 90P209Q (CCG→CAG) ubiH2‑octaprenyl‑6‑methoxyphenol hydroxylase
Reads supporting (aligned to +/- strand):  ref base G (41/44);  new base T (3/2);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

GGTAATTTTCCCAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATCAGCTGCT  >  NC_000913/3052745‑3053042
                                                                                                                                                    |                                                                                                                                                     
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:698176/150‑1 (MQ=255)
           cAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCg                                                                                                                                           <  1:1035872/150‑1 (MQ=255)
           cAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCg                                                                                                                                           >  2:102275/1‑150 (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                                                                              <  1:1643160/150‑1 (MQ=255)
                                                                        cacCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGc                                                                              >  2:83673/1‑150 (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:996907/150‑1 (MQ=255)
                                                                                      gTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAg                                                                <  2:1899197/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:50967/150‑1 (MQ=255)
                                                                                                    cacCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAg                                                  >  2:1423800/1‑150 (MQ=255)
                                                                                                    cacCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAg                                                  <  2:1490858/150‑1 (MQ=255)
                                                                                                    cacCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAg                                                  <  2:1080182/150‑1 (MQ=255)
                                                                                                    cacCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAg                                                  >  1:288255/1‑150 (MQ=255)
                                                                                                     accagaccagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGt                                                 <  2:1379313/150‑1 (MQ=255)
                                                                                                      ccagaccagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTGCTGCTGCCAGTc                                                <  1:476484/150‑1 (MQ=255)
                                                                                                      ccagaccagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTc                                                >  2:1069493/1‑150 (MQ=255)
                                                                                                         gaccagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAAc                                             <  2:15755/150‑1 (MQ=255)
                                                                                                         gaccagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAAc                                             >  2:552630/1‑150 (MQ=255)
                                                                                                         gaccagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAAc                                             >  1:1114648/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:661332/150‑1 (MQ=255)
                                                                                                                       gtccgtcAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTgg                               <  1:921911/150‑1 (MQ=255)
                                                                                                                         ccgtcAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCt                             >  2:429119/1‑150 (MQ=255)
                                                                                                                            tcAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAAc                          <  2:1717185/150‑1 (MQ=255)
                                                                                                                            tcAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAAc                          <  2:1710263/150‑1 (MQ=255)
                                                                                                                              aGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGc                        >  1:659620/1‑150 (MQ=255)
                                                                                                                               gctttacgCAACATGGCCCGCTGGCGATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCt                       >  2:1211155/8‑150 (MQ=255)
                                                                                                                               gctttacgCAACATGGCCCGCTGGCGATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCt                       <  1:1211155/143‑1 (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)
                                                                                                                                               ccAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATCAg       <  2:1346026/150‑1 (MQ=255)
                                                                                                                                                 aGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATCAgct     >  2:258982/1‑150 (MQ=255)
                                                                                                                                                    gggCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATCAgctgct  >  1:139182/1‑150 (MQ=255)
                                                                                                                                                    gggCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATCAgctgct  <  1:83673/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGTAATTTTCCCAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATCAGCTGCT  >  NC_000913/3052745‑3053042

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: ‑