Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 508,754 T→G 5.4% intergenic (+195/+121) ybaQ → / ← copA putative DNA‑binding transcriptional regulator/copper transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913508,7540TG 5.4% 255.4 / 4.0 74intergenic (+195/+121)ybaQ/copAputative DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base T (39/31);  new base G (2/2);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAAC  >  NC_000913/508607‑508902
                                                                                                                                                   |                                                                                                                                                    
gcatccggcatccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCtt                                                                                                                                                     <  2:537315/149‑1 (MQ=255)
gcatccggcatccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCtt                                                                                                                                                     <  4:347792/149‑1 (MQ=255)
    ccggcatccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCa                                                                                                                                                 <  2:387980/149‑1 (MQ=255)
      ggcatccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAgg                                                                                                                                               >  3:295248/1‑149 (MQ=255)
                    aTTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCtgtg                                                                                                                                 >  1:460356/1‑149 (MQ=255)
                     ttGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGc                                                                                                                                <  2:534960/149‑1 (MQ=255)
                     ttGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGc                                                                                                                                <  4:295248/149‑1 (MQ=255)
                      tGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGcc                                                                                                                               <  2:85487/149‑1 (MQ=255)
                          tGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGaa                                                                                                                           >  2:296746/1‑149 (MQ=255)
                          tGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGaa                                                                                                                           >  2:527957/1‑149 (MQ=255)
                           gATGCGACGCTTACCCGTCTTACCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAAc                                                                                                                          <  3:324109/149‑1 (MQ=255)
                                   gCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGcc                                                                                                                  >  3:423629/1‑149 (MQ=255)
                                       aCCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGAt                                                                                                              <  2:460356/149‑1 (MQ=255)
                                        cccGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATa                                                                                                             >  4:342482/1‑149 (MQ=255)
                                         ccGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATaa                                                                                                            >  4:178438/1‑149 (MQ=255)
                                         ccGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATaa                                                                                                            <  1:711251/149‑1 (MQ=255)
                                          cGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAg                                                                                                           >  3:433814/1‑149 (MQ=255)
                                           gTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAgg                                                                                                          <  2:387399/149‑1 (MQ=255)
                                           gTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAgg                                                                                                          >  1:33270/1‑149 (MQ=255)
                                            tCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGc                                                                                                         <  4:423629/149‑1 (MQ=255)
                                              ttATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGt                                                                                                       >  1:497427/1‑149 (MQ=255)
                                               tATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGtt                                                                                                      <  3:613168/149‑1 (MQ=255)
                                                aTCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGttt                                                                                                     >  4:270152/1‑149 (MQ=255)
                                                   aGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATc                                                                                                                                                  <  4:122880/101‑1 (MQ=255)
                                                   aGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATc                                                                                                                                                  >  3:122880/1‑101 (MQ=255)
                                                     gCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                <  4:369834/149‑1 (MQ=255)
                                                     gCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                >  2:561900/1‑149 (MQ=255)
                                                     gCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                >  1:455783/1‑149 (MQ=255)
                                                        tACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCa                                                                                             >  4:251312/1‑149 (MQ=255)
                                                        tACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCa                                                                                             >  4:29785/1‑149 (MQ=255)
                                                        tACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCa                                                                                             >  1:154794/1‑149 (MQ=255)
                                                        tACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCa                                                                                             >  2:502961/1‑149 (MQ=255)
                                                         aCAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCAt                                                                                            >  2:385265/1‑149 (MQ=255)
                                                              cTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                       >  3:634663/1‑149 (MQ=255)
                                                               ttGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                      >  3:569780/1‑149 (MQ=255)
                                                               ttGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                      >  1:695993/1‑149 (MQ=255)
                                                                   gCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                                         <  1:256902/142‑1 (MQ=255)
                                                                   gCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                                         >  2:256902/1‑142 (MQ=255)
                                                                     cTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATg                                                                                >  3:460517/1‑149 (MQ=255)
                                                                      tGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAgg                                                                                                          <  4:691842/122‑1 (MQ=255)
                                                                      tGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAgg                                                                                                          >  3:691842/1‑122 (MQ=255)
                                                                       gAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCg                                                                              >  1:764090/1‑149 (MQ=255)
                                                                         aCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTc                                                                            <  1:11145/149‑1 (MQ=255)
                                                                           cGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGa                                                                          <  3:342482/149‑1 (MQ=255)
                                                                            gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGAt                                                                         <  2:473565/149‑1 (MQ=255)
                                                                                         ggCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCa                                                            >  4:408054/1‑149 (MQ=255)
                                                                                         ggCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCa                                                            >  1:692364/1‑149 (MQ=255)
                                                                                             ttttCGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCt                                                        >  2:329013/1‑149 (MQ=255)
                                                                                                     gcatccggcatccgAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAagcagc                                                <  2:318329/149‑1 (MQ=255)
                                                                                                     gcatccggcatccgAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAagcagc                                                <  2:763409/149‑1 (MQ=255)
                                                                                                           ggcatccgAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATcc                                          >  3:17918/1‑149 (MQ=255)
                                                                                                           ggcatccgAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATcc                                          >  3:643020/1‑149 (MQ=255)
                                                                                                            gcatccgAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCg                                         <  1:561900/149‑1 (MQ=255)
                                                                                                                ccgAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAat                                     <  1:527957/149‑1 (MQ=255)
                                                                                                                   aGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAatgat                                  <  1:449600/149‑1 (MQ=255)
                                                                                                                    gCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAatgatg                                 <  3:270152/149‑1 (MQ=255)
                                                                                                                        cATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACt                             <  3:251312/149‑1 (MQ=255)
                                                                                                                         aTTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCAcc                                                          >  2:332803/1‑119 (MQ=255)
                                                                                                                         aTTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCAcc                                                          <  1:332803/119‑1 (MQ=255)
                                                                                                                         aTTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACtt                            >  4:419453/1‑149 (MQ=255)
                                                                                                                           ggCCTGATAAGACGCGTAATCGTCGCATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAagcagc                                                >  4:533553/2‑127 (MQ=21)
                                                                                                                           ggCCTGATAAGACGCGTAAGCGTCGCATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAagcagc                                                <  3:533553/126‑1 (MQ=21)
                                                                                                                              ctgataagACGCGTAAGCGTCGCATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTcc                       >  3:438536/8‑149 (MQ=37)
                                                                                                                                  tGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCg                   <  1:276351/149‑1 (MQ=255)
                                                                                                                                   gCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCgg                  <  1:565132/149‑1 (MQ=255)
                                                                                                                                     gACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGtt                >  2:494/1‑149 (MQ=255)
                                                                                                                                     gACGCGTAAGCGTCGCATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGt                 <  4:531713/148‑1 (MQ=37)
                                                                                                                                        gCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCAACCGCAATGATGTACTTATTCCTTCGGTTTaa             <  1:385265/149‑1 (MQ=255)
                                                                                                                                          ttACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAAc           <  1:757529/149‑1 (MQ=255)
                                                                                                                                          ttACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAAc           <  3:408054/149‑1 (MQ=255)
                                                                                                                                           tACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAAcc          <  2:695993/149‑1 (MQ=255)
                                                                                                                                            aCGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCg         >  3:339019/1‑149 (MQ=255)
                                                                                                                                              gcgTCTTATCAGGCCTACAAACCTGTGCCTGACCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCgca       >  4:701147/1‑149 (MQ=255)
                                                                                                                                                 tCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCgcagca    >  2:79439/1‑149 (MQ=255)
                                                                                                                                                 tCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCgcagca    >  2:310280/1‑149 (MQ=255)
                                                                                                                                                 tCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCgcagca    >  2:106573/1‑149 (MQ=255)
                                                                                                                                                   ttATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAAc  >  1:195367/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAAC  >  NC_000913/508607‑508902

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