Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,239,609 T→G 8.5% intergenic (+64/‑335) alx → / → sstT putative membrane‑bound redox modulator/sodium:serine/threonine symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,239,6090TG 8.5% 210.8 / 12.9 71intergenic (+64/‑335)alx/sstTputative membrane‑bound redox modulator/sodium:serine/threonine symporter
Reads supporting (aligned to +/- strand):  ref base T (32/33);  new base G (3/3);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-01

TGGCGGGCGTAGCAGAGCGTTTCTCGATGCTCAAATATGGCCTGGCGGTGATTCTGGTGTTTATCGGTATCAAGATGCTGATTGTCGACTTCTACCATATTCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTAAAATCCAGCATTTTTCGCTTCCCGAAGCTGTAACTTTCCTTATACTCGAC  >  NC_000913/3239343‑3239859
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          
tgGCGGGCGTAGCAGAGCGTTTCTCGATGCTCAAATATGGCCTGGCGGTGATTCTGGTGTTTATCGGTATCAAGATGCTGATTGTCGACTTCTACCATATTCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGTGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaa                                                                                                                                                                                                                                                    >  1:98770/1‑274 (MQ=255)
gggggggcgTTGCAGAGCGGTTCTCGATGCTCAAATATGGCCTGGCGGTGATTCTGGTGTTTATCGGTATCAAGATGCTGATTGTCGACTTCTACCATAGTCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGga                                                                                                                                                                                                                                                     <  4:476743/270‑1 (MQ=255)
 tgcgggcgTAGCCGAGCGTTTCTCGATGCTCCAATATGGCCTTGCGTGTATTCTGTGGTTTATCGGTATCAAGATGTCTATTGTCGACTTCTACCATATTCCACTCGCCGTCTCGCTGGGCTTTGTGTTTGGCATTCTTGTGATGACTTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                                                                                                                                                                                                                   <  4:509060/274‑1 (MQ=255)
       cgTAGCAGAGCGTTTCTCGATGCTCAAATATGGCCTGGCGGTGATTCTGGTGTTTATCGGTATCAAGATGCTGATTGTCGACTTCTACCATATTCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCg                                                                                                                                                                                                                                             >  2:22136/1‑275 (MQ=255)
               agCGTTTCTCGATGCTCAAATATGGCCTGGCGGTGATTCTGGTGTTTATCGGTATCAAGATGCTGATTGTCGACTTCTACCATATTCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGGGAg                                                                                                                                                                                                                                      >  2:524444/1‑274 (MQ=255)
                      ctcGATGCTCAAATATGGCCTGGCGGTGATTCTGGTGTTTATCGGTATCAAGATGCTGATTGTCGACTTCTACCATATTCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTTTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGCGGATCGGGGTGAGGGTaa                                                                                                                                                                                                                                 >  4:305973/1‑272 (MQ=255)
                                                 gtttttGGTTTTTTTCTGTTTCAAGTTGCTTATTGTGTCTTTCTCCCTGATTCCAATCGCCGTCTCGTTGGGCGTGGTGTTTGGAATTCTGGGGATGAGGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTTGTTAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCt                                                                                                                                                                                                   <  4:524443/270‑1 (MQ=255)
                                                 gcttttGGTGTTTATCGGTACCAAGATGCTGATTGTCGACTTCTACCATATTCCAATCGCCGTCTCGCTTGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGTCTGGGTTAATTATCGGCATGTTAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGGGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCt                                                                                                                                                                                                   <  4:281940/270‑1 (MQ=255)
                                                             tATCGGTATCAAGATGCGGTTTGTCGACTTCTACCATATTCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAGCGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAccc                                                                                                                                                                                        <  2:375976/274‑1 (MQ=255)
                                                                                           cTACCATATTCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTc                                                                                                                                                                                                       >  4:420502/1‑229 (MQ=255)
                                                                                               cATATTCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCt                                                                                                                                                                                                   <  2:420503/229‑1 (MQ=255)
                                                                                                        aTCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGGGGGTAATTTTTCGCACCGATGCGGGCCTGGTCCCCTCACCCTAACCCTCCCCCCAAAAGGGGCGAGGGGACTGCCCCATTCCtttttgtaa                                                                                                                                            >  4:49151/1‑271 (MQ=255)
                                                                                                         ccGCCTTGTCGGTGGTGGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTTGATAATTTTTAATCTGCCTAAGCCGGGGACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGGGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtg                                                                                                                                           <  4:151750/274‑1 (MQ=39)
                                                                                                                ttcgcgggtgtgggtgTTTGTACTTCTGGGAATGGCGTTTATTCTCTACGCGTGTGTGAGTTATCGGCATGATAAGCAGGGGCGTGGGTAATTTTTGATCGGCCTATTCCGTGTACCCTGTCATTAACGTGAGAAACAGTTTCTCCCTCTCCTTTCCATGGAGAGGGTCGAGGTGTGGGTAATTTTTCGCACCGATGATGGCCTGATCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtcat                                                                                                                                    <  4:356619/263‑1 (MQ=255)
                                                                                                                                   ggCATTCTGTTGATGACGTTTATTATCAACGCCTGGGGGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACg                                                                                                                  <  2:330346/274‑1 (MQ=255)
                                                                                                                                                      ttattaTCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                                                                                                                                                                                                                      <  5:529658/123‑1 (MQ=255)
                                                                                                                                                      ttattaTCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAggg                                                                                                                                                                                                                                                      >  6:529658/1‑123 (MQ=255)
                                                                                                                                                       tattatCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCaa                                                                                             <  2:250326/275‑1 (MQ=255)
                                                                                                                                                       tattatCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGgaga                                                                                                                                                                                                                                                   >  7:431859/1‑125 (MQ=255)
                                                                                                                                                           atCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGgag                                                                                                                                                                                                                                                    >  8:200763/1‑120 (MQ=255)
                                                                                                                                                           atCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTGGAAGGgag                                                                                                                                                                                                                                                    <  7:200763/120‑1 (MQ=255)
                                                                                                                                                               aCGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                                                                                                                                                                                                                            >  5:48770/1‑124 (MQ=255)
                                                                                                                                                               aCGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                                                                                                                                                                                                                            <  6:48770/124‑1 (MQ=255)
                                                                                                                                                                            ttATCGGCATGAGAAGCAGCGGGGTGGATAATTTTTAGTCTGCCTAAGCCTTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGttt                                                                        <  2:49152/275‑1 (MQ=255)
                                                                                                                                                                             tATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTc                                                                                                                                                                                                                                                           <  7:139215/95‑1 (MQ=255)
                                                                                                                                                                             tATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTc                                                                                                                                                                                                                                                           >  8:139215/1‑95 (MQ=255)
                                                                                                                                                                             tATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                                                                                                                                                                                                              >  7:25937/1‑108 (MQ=255)
                                                                                                                                                                             tATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                                                                                                                                                                                                              <  8:25937/108‑1 (MQ=255)
                                                                                                                                                                               tCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTcc                                                                                                                                                                                                                                                          >  6:979683/1‑94 (MQ=255)
                                                                                                                                                                               tCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTcc                                                                                                                                                                                                                                                          <  5:979683/94‑1 (MQ=255)
                                                                                                                                                                                  gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGagag                                                                                                                                                                                                                                                  <  7:57485/99‑1 (MQ=255)
                                                                                                                                                                                  gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGagag                                                                                                                                                                                                                                                  >  8:57485/1‑99 (MQ=255)
                                                                                                                                                                                                    tGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCCCCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGTCGACCGTGGGCTTTATTTTTAGTTCAAAGAAATGTTCCCATTAAGTTTCCAAATTTGACCTTTCTTCaaatt                                                >  4:168552/1‑273 (MQ=38)
                                                                                                                                                                                                           tttttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATTTTGTCATCAGTCTGGAAGCCGCCCGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAATTATTACCTCTCTTTAAAATCCAACAt                                         >  4:492715/1‑275 (MQ=255)
                                                                                                                                                                                                                   cTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                                                                                                                                                                                                                   <  6:990253/65‑1 (MQ=255)
                                                                                                                                                                                                                   cTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaga                                                                                                                                                                                                                                                   >  5:990253/1‑65 (MQ=255)
                                                                                                                                                                                                                   cTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                                                                                                                                                                                                              <  8:322657/70‑1 (MQ=255)
                                                                                                                                                                                                                   cTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTc                                                                                                                                                                                                                                              >  7:322657/1‑70 (MQ=255)
                                                                                                                                                                                                                         aaGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTc                                                                                                                                                                                                                           >  5:950831/1‑83 (MQ=255)
                                                                                                                                                                                                                         aaGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTc                                                                                                                                                                                                                           <  6:950831/83‑1 (MQ=255)
                                                                                                                                                                                                                           gCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                                                                                                                                                         <  7:252187/51‑1 (MQ=255)
                                                                                                                                                                                                                           gCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                                                                                                                                                                                         >  8:252187/1‑51 (MQ=255)
                                                                                                                                                                                                                               tgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACg                                                                                                                                                                            >  6:453659/1‑124 (MQ=255)
                                                                                                                                                                                                                               tgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACg                                                                                                                                                                            <  5:453659/124‑1 (MQ=255)
                                                                                                                                                                                                                                gtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCCAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTAAAATCCAGCATTTTTCGCTTCCCGAAGCTGt                     >  2:492509/1‑274 (MQ=255)
                                                                                                                                                                                                                                     cTGTCATTAACATGAGCCCCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTcccc                                                                                                                                                                                                    <  6:582081/94‑1 (MQ=255)
                                                                                                                                                                                                                                     cTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTcccc                                                                                                                                                                                                    >  5:582081/1‑94 (MQ=255)
                                                                                                                                                                                                                                         cATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTAAAATCCAGCATTTTTCGCTTACCGAAGCTGTAACTTTCCtt           >  2:454875/1‑275 (MQ=255)
                                                                                                                                                                                                                                           ttAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                                                                                                                                                                                                                            >  6:10480/1‑48 (MQ=255)
                                                                                                                                                                                                                                           ttAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                                                                                                                                                                                                                            <  5:10480/48‑1 (MQ=255)
                                                                                                                                                                                                                                            tAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTcag                                                                                                                                                                                                                                            >  7:384171/1‑45 (MQ=255)
                                                                                                                                                                                                                                            tAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCgg                                                                                                                                                                                                                                            <  8:384171/47‑1 (MQ=255)
                                                                                                                                                                                                                                            tAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGa                                                                                                                                             <  6:348898/142‑1 (MQ=255)
                                                                                                                                                                                                                                            tAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGa                                                                                                                                             >  5:348898/1‑142 (MQ=255)
                                                                                                                                                                                                                                               cATGAGCACCGTTTTCTCCCTCTCCCTTCCGGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                                                                                                                               <  8:58004/137‑1 (MQ=255)
                                                                                                                                                                                                                                               cATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCtttttt                                                                                                                                               >  7:58004/1‑137 (MQ=255)
                                                                                                                                                                                                                                                  gAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTAAAATCCAGCATTTTTCGCt                                 >  2:152519/1‑244 (MQ=255)
                                                                                                                                                                                                                                                  gAGCACCGTTTTCTCCCGCTCCCTTCCAGGGAGAGGGTCGGGGGGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTAAAATCCAGCATTTTTCGCTTCCCGAAGCTGTAACTTTCCTTATACTCGAc  <  2:492716/275‑1 (MQ=255)
                                                                                                                                                                                                                                                    gCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                                                                                                                                               <  7:185368/84‑1 (MQ=255)
                                                                                                                                                                                                                                                    gCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAcc                                                                                                                                                                                               >  8:185368/1‑84 (MQ=255)
                                                                                                                                                                                                                                                      aCCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTg                                                                                                                             >  5:310786/1‑148 (MQ=255)
                                                                                                                                                                                                                                                      aCCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTAAAATCCAGCATTTTTCGCTTccc                             <  4:152518/244‑1 (MQ=255)
                                                                                                                                                                                                                                                       ccGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGa                                                                                                                                                         <  7:788292/119‑1 (MQ=255)
                                                                                                                                                                                                                                                       ccGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGa                                                                                                                                                         >  8:788292/1‑119 (MQ=255)
                                                                                                                                                                                                                                                                ccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGc                                                                                                                                                                                                                 >  7:22262/1‑54 (MQ=255)
                                                                                                                                                                                                                                                                ccctctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGc                                                                                                                                                                                                                 <  8:22262/54‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAc                                                                                                                                                                                          <  8:1051634/75‑1 (MQ=21)
                                                                                                                                                                                                                                                                  ctctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAc                                                                                                                                                                                          >  7:1051634/1‑75 (MQ=21)
                                                                                                                                                                                                                                                                   tctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCctct                                                                                                                                                                                     >  6:615712/1‑79 (MQ=255)
                                                                                                                                                                                                                                                                   tctccctTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCctct                                                                                                                                                                                     <  5:615712/79‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tctccctGGAAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGa                                                                                                                                                         <  6:84828/107‑1 (MQ=255)
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          
TGGCGGGCGTAGCAGAGCGTTTCTCGATGCTCAAATATGGCCTGGCGGTGATTCTGGTGTTTATCGGTATCAAGATGCTGATTGTCGACTTCTACCATATTCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTAAAATCCAGCATTTTTCGCTTCCCGAAGCTGTAACTTTCCTTATACTCGAC  >  NC_000913/3239343‑3239859

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