Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,612,331 C→T 7.1% intergenic (+42/‑80) ytjA → / → yjjU uncharacterized protein/putative patatin‑like family phospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,612,3310CT 7.1% 209.4 / 6.6 71intergenic (+42/‑80)ytjA/yjjUuncharacterized protein/putative patatin‑like family phospholipase
Reads supporting (aligned to +/- strand):  ref base C (34/32);  new base T (2/3);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.19e-01

TCGCATTGAGTGATTTTCACGCTCATATCAAGCGGTGACATTAACTATGGTAAAGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGACGAATAGCTCTGGTGTGCGAAGGCGGCGGACAGCGTGGAATTTTCACGGCTGGCGTGCTGGATGAGTTTATGCGCGCGCAGTTTAATCCTTTCGATCTTTATCTCGG  >  NC_000913/4612064‑4612586
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                               
tCGCATTGAGTGATTTTCACGCTCATATCAAGCGGTGACATTAACTATGGTAAAGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTgg                                                                                                                                                                                                                                                          >  2:133722/1‑275 (MQ=255)
              tttCACGCTCATATCAAGCGGTGACATTAACTATGGTAAAGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTa                                                                                                                                                                                                                                            <  2:146337/275‑1 (MQ=255)
                  aCGCTCATATCAAGCGGTGACATTAACTATGGTAAAGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCNAAACTTGTCTTTGTCTTCGGGATTATTTTGTTCCTGGTGAGTTTGTTCATGGTCCGAAAACGACCCTAGATTTCAGTCCCNTAATACATTTTTATGCAAACCACATNCCGCGCACAGGCTCTTTCGGtttttgagt                                                                                                                                                                                                                                        >  3:103890/1‑270 (MQ=255)
                  aCGCTCATATCAAGCGGTGACATTAACTATGGTAAAGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATAGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGa                                                                                                                                                                                                                                        >  2:369331/1‑275 (MQ=255)
                              aGCGGTGACATTAACTATGGTAAAGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTTCAGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGACAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAGTACTAAATTGc                                                                                                                                                                                                                            >  2:267853/1‑275 (MQ=255)
                               gcgtggACATTAACTTAGTGAAAGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCGGGTTATCGCGTTAATCGCCGCCGTACTGGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGc                                                                                                                                                                                                                            <  4:240550/269‑1 (MQ=255)
                                cGGTGACATTAACTATGGTAAAGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAAGGACCCTCGATTTCAGACAATTAAGAAAATTTTATCCAAATCCAGTCCATCGGACTGGCTTTTGCGGTTTTAGCAAATAATAAATTACGt                                                                                                                                                                                                                          >  2:430288/1‑275 (MQ=255)
                                  gTGACATTAACTATGGTAAAGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTACTCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTa                                                                                                                                                                                                                        >  2:108928/1‑275 (MQ=255)
                                  gTGACATTAACTATGGTAAAGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGCTTTCAGACCATTACGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGt                                                                                                                                                                                                                          >  2:240551/1‑273 (MQ=255)
                                    gACATTAACTATGGTAAAGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTa                                                                                                                                                                                                                        <  4:133721/273‑1 (MQ=255)
                                                ggTAAAGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGc                                                                                                                                                                                                                            >  4:370482/1‑257 (MQ=255)
                                                    aaGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTa                                                                                                                                                                                                                        <  2:370483/257‑1 (MQ=255)
                                                    aaGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTcaca                                                                                                                                                                                                      >  4:262167/1‑275 (MQ=255)
                                                     aGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGATATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCCGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTCTGTTCATGGGCCGAAAACGACCCTAGATTTCAGCCCATAAAGAAAATTTTATCCACAGCCAGTCCAGCGCACTGGCCTTTCCGGTTTTCGCTAATATTAAATTGCGTTACTTTgccccttctcctta                                                                                                                                                                                                      >  4:483336/1‑261 (MQ=255)
                                                         gaCGCTTATGTTTCGTTGGGGCATCATACTTCTGATGATCGCGTTAATCGCCGCCGCGCTTGGGTTTGGTGATCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCg                                                                                                                                                                                                                           <  2:483337/249‑1 (MQ=255)
                                                          aCGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTAACCAAAGCCAGTCCAGCGGACTGGCTTTTGCCGGTTTAGCAAATAATAAATTGCGTTACTTTTTCATTTTTTCacaacc                                                                                                                                                                                                   >  2:435024/1‑271 (MQ=255)
                                                           tgtttATGTTTCGTTGGGGCTTCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACaaaaa                                                                                                                                                                                               <  4:108927/272‑1 (MQ=255)
                                                           cGCTTATGTTTCGTTGGGGCCTCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGATCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTGCTAGCGAATACTAAATTGCGCTACTTTGTTATTCTTTAACAACaaaca                                                                                                                                                                                               >  2:93627/1‑273 (MQ=255)
                                                           cGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTcacaacacaac                                                                                                                                                                                               >  1:59548/1‑274 (MQ=255)
                                                            gCTTATGTGTCTGTGGGGCATCATGTTTCTGGTTATCGCGTTGAGCGCCGCCGCGCTTGTGTTTGGTGGGCTGGCCGGTACCGCGTCAGGCGCAGCTAAAATTTGCTTTGTCGTCGGGATTATTCTGTTCCTGGGGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATGTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACaaaaa                                                                                                                                                                                               <  4:435023/274‑1 (MQ=255)
                                                                 tGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGaa                                                                                                                                                                                         >  2:154659/1‑275 (MQ=255)
                                                                      cGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTATGTTCCTGGTGAGTTTGTTCATGGGCCGAAAGCGACCCTAGCTTTCAGACCATTAAGAAAATTATATCCACAGCCAGTCAAGCGTAATGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCCTAACaaacaaatgcatgcag                                                                                                                                                                                    >  4:209256/1‑262 (MQ=255)
                                                                        ttGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTCTCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATACATTGCGTTACTTTTTCATTTTTTCACCACAAAAAAAGGAAGGCagag                                                                                                                                                                                  >  4:396298/1‑275 (MQ=255)
                                                                          ggggCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGa                                                                                                                                                                                                                                                              <  1:77681/197‑1 (MQ=255)
                                                                              caacgtagttCTGGTGTTGGCGTTGGTCGCCGGCGCACTTGGGTTTGGTGGTCGGGCCGGTATCGTTGCAGGCGCAGCTAAAATTGTTATTGTCGTTGGGATTATTATGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCTAGCGGACTGGCTTTTGCGGTTTTAGAGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGAAAGGCAGAGGTgggg                                                                                                                                                                            <  4:494703/267‑1 (MQ=255)
                                                                                        cTGGTTACCGCGTTAATCGCCGCGGCACTCGGGTTTGGTGTTCTGGCCGGTACCGCTGCAGGCGCGGCTAAAATTGTCGTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGTTTTCAGACCATTACGCACATCTTAGCCAAAGCCAGTCCAGCGGACTGGCTTGTGCGGTTTTAATGAATAATACATTTCGTTACTTTTTCAATCTTTCACAAATAAAAAGGGAAGGCAATGGAGGTGCAGCaaaatc                                                                                                                                                                  >  4:261285/1‑272 (MQ=255)
                                                                                                ggcgTTAATCGCCGACGGACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACg                                                                                                                                                          <  4:92127/274‑1 (MQ=255)
                                                                                                  cgTTAATGGCCGCCGCACTGTGGTTTGGTGGTCGGGCCGGGACCGCTGCAGGCGCAGCTAAAACGGTCTTTGTCGTCTGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACACTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTCGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCt                                                                                                                                                        <  4:56256/275‑1 (MQ=255)
                                                                                                  cgTTAATCGCCGCCGCACTTGTGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCt                                                                                                                                                        <  2:322045/275‑1 (MQ=255)
                                                                                                  cgTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCt                                                                                                                                                        <  4:370582/275‑1 (MQ=255)
                                                                                                     tAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTgg                                                                                                                                                     >  2:163413/1‑275 (MQ=255)
                                                                                                               cgcACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgcg                                                                                                                                           <  2:262168/275‑1 (MQ=255)
                                                                                                                    tttggTGTGGTGGTCTGGCCGGTACCGCTTAAGGCGCAGCTAAAATTGTCTTTGGTGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGtt                                                                                                                                      <  2:209257/272‑1 (MQ=255)
                                                                                                                    ttGGGTTTGGTGGTCTGGCCGGTACCGCTTCAGGCGAAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGtt                                                                                                                                      <  1:103890/275‑1 (MQ=255)
                                                                                                                      gggTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTg                                                                                                                                     >  4:294608/1‑274 (MQ=255)
                                                                                                                                    ggtggttACCGTTGCAGGGTAAGCTAAANTTTTTTTTGTCGTCGGGTTTTTTCANTTCCTGGTGAGTTTTTTCCTGTGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGNGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGtt                                                                                                                       <  3:89775/268‑1 (MQ=255)
                                                                                                                                                                                                ggTGAGTTTGTTCATGGGCCGAAAACGAACCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTCCTTTGTCATTCTTGCACACCAAAACGAGGAAGGCAGAGGTGCGGCAGCGAATACCTGTAACGCTTGGTAATCTTGCGCCGTTGACGCTACAGCCGTTCACGCCTGGACGATTCTCCCTGTTTTGCGAAGCAGGCGGACAGCAGGGAATTTTCAAGGCt                                                          >  4:39046/1‑275 (MQ=255)
                                                                                                                                                                                                      tttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGACGAATACCTCTGGTGTGCGAAGGCGGCGGACCGCGTGGAATTTTCCCGGCTGGCggc                                                    >  4:221818/1‑273 (MQ=255)
                                                                                                                                                                                                       ttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACt                                                                                                                                                                                                                                                            <  7:998669/74‑1 (MQ=255)
                                                                                                                                                                                                       ttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACt                                                                                                                                                                                                                                                            >  8:998669/1‑74 (MQ=255)
                                                                                                                                                                                                            cATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGACGAATAGCTCTGGTGTGCGAAGGCGGCGGACAGCGTGGAATTTTCACGGCTGGCGTGCTGGAt                                              <  2:261286/275‑1 (MQ=255)
                                                                                                                                                                                                            cATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCAAAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGACGAATAGCTCTGGTGTGCGAAGGCGGAGGACAGCGTGGAATTTTCACGGCTGGCGTGCTGGAt                                              <  4:154658/275‑1 (MQ=255)
                                                                                                                                                                                                                ggCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTTCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGACGAATAGCTATGGTGTGCGAAGGCGGCGGACGGCGTGGATTTTTAACGGCTGGCGTTCGGGAAGAGt                                          >  4:39658/1‑275 (MQ=255)
                                                                                                                                                                                                                    gAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGt                                                                                                                                                                                >  2:484197/1‑137 (MQ=255)
                                                                                                                                                                                                                        aCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTgggg                                                                                                                                                                            <  4:484196/137‑1 (MQ=255)
                                                                                                                                                                                                                           aCCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGTTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCGTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGTTTGGTAGTATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGACGAATAGCTCTGGTGTGCGAAGGCGGCGGACAGCGTGGAATTTTCACGGCTGGCGTGCTGGATGAGTTTATgcgcgcg                               <  2:288774/275‑1 (MQ=255)
                                                                                                                                                                                                                           aCCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGaataat                                                                                                                                                                                                                                   <  6:374439/79‑1 (MQ=255)
                                                                                                                                                                                                                           aCCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGaataat                                                                                                                                                                                                                                   >  5:374439/1‑79 (MQ=255)
                                                                                                                                                                                                                           aCCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGACGAATAGCTCTGGTGTGCGAAGGCGGCGGACAGCGTGGAATTTTCACGGCTGGCGTGCTGGATGAGTTTATgcgcgcg                               <  2:294609/275‑1 (MQ=255)
                                                                                                                                                                                                                           aCCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGACGAATAGCTCTGGTGTGCGAAGGCGGCGGACAGCGTGGAATTTTCACGGCTGGCGTGCTGGATGAGTTTATgcgcgcg                               <  4:163412/275‑1 (MQ=255)
                                                                                                                                                                                                                                   tttCAGACCATTAAGAAAAATTTATCCAAAGCCAGTCCACCTGACTGGCTTTTGCGGTGTTTGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGGGGGGCAGCGAATACCTGTAACGCGTGGAAATAGTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGACGAATAGCTCTGGTGTGCGAAGGCGGCGGACAGCGTGAAATTTTCACGGCTGGCGTGCTGGATGAGTTTATGCGCGCGCAGTTTaa                       <  4:93626/275‑1 (MQ=255)
                                                                                                                                                                                                                                      cAGACCATTAAGAATATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAACAAATTGCGTTACTTTGTCATTCTTTCATAACAAAAACAGGAAGTCAGAGGTGGGGCAGCGAATACCTGTAACGCCTGGTAATATTGCCCCGTTGTCGCTAAGGCCGTTCAAGCTTGGACGAATAGCTCTGGGGTGCGTAGGTGGCGGACAGCGTGGAATTTTCACGGCCGGCGTGCTGGAAGGGTTGATGCCCCCCCAGTTTAATc                     >  4:32541/1‑274 (MQ=255)
                                                                                                                                                                                                                                      cAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGc                                                                                                                                                                           >  2:533054/1‑124 (MQ=255)
                                                                                                                                                                                                                                       aGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGACGAATAGCTCTGGTGTGCGAAGGCGGCGGACAGCGTGGAATTTTCACGGCTGGCGTGCTGGATGAGTTTATGCGCGCGCAGTTTAATcc                    >  2:378593/1‑274 (MQ=255)
                                                                                                                                                                                                                                         aCCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGTCTTTTGCGGTTTAAGCGAATAATAAATTGCGTTAATTTGTGATTCTTTCACCACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCGGGAACGATCGGTAGTATTGCGCCGTTGCGACTACGGACGTCCGAGCCTGGCCGAACAGCTAGGGTGTGAGAAGGCAGCGGAGAGCGTGGATATATCACGGCTGGAGTTCGGGTTGGCGTTATgcgggcactaggaaagct                    >  4:117939/1‑260 (MQ=255)
                                                                                                                                                                                                                                          ccATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCg                                                                                                                                                                       <  4:533053/124‑1 (MQ=255)
                                                                                                                                                                                                                                              tAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGACGAATAGCTCTGGTGTGCGAAGGCGGCGGACAGCGTGGAATTTTCACGGCTgg                                                        >  2:413448/1‑231 (MQ=255)
                                                                                                                                                                                                                                                  aaaaTTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGACGAATAGCTCTGGTGTGCGAAGGCGGCGGACAGCGTGGAATTTTCACGGCTGGCGTg                                                    <  4:413447/231‑1 (MQ=255)
                                                                                                                                                                                                                                                        ttATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGGCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGACGAATAGCTCTGGTGTGCGAAGGCGGCGGACAGCGTGGAATTTTCACGGCTGGCGTGCTGGATGAGTTTATGCGCGCGCAGTTTAATCCTTTCGATCTTTATCTCgg  <  1:18210/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                aGCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCaga                                                                                                                                                                                   <  6:545504/90‑1 (MQ=255)
                                                                                                                                                                                                                                                                aGCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCaga                                                                                                                                                                                   >  5:545504/1‑90 (MQ=255)
                                                                                                                                                                                                                                                                 gCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTc                                                                                                                                                                                                         >  7:30047/1‑67 (MQ=38)
                                                                                                                                                                                                                                                                 gCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTc                                                                                                                                                                                                         <  8:30047/67‑1 (MQ=38)
                                                                                                                                                                                                                                                                  ccAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTa                                                                                                                                                                                                                        <  7:611400/51‑1 (MQ=25)
                                                                                                                                                                                                                                                                      tCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCt                                                                                                                                                        >  7:301529/1‑111 (MQ=255)
                                                                                                                                                                                                                                                                      tCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCt                                                                                                                                                        <  8:301529/111‑1 (MQ=255)
                                                                                                                                                                                                                                                                       ccAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCAt                                                                                                                                                                                                               >  7:632092/1‑55 (MQ=255)
                                                                                                                                                                                                                                                                       ccAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCAt                                                                                                                                                                                                               <  8:632092/55‑1 (MQ=255)
                                                                                                                                                                                                                                                                       ccAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGt                                                                                                                                       >  5:970037/1‑127 (MQ=255)
                                                                                                                                                                                                                                                                       ccAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGt                                                                                                                                       <  6:970037/127‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCag                                                                                                                                                                                    >  7:385319/1‑80 (MQ=255)
                                                                                                                                                                                                                                                                         aGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCag                                                                                                                                                                                    <  8:385319/80‑1 (MQ=255)
                                                                                                                                                                                                                                                                           cGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACa                                                                                                                                                                                             >  6:808140/1‑69 (MQ=255)
                                                                                                                                                                                                                                                                           cGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACa                                                                                                                                                                                             <  5:808140/69‑1 (MQ=255)
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                               
TCGCATTGAGTGATTTTCACGCTCATATCAAGCGGTGACATTAACTATGGTAAAGGAGACGCTTATGTTTCGTTGGGGCATCATATTTCTGGTTATCGCGTTAATCGCCGCCGCACTTGGGTTTGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTGGACGAATAGCTCTGGTGTGCGAAGGCGGCGGACAGCGTGGAATTTTCACGGCTGGCGTGCTGGATGAGTTTATGCGCGCGCAGTTTAATCCTTTCGATCTTTATCTCGG  >  NC_000913/4612064‑4612586

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