Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 138,797 A→G 6.4% intergenic (+164/+38) cueO → / ← gcd multicopper oxidase (laccase)/glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913138,7970AG 6.4% 316.7 / 8.3 93intergenic (+164/+38)cueO/gcdmulticopper oxidase (laccase)/glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (46/41);  new base G (4/2);  total (50/43)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

CCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTA  >  NC_000913/138667‑138934
                                                                                                                                  |                                                                                                                                         
ccGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGTGCGCGtt                                                                                                                          >  1:185969/1‑148 (MQ=255)
    gCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGTGCGCGTTGttc                                                                                                                      <  2:185969/148‑3 (MQ=255)
     cATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaacaa                                                                                                                                    >  2:396684/1‑133 (MQ=255)
     cATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaacaa                                                                                                                                    <  1:396684/133‑1 (MQ=255)
              gatgaCAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAg                                                                                                            >  4:493149/1‑148 (MQ=255)
                   cAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGt                                                                                                       >  3:813553/1‑148 (MQ=255)
                       cGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTa                                                                                                   <  3:493149/148‑1 (MQ=255)
                            aTCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTca                                                                                              <  4:813553/148‑1 (MQ=255)
                                  cTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGt                                                                                                       >  1:982320/1‑133 (MQ=255)
                                  cTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGt                                                                                                       <  2:982320/133‑1 (MQ=255)
                                   tGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCGTACTTCTCatcatc                                                                                       >  1:802451/1‑148 (MQ=255)
                                     aTGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATg                                                                                                                                         <  4:761512/96‑1 (MQ=255)
                                     aTGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATg                                                                                                                                         >  3:761512/1‑96 (MQ=255)
                                     aTGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaacaa                                                                                                                                    >  3:395151/1‑101 (MQ=255)
                                     aTGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaacaa                                                                                                                                    <  4:395151/101‑1 (MQ=255)
                                          cTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTc                                                                                               <  3:1041865/133‑1 (MQ=255)
                                          cTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTc                                                                                               >  4:1041865/1‑133 (MQ=255)
                                           tACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAcg                                                                                                                                  >  3:1077968/1‑97 (MQ=255)
                                           tACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAcg                                                                                                                                  <  4:1077968/97‑1 (MQ=255)
                                            aCGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCaa                                                                                                         <  4:1068818/121‑1 (MQ=255)
                                            aCGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCaa                                                                                                         >  3:1068818/1‑121 (MQ=255)
                                               cTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTg                                                                                                              >  3:821869/1‑113 (MQ=255)
                                               cTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTg                                                                                                              <  4:821869/113‑1 (MQ=255)
                                                 tATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatca                                                                                         <  4:210064/132‑1 (MQ=255)
                                                 tATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatca                                                                                         >  3:210064/1‑132 (MQ=255)
                                                      ggCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCcaca                                                                    <  4:319778/148‑1 (MQ=255)
                                                      ggCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCcaca                                                                    >  3:319778/1‑148 (MQ=255)
                                                           aCGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGt                                                                                                                           >  3:1049601/1‑88 (MQ=255)
                                                           aCGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGt                                                                                                                           <  4:1049601/88‑1 (MQ=255)
                                                            cGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCa                                                                                                                              <  4:156820/84‑1 (MQ=255)
                                                            cGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCa                                                                                                                              >  3:156820/1‑84 (MQ=255)
                                                              cTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTcac                                                                                             <  3:242870/115‑1 (MQ=255)
                                                              cTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTcac                                                                                             >  4:242870/1‑115 (MQ=255)
                                                               tAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCAGCCGGCAGCGCATAATCCACAATATAGTCg                                                           <  2:802451/148‑1 (MQ=255)
                                                                aGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTc                                                                                                      >  2:463357/1‑104 (MQ=255)
                                                                aGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTc                                                                                                      <  1:463357/104‑1 (MQ=255)
                                                                    gtgCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCcaca                                                                    <  2:554256/134‑1 (MQ=255)
                                                                    gtgCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCcaca                                                                    >  1:554256/1‑134 (MQ=255)
                                                                         aTCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCg                                                 >  4:526541/1‑148 (MQ=255)
                                                                          tCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGt                                                >  2:108425/1‑148 (MQ=255)
                                                                           ccATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgc                                                                                >  3:903796/1‑115 (MQ=255)
                                                                           ccATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgc                                                                                <  4:903796/115‑1 (MQ=255)
                                                                           ccATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATa                                                                           >  2:54766/1‑120 (MQ=255)
                                                                           ccATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATa                                                                           <  1:54766/120‑1 (MQ=255)
                                                                           ccATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGt                                                             >  2:326867/1‑134 (MQ=255)
                                                                           ccATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGt                                                             <  1:326867/134‑1 (MQ=255)
                                                                            cATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatcat                                                                                        <  2:883753/106‑1 (MQ=255)
                                                                            cATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatcat                                                                                        >  1:883753/1‑106 (MQ=255)
                                                                              ttGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTca                                                                                              <  3:1132739/98‑1 (MQ=255)
                                                                              ttGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTca                                                                                              >  4:1132739/1‑98 (MQ=255)
                                                                              ttGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCAt                                                      >  2:1111630/1‑138 (MQ=255)
                                                                              ttGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCAt                                                      <  1:1111630/138‑1 (MQ=255)
                                                                               tGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGt                                                                                                       >  3:154027/1‑88 (MQ=255)
                                                                               tGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGt                                                                                                       <  4:154027/88‑1 (MQ=255)
                                                                               tGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCaa                                           <  3:526541/148‑1 (MQ=255)
                                                                                  ttttGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAatat                                                                <  2:270172/124‑1 (MQ=255)
                                                                                  ttttGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAatat                                                                >  1:270172/1‑124 (MQ=255)
                                                                                   tttGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTg                                                                                                              <  1:1109900/77‑1 (MQ=255)
                                                                                   tttGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTg                                                                                                              >  2:1109900/1‑77 (MQ=255)
                                                                                   tttGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCGCATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGa                                       >  2:867861/1‑148 (MQ=255)
                                                                                        acaATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCtt                                                                                                    >  1:296812/1‑82 (MQ=255)
                                                                                        acaATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCtt                                                                                                    <  2:296812/82‑1 (MQ=255)
                                                                                                    gCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatc                                                                                          <  1:805697/80‑1 (MQ=255)
                                                                                                    gCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAatat                                                                <  3:524212/106‑1 (MQ=255)
                                                                                                    gCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAatat                                                                >  4:524212/1‑106 (MQ=255)
                                                                                                    gCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCgg                      >  3:618921/1‑148 (MQ=255)
                                                                                                    gCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCgg                      <  4:618921/148‑1 (MQ=255)
                                                                                                    gCTGGAGAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatc                                                                                          >  2:805697/1‑80 (MQ=255)
                                                                                                          tAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTACGGAGATca                >  2:255189/1‑148 (MQ=255)
                                                                                                              gCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTcac                                                                                             >  1:417489/1‑67 (MQ=255)
                                                                                                              gCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTcac                                                                                             <  2:417489/67‑1 (MQ=255)
                                                                                                               cGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACcaca           <  1:108425/148‑1 (MQ=255)
                                                                                                                 ttCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGt                                  >  4:430575/1‑123 (MQ=255)
                                                                                                                 ttCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGt                                  <  3:430575/123‑1 (MQ=255)
                                                                                                                 ttCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATca                <  3:6247/141‑1 (MQ=255)
                                                                                                                 ttCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATca                >  4:6247/1‑141 (MQ=255)
                                                                                                                 ttCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATa         >  4:761332/1‑148 (MQ=255)
                                                                                                                  tCATGCCGGATGCGGCGTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAg                                                              <  3:1097504/94‑1 (MQ=37)
                                                                                                                  tCATGCCGGATGCGGCGTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAg                                                              >  4:1097504/1‑94 (MQ=37)
                                                                                                                   catgccggATGCGGCGTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGaa                                      >  2:686477/4‑117 (MQ=39)
                                                                                                                   catgccggATGCGGCGTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTg                         <  3:210769/127‑1 (MQ=39)
                                                                                                                   catgccggATGCGGCGTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTg                         >  4:210769/4‑130 (MQ=255)
                                                                                                                    atgccggATGCGGCGTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATAc        >  1:637765/3‑146 (MQ=255)
                                                                                                                     cgccgcATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGc     >  1:560743/1‑148 (MQ=255)
                                                                                                                     cgccgcATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGc     >  2:1042054/1‑148 (MQ=255)
                                                                                                                     cgccgcATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGc     >  2:630606/1‑148 (MQ=255)
                                                                                                                        cgcATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTa  <  1:867861/148‑1 (MQ=255)
                                                                                                                          cATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGgag                    >  1:1111764/1‑128 (MQ=255)
                                                                                                                          cATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGgag                    <  2:1111764/128‑1 (MQ=255)
                                                                                                                             ccGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCgg                      >  3:824803/1‑123 (MQ=255)
                                                                                                                             ccGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCgg                      <  4:824803/123‑1 (MQ=255)
                                                                                                                                gCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCtt                                                   >  3:684914/1‑91 (MQ=255)
                                                                                                                                gCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCtt                                                   <  4:684914/91‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                         
CCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTA  >  NC_000913/138667‑138934

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