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.5% 207.0 / NA 66S683Y (TCC→TAC) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (37/28);  total (38/28)
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

TAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAGACTGT  >  NC_000913/2145052‑2145581
                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                            
tAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCCGCTGAGCTACAGCGCCTACCATGATGCACt                                                                                                                                                                                                                                                                  >  2:31105/1‑274 (MQ=255)
           cTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTACCGCATCTCGc                                                                                                                                                                                                                                                      >  4:67563/1‑275 (MQ=255)
                gACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGACTACCATGATGCACGGACGCATCTCGCCTAcg                                                                                                                                                                                                                                                 >  2:247485/1‑274 (MQ=255)
                   gtgtttGCATATTACTTTTTGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCCCGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGAATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACcgcg                                                                                                                                                                                                                                              <  4:17966/271‑1 (MQ=255)
                               tACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGGCCCAGCTGCGCTACAGAGCCTACCATGATGCACTGACGCATCTCCCCAAACTCGCCAGCTATgag                                                                                                                                                                                                                                   >  1:80521/1‑274 (MQ=255)
                               tACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTTCAGCTCCTACTATGATGCACTCACGAATCTCGCAAACCGCCCCAGCTTTgag                                                                                                                                                                                                                                   >  1:8902/1‑274 (MQ=255)
                                                    agtAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGGTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTGCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCtgctgc                                                                                                                                                                                                             <  1:186657/273‑1 (MQ=255)
                                                    aTTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCtgctgc                                                                                                                                                                                                             <  1:183542/275‑1 (MQ=255)
                                                           gagaACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGt                                                                                                                                                                                                      >  1:159762/1‑275 (MQ=255)
                                                             gaACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTTTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTAAGGACGTTCCCGAATCAAGCAAAATGCTACGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGATTTGAGAACCAACTTCGTAACCTGCTGCCAACGGCaa                                                                                                                                                                                                    >  3:191437/1‑275 (MQ=255)
                                                             gaACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGACTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCCCTGACCCACCTCCCCAACCCCGCCCGTTTTGAGCAACAACTTCGTACCCTGCTGCCACCGGTa                                                                                                                                                                                                     >  1:56021/1‑274 (MQ=255)
                                                               aactgtgcAGTGGCGACGACACACGTTCCGCGATCGAGCAAGATGTGGTTTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATAGTGGTTCGGTTCTGGTGATTCATGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTaaa                                                                                                                                                                                                   <  2:202048/267‑1 (MQ=255)
                                                                   ttACAGTGCCGCCACCTCACGTTCCGCGATCGAGCAAGATGGGGTGTTGAACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTa                                                                                                                                                                                              <  2:33579/275‑1 (MQ=255)
                                                                        gTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTCCGCCAGCTGAGCTACAGCGCCGACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTTCTGCAACAGCTAAACAGTaca                                                                                                                                                                                            >  3:131753/1‑272 (MQ=255)
                                                                           ccGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAACCAAACGGGTATCCTGCTGAAAACGGTAAACAGCACACATCAg                                                                                                                                                                                      >  3:265974/1‑275 (MQ=255)
                                                                           ccGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCGTGCTGCAAACGGTAAACAGTACACATCAg                                                                                                                                                                                      >  2:225139/1‑275 (MQ=255)
                                                                              acacCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCatg                                                                                                                                                                                                                                                                         <  4:35047/189‑1 (MQ=255)
                                                                               cacCTCACGTTCCGGGTGCGGAGAAGTGGTGGTGTTGACATGTAGGGGCGGCGGGAGTTACGACTGGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGAc                                                                                                                                                                                  <  3:267064/275‑1 (MQ=255)
                                                                               aacCGCACGTTCGTCGATCGCGTAAGATGGGGTGTTGCACGGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGAAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGAc                                                                                                                                                                                  <  2:82528/274‑1 (MQ=255)
                                                                                   tCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACTGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACAATGATGCACAGACGCAGCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGATGCAAACGGTACACAGTAAACACCAGCGCCAtac                                                                                                                                                                              >  2:74764/1‑273 (MQ=255)
                                                                                   tCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCAAACCGAGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACAACTACGCGACATg                                                                                                                                                                               >  1:195929/1‑274 (MQ=255)
                                                                                          ccGCGATCGAGCAAGGAGGGGTGGTGCACTGTCTGTGCGTCGGCAGTTACGACGTGCATTACAGTATTGCGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGgt                                                                                                                                                                        <  3:152931/274‑1 (MQ=255)
                                                                                               gtgcaGCAAGGTGTGTGTGTGCACTGCCGGGGCGTCGGCCGTTACGAGGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAAGATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCCCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTAGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTAt                                                                                                                                                                  <  2:205688/271‑1 (MQ=255)
                                                                                                         gtgtgGTGTTGCACTTGCGGAGCGGCGGCAGTTACGACGTGAAGTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTGCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGAt                                                                                                                                                        <  4:239999/274‑1 (MQ=255)
                                                                                                                   tGCACTGTCGGGACGGCGGCAGTTACGACGTTCATTACAGTATTAAGCCGTTAAGTACTCTGGAAGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCCTTaa                                                                                                                                                <  3:195928/274‑1 (MQ=255)
                                                                                                                     cACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGCCGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGTCAACCGCGCCCGTGTTGAGAAATAACTGCGTATGCTGCTGCCAACGGTAAACAGTACACATCAGCGACATGAACTGGTGTTAATCGATCTTGATCCTTTTAAAg                                                                                                                                              >  2:265353/1‑274 (MQ=255)
                                                                                                                     cACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAg                                                                                                                                              >  3:170212/1‑274 (MQ=255)
                                                                                                                     cACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATAAGCGAAATGCCCTGGTGTTTATCGATCTTGACCGCTTTAAAg                                                                                                                                              >  3:30524/1‑274 (MQ=255)
                                                                                                                          tCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCCGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCACACCGCGCCAGTTTTGAGAAACAACTGCGTCTCCTGCTGCAACCGGAAAACAGTACACATCAGCGACATGCCCTGGTGTTGATCGATTTTGATCGGTTGAAAGAGGTGa                                                                                                                                        >  2:120528/1‑275 (MQ=255)
                                                                                                                          tCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCTTATCCTGCTGCAAACGGTAAACAGTACACATAAGCGACATGCCCTGGTGTTTATCGATCTTGCTCGCTTTAAAGCGGTg                                                                                                                                         >  2:237939/1‑274 (MQ=255)
                                                                                                                          tCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGAATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGcc                                                                                                                                                                             >  1:84836/1‑238 (MQ=255)
                                                                                                                              aGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGAATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTgg                                                                                                                                                                         <  3:84836/238‑1 (MQ=255)
                                                                                                                               gcggcgGCACGTACGACGGGCATTACATTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTAGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGAc                                                                                                                                   <  3:8902/275‑1 (MQ=255)
                                                                                                                                    ggcAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATCTTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCCCTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACCTGCCCTGGTGTTTATCGATCTTGATTGCTTTAAAGCGGTGACAGGCAgcgc                                                                                                                              >  4:82251/1‑275 (MQ=255)
                                                                                                                                    ggcAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCTCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTCAGAAACAACTGCGTATCCTGCTGCAAACGGTCAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTTCTCGCTTTAAAGCGGTGAATGACAgcgc                                                                                                                              >  2:107621/1‑275 (MQ=255)
                                                                                                                                    ggcAGTTACAACGTGAATTAAAGTTTTACGCCGTTAAGTACTCGTGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGAAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCAGCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAgcgc                                                                                                                              <  3:219372/275‑1 (MQ=255)
                                                                                                                                                 tGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCCGGTGTTCATCGACCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATgcc                                                                                                                    >  2:13830/1‑271 (MQ=255)
                                                                                                                                                        cAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACg                                                                                                          >  3:11524/1‑275 (MQ=255)
                                                                                                                                                           tATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTTAGAAACAACTCCGTATCCTGCTGCAAACGGTAACCAGTACACATCCGCGCCATTCCCTGGTGTTTATAGATCTTGATCGCTTTCAAGCGGTGAATGACAGCGCCTGGCATGCGGGCTGTTACGCtt                                                                                                       >  4:24409/1‑275 (MQ=255)
                                                                                                                                                                cGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCCAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAACGGTGAATGACAGCGCCGGGCATGCGGCGGGTGAGGCTAtgctt                                                                                                  >  1:113980/1‑274 (MQ=255)
                                                                                                                                                                  ccGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAAc                                                                                                                                                                                                                                                  >  2:97499/1‑129 (MQ=255)
                                                                                                                                                                   cGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGGGAATGACCGCGCCGGGCATGCGGCGGGTGCGGCTTtgcgtcgc                                                                                               >  3:104575/1‑270 (MQ=255)
                                                                                                                                                                     ttAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGCCGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCACACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGAGGTGAATGACAGCTCCTGGCATGCGGCGGGTGCCGCTTTGATGCGCGa                                                                                             >  2:147327/1‑275 (MQ=255)
                                                                                                                                                                     ttAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAACCGGTAAACAGTACACATCAGCGTAATGCCCTGGCGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTACTGCGCGa                                                                                             >  1:151087/1‑275 (MQ=255)
                                                                                                                                                                      tAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACcgcg                                                                                                                                                                                                                                              <  4:97499/129‑1 (MQ=255)
                                                                                                                                                                       aaGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACgg                                                                                                                                                                                                       >  3:221188/1‑167 (MQ=255)
                                                                                                                                                                           accctGCACGGAAGCAATATTGGTTCGGTTCTGGGGCTTCAGGACGTTACCGAATCACGCAAAAGGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGGCAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGc                                                                                       <  4:107621/272‑1 (MQ=255)
                                                                                                                                                                           aCTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTaaa                                                                                                                                                                                                   <  1:221189/167‑1 (MQ=255)
                                                                                                                                                                             tctGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTTTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGCCAACGCCGGGCATGCGGCGGGTTCGGTTTTGCTGGTCGGACTgggct                                                                                     >  4:187928/1‑272 (MQ=255)
                                                                                                                                                                               tGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGGGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACGGGcgtcg                                                                                   >  1:219546/1‑275 (MQ=255)
                                                                                                                                                                               tGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGGGCGAACTGGagtcg                                                                                   >  2:186501/1‑275 (MQ=255)
                                                                                                                                                                                  aaggaagaaaTATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTa                                                                                <  4:13830/267‑1 (MQ=255)
                                                                                                                                                                                    ggccggaaTCTTGGTTCGGTTCTGGTGATGCAGAGCGTTACCGAATCACGCAAAATGCTGCGCCCGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAAt                                                                              <  4:31105/269‑1 (MQ=255)
                                                                                                                                                                                         caATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTTACGCTTTGCTGCGCGAACTGGCGGCGTTAATGCTg                                                                          >  1:78615/1‑274 (MQ=255)
                                                                                                                                                                                            tatTGGTTCGGTTCTGGTGATTCAGGACGTTAACGAATCACGCAAAATGCTGCGCCAGCTGGGCTACAGCGCCTACCAGGATGCACTGACGCATCTCGCCAACAGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTa                                                                      <  1:131754/275‑1 (MQ=255)
                                                                                                                                                                                                 gttcggttcTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGCATGACAGCACCGGGCATGCGGCAGGTGACGCTTTGCTGCGCGAACCGGCGGCGTTAATGCCGAGTGtgctg                                                                 >  4:60012/1‑275 (MQ=255)
                                                                                                                                                                                                           ggTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCg                                                       <  3:50374/275‑1 (MQ=255)
                                                                                                                                                                                                                 tCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATTCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCTGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTCATGCTGATTATGCTGCGTTCCAGctgcgttc                                                 >  4:160077/1‑268 (MQ=255)
                                                                                                                                                                                                                                    cGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTggg                                >  1:220122/1‑273 (MQ=255)
                                                                                                                                                                                                                                              tgcGCCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGt                    <  3:220121/275‑1 (MQ=255)
                                                                                                                                                                                                                                               gcgcCAGCTGAGCTACAGCGCCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTGCCAGCGACGTCCTGGCGCGACTCGGGGGGGATGAAGTTGGTc                   >  1:55278/1‑275 (MQ=255)
                                                                                                                                                                                                                                                  ccAGCTGCGCTACAGCGCCTACCATGATGCACGGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTTTCCTGCTGCAAACGGTACACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTctgc                <  1:265975/275‑1 (MQ=255)
                                                                                                                                                                                                                                                  ccAGCTGCGCTACAACGCCTACCAGGATGCCCTGACGCATCTCGCCAACCGCGCCAGTGTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACCCATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGCACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTctgc                <  1:30524/275‑1 (MQ=255)
                                                                                                                                                                                                                                                  ccACCTGAGCTACAGCCCCTACCATGATGCACTGACGAATCTCGCCAACCGCGCCAGTTTTTAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTctgc                <  3:219545/275‑1 (MQ=255)
                                                                                                                                                                                                                                                             aCAGCGCCTACCATGATGCACTGCTGCATCTCGCCAGCCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACAGCAGCGACATGCCCTGGTGTGTAGCGAGCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAGa      <  1:209838/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                gcggCTACCATGATGCACTGACGCATCTCGCCAACCGCGCCAATTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAGACTGt  <  3:29044/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                            
TAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGCAGCAATATTGGTTCGGTTCTGGTGATTCAGGACGTTACCGAATCACGCAAAATGCTGCGCCAGCTGAGCTACAGCGCCTCCCATGATGCACTGACGCATCTCGCCAACCGCGCCAGTTTTGAGAAACAACTGCGTATCCTGCTGCAAACGGTAAACAGTACACATCAGCGACATGCCCTGGTGTTTATCGATCTTGATCGCTTTAAAGCGGTGAATGACAGCGCCGGGCATGCGGCGGGTGACGCTTTGCTGCGCGAACTGGCGTCGTTAATGCTGAGTATGCTGCGTTCCAGCGACGTGCTGGCGCGACTCGGTGGGGATGAATTTGGTCTGCTGCTACCAGACTGT  >  NC_000913/2145052‑2145581

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