Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,145,313 C→A S683Y (TCC→TAC)  yegE → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,145,3130CA98.2% 179.5 / NA 56S683Y (TCC→TAC) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (27/28);  minor base G (0/1);  total (27/29)
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.00e+00

AAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAGACTGTAATGTTGAAAGCGCGCGTTTTATCGC  >  NC_000913/2145021‑2145607
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                      
aaaaTGAGCGGCTGAACGCAGGAAGAAGCTTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCatgatg                                                                                                                                                                                                                                                                                                <  2:357836/298‑1 (MQ=255)
  gaTGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCa                                                                                                                                                                                                                                                                                              >  1:473659/1‑301 (MQ=255)
         ggtgggACGCAGGAAGAAGCTTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACGGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACg                                                                                                                                                                                                                                                                                        <  2:473649/296‑1 (MQ=255)
               aCGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAAGATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCCGCTGGCCTACAGCGCCTACCATGATGCACTGACGCCtctc                                                                                                                                                                                                                                                                                  >  2:126992/1‑300 (MQ=255)
                         atgCGTTCGTTGTGCCGCTCCTGATGGTGTTGCATATTGCTTTGGGCGACAACGGACCATTAATGGGGAACAGTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTTGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACcgcg                                                                                                                                                                                                                                                                        <  2:170229/298‑1 (MQ=255)
                           gCGTTAGGTGCTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCCCTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGTTTCGGTTCTGGTGATTCAGGACGTTACCGCATCACGCAAAATGCTGCGCCAGCTGATCTCCATCGCCTACCATTATGACCTGACGGAtcgaccaaaccgcccc                                                                                                                                                                                                                                                                      >  2:126819/1‑286 (MQ=255)
                                    gtgcggcgtgtggaggTGTGGCATATTACTTTTGGCGCAAACGGACCAGTAATGGAGAACTTTTACAGGGCCGAGGCCTCACGTGCCGCGGTCGGGCAAGGAGTGGTGTTGCCCTGTCGGAGCGGAGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTACTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTgag                                                                                                                                                                                                                                                             <  2:299079/286‑1 (MQ=255)
                                        cGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGTTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAaaca                                                                                                                                                                                                                                                        <  1:209560/301‑1 (MQ=255)
                                               gaaggggtTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCg                                                                                                                                                                                                                                                  <  2:198966/294‑1 (MQ=255)
                                               gaaggTGTGGCATATTACTTTGGGAGAAAACGGACCATTAATGGAGAAAATTTACAGGGCCGACACCTCACGTTCAGCGATCGAGCAAGATGTGTGGTTGCACTGGCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAAGATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCg                                                                                                                                                                                                                                                  <  2:342205/297‑1 (MQ=255)
                                                  tgtgttacaatgtacttgagataccAACGGAACATTCATGAAGAACAGTTACCGTGCCGCCACCGCCAGTTCCGCGATGGAGAAAGATGTGTTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCTTTACGGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGCCGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTAt                                                                                                                                                                                                                                               <  2:280287/276‑1 (MQ=255)
                                                     gtTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACCGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAATTGCTGCGCCATCTGAGCTCCCGCGCCTACCATGATGCACTGACGCATCTCGCCAACCACGCAAGTTTTGAGAAACAACTGCGTAGCct                                                                                                                                                                                                                                            >  2:42983/1‑300 (MQ=255)
                                                                     ggCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGGTTTGAGAAACACATGCGGATCCTGCTCCAAACGGAAAACa                                                                                                                                                                                                                           >  2:166954/1‑301 (MQ=255)
                                                                                     tAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCTAAATGCTGCGTCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCACTCCGCCAACCGCCCCAGTTTGTACAACCAACGTCGTATCATGCTGCACACGGCAAAAAGCAAACAtaaacacaa                                                                                                                                                                                                           >  2:190926/1‑293 (MQ=255)
                                                                                      aaTGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCACTATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACGCCCCCAGTTTTGAGCAACAACTGCCTCTCCTTCTGCAAACCGTAAACAGtcaaaaacatgtag                                                                                                                                                                                                            >  2:444605/1‑286 (MQ=255)
                                                                                             gacagtgaccgggCCGACACCCACCGTTCCGCGATGGCGCAGTGTGTGGTGTTGCCAGGGCGGGGCGGCGGCGTTTACGACGTGCATTACAGAATTCGGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTAGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGTGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCt                                                                                                                                                                                                     <  2:367304/287‑1 (MQ=255)
                                                                                               cATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTTGCCAACCGCGCAAGCTTTGAGAAACAAGTGCGTATCCTGCCGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGgtg                                                                                                                                                                                                 >  1:184348/1‑301 (MQ=255)
                                                                                                   tACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCCCCAGTTCTGAGAAACAACTGCGTATCCTGCCTCAAACGGTAAACAGTACACAACAGCGACATACCCTGGTGtttt                                                                                                                                                                                             >  2:174538/1‑300 (MQ=255)
                                                                                                           caccgccgcccattcgccgATCAGGCAAGGTGTGGTGTTGCGCTGTCGGGGAGACTGCGCTTACGACTTGCATTACCGGATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCtt                                                                                                                                                                                     <  2:184346/284‑1 (MQ=255)
                                                                                                                      gTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGACTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTACACAGTACACATCAGCGACATGCCCTGGTGTTTATCGACCTTGATCGCTTTaa                                                                                                                                                                          >  1:89640/1‑301 (MQ=255)
                                                                                                                       ttCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTAGCCAACCTCGCCAGTTTTGAGAAACCACTACCTAACCTCCTGCAAACGGTAAACAGTACACATCCGCCAACTCCCCTTGCGTTTATTCCTCTTGATCCCTTTaaa                                                                                                                                                                         >  2:262621/1‑301 (MQ=255)
                                                                                                                         ccGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACCGTACACATCAGCGACATGCCCTGGTGTTTATCGACCTTGATCGCTTTaaac                                                                                                                                                                        >  2:60809/1‑299 (MQ=255)
                                                                                                                              cgagggaacggtgtgATGTTGAACAGTGCGTGCGGAGGCCATTTCGACGTGAATTACGGTTATACACCGATAAGTCAGCTGGACGGCAGCAATCTTGGTTCGGTTCTGGGAAGTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCAGGGCTACAGCGCCGACCATGATGCAATGACGCATCTCGCCAACCGCGCTAGTGTTGAGAAACAAATGCGTATCCTGCTGCAAACGGTAAACGGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTg                                                                                                                                                                   <  2:519592/289‑1 (MQ=255)
                                                                                                                                cGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTTATAGCTTTAAAGCGGTGAAt                                                                                                                                                                >  1:119082/1‑301 (MQ=255)
                                                                                                                                         tgtgGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGcc                                                                                                                                                       <  1:469753/301‑1 (MQ=255)
                                                                                                                                           tggtgTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTACCCTGCTGCAAACGGTAAACAGTAAAAATCAGCCAAATGCCCTGGTGTTTATCGATCTTAATCGCTTTAAATGGGTGAATTACATCGCCgg                                                                                                                                                     >  2:455260/1‑301 (MQ=255)
                                                                                                                                            ggttgtGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGGATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCgg                                                                                                                                                     <  2:83837/295‑1 (MQ=255)
                                                                                                                                             gtgtTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCggg                                                                                                                                                    <  1:328443/300‑1 (MQ=255)
                                                                                                                                                          gtgggggcggcCGTAAAGACGTTGAATACATTATTACACCCTTAAGTACCCTGGAGCGCAGCAATATTGGTTGGGTTCTGGTGATTCAGGACGTTACCCAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCAGCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTg                                                                                                                                       <  2:248260/295‑1 (MQ=255)
                                                                                                                                                                      gTTACGACGTGCATTACAGTATTACGACGCTAAGTACTATGGATGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAACCACGCAAAATGCTGCGCTAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCCTCTCGCCAACCGCGCCAGCTTTCAGAAATAAATGCGTATCACGTTGCAAACAGTAAACAGTACAAAACAACGACAAGCGCTGGTGTTTACCGATATTGATCGCTTTAAAGCGGTGAAGCACAGCGCCGGTCATGCGGAGGGGGACGCTTTgcggcg                                                                                                                          >  2:466446/1‑301 (MQ=255)
                                                                                                                                                                       ttACGACGTGCATTACAGTATTACGCCGGTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTgcgc                                                                                                                         <  1:455270/301‑1 (MQ=255)
                                                                                                                                                                                         ttcagacgtgtgctcttccgatctcGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAAt                                                                                                        <  1:108512/276‑1 (MQ=255)
                                                                                                                                                                                            tACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGCATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCtt                                                                                                    >  1:381339/1‑300 (MQ=255)
                                                                                                                                                                                                  gTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCCACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATACTGCTGCAAACGGGAAACAGTACACATCAGCGACATTCCAAGGTGTTTATCGCTTCTCACGTCCTTAAAGCGGTGAATGACAGCGCCTGCCAACAGGCGGGTGATGGTTTCCGTCGCGCACTGGGGTTGTTACTCATGATTAtg                                                                                              >  2:293445/1‑301 (MQ=255)
                                                                                                                                                                                                         aCTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGtt                                                                                       >  1:148707/1‑301 (MQ=255)
                                                                                                                                                                                                           tctGGACGGCCGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTAAAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAACCGGTAAACAGTACCCATAAGCCACATGCCCGGCTGTTTATCCATCTTGATCGCTTTAAAGCGGTGCAGGACAGCGACGGACATGTGGCGGTGGACGCTTTGCTGCNGGGCCTGGCGTCGTTCATggctggtgtgctggcgtc                                                                                      >  2:70419/1‑283 (MQ=255)
                                                                                                                                                                                                            ctGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACAACAGCGACATGCCCTGGTGTTTATCGATCTTGCTCGCTTTAAAGCGGTGAATGACCGCGCCGGGCATGCGGCGGGAGACGCTTGTCTGCGCGAACACGCGTCGTTACTGCTGTGTATGCTGCGtgtca                                                                                    >  2:138121/1‑297 (MQ=255)
                                                                                                                                                                                                                 cGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCCCGAACTGGCGTCGTTattatgtcggacgagcgtcgtgtggggg                                                                               >  2:108510/1‑274 (MQ=255)
                                                                                                                                                                                                                   gcagcaATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCCGCGACGt                                                                             >  1:410595/1‑301 (MQ=255)
                                                                                                                                                                                                                            tgtgttcgagtcTGGTGATGAAGGACAGTACCGAAGCAGTCAACAGGCTGCGCCAGCTGAGCTACAGCGCCTACCATGGTGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAAAAACTGCGTGTCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGgcgc                                                                     <  2:89639/297‑1 (MQ=255)
                                                                                                                                                                                                                                    gttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTg                                                            >  1:387717/1‑301 (MQ=255)
                                                                                                                                                                                                                                        tGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGGCGCTAATCCTGCGTATGCTGCGTTCCCGCGACGTGCTGGCCCGACTCGggggggg                                                        >  2:470278/1‑300 (MQ=255)
                                                                                                                                                                                                                                        gcgtgattcggtaacgtcctgAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGGCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGggggggg                                                        >  1:79120/21‑300 (MQ=255)
                                                                                                                                                                                                                                          gTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTgggg                                                         >  1:135448/1‑298 (MQ=255)
                                                                                                                                                                                                                                                  ggACGTTACCGAATCACGCAAAATCCGCCGCAAGCTGACGTACAGCGCCTACCATGTTGCACTGACGCATCTCGCCAACCGGGCCAGTTTTGGGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTgg                                               <  2:148705/300‑1 (MQ=255)
                                                                                                                                                                                                                                                       ttACCGAATAACGAAAAATGCTTCGCCAGCTGAGCTCCAGCGCCGACAATGATGCACTGACGCAGCTCGACAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAGAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTctgc                                          <  2:410586/300‑1 (MQ=255)
                                                                                                                                                                                                                                                            gAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTa                                     >  1:270403/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                    cAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAGACTGt                            <  1:18878/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                           ctgcGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATTTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAGACTGTAATGTTg                     <  1:444615/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                           ctgcGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAGACTGTAATGTTg                     <  1:126994/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                             gcgcCAGATGGGGGACGGCGCCTACCATGATGAACTGACGCAGCTCGCCAACCGGGGCAGTTTTGAGAAACAACTGCGTATCCGGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAGACTGTAATGTTGa                    <  2:135446/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                              cgcCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAGACTGTAATGTTGaa                   <  1:367184/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    cTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAGACTGTAATGTTGAAAgcgcg             <  1:470288/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        gCTAAAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAGACTGTAATGTTGAAAGCGCGCGtt         <  1:262625/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         cTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAGACTGTAATGTTGAAAGCGCGCGttt        <  1:242162/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                               cgcCTGCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAGACTGTAATGTTGAAAGCGCGCGTTTTATCGc  <  1:70420/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                      
AAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAGACTGTAATGTTGAAAGCGCGCGTTTTATCGC  >  NC_000913/2145021‑2145607

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