Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,627,539 A→G P75P (CCA→CCG)  nadR → DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,627,5390AG91.3% 112.6 / 2.9 46P75P (CCA→CCG) nadRDNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (16/26);  total (18/28)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.61e-01

AGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTAGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCA‑GCCGACCGTG‑CCGGATCGTCTGCGTTGGT‑TATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGG  >  NC_000913/4627307‑4627785
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         
agggagATATGTAGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACatatt                                                                                                                                                                                                                                         >  1:407811/1‑250 (MQ=255)
    agaTATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGACTTGCGTTTCCCCGGCAGAAGAAAACCATTGGGGTGGTATTCGGTCATTTTTACCCGCTGCAAACCGGACATATCTAcc                                                                                                                                                                                                                                     >  2:279793/1‑251 (MQ=255)
           tcgtcATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTAACTTATcc                                                                                                                                                                                                                               >  1:493587/1‑250 (MQ=255)
                 tttGATTTCCTGACAACTGCCATCAAACAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGAGATTTAAGCCAGTTACTGAATGCCAAAATCATAAGCCCCAGCGCGCAAAAGCGGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTGTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATAGCTGCCTTATCCAgcgct                                                                                                                                                                                                                         <  2:14946/250‑2 (MQ=255)
                         ccTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGCATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGcc                                                                                                                                                                                                                >  2:20338/1‑251 (MQ=255)
                                   gCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGa                                                                                                                                                                                                      <  1:14300/251‑1 (MQ=255)
                                     cATCAAGCAACAGGGCTGCACGCTACAGCAGGCAGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCCAGCGCGCCTATAGCCAGGTTGACGAGc                                                                                                                                                                                                    >  2:101090/1‑251 (MQ=255)
                                             aaCAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATc                                                                                                                                                                                            <  1:437141/251‑1 (MQ=255)
                                               aaGGGCTGCACGCAACAGCAGGTGGCTGATGCCAGCGGGATGACCAAAGGGACTTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGAAAAAGCTGGAGGCGTTGAACCGTTTTTTGAGGCTTGAGTTTCCACGGCCGAAGTAAACCATTGGTGTGCTATTCGGAAAATTCTACCCGCTGCTTACCGGACATATCTACCTTATCCAGCGCGCCTGTCGCCAGGCAGACGAACTGAATATCAt                                                                                                                                                                                          >  2:268277/2‑251 (MQ=255)
                                                               agcagGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGAt                                                                                                                                                                          >  1:414271/1‑251 (MQ=255)
                                                                    gTAGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCg                                                                                                                                                                     >  1:139696/1‑251 (MQ=255)
                                                                      ggCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTGGACGAACTGCATATCATTATGTGTTTTGACGATACCCgc                                                                                                                                                                    <  1:135525/249‑2 (MQ=255)
                                                                             gCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATcgcgc                                                                                                                                                            <  1:379004/251‑1 (MQ=255)
                                                                                    gTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCGGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCgttgttt                                                                                                                                                     <  2:208033/251‑2 (MQ=255)
                                                                                      aTGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCg                                                                                                                                                    <  2:139696/250‑1 (MQ=255)
                                                                                        gACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAg                                                                                                                                                 <  1:273660/251‑1 (MQ=255)
                                                                                               gggTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGAGCGTTGTTTGAAGACAgcgc                                                                                                                                          >  1:326613/1‑248 (MQ=255)
                                                                                                    tttAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCTAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTAt                                                                                                                                                                                       >  2:331888/1‑201 (MQ=255)
                                                                                                    tttAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTAt                                                                                                                                                                                       <  1:331888/201‑1 (MQ=255)
                                                                                                     ttAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGt                                                                                                                                     <  1:142371/250‑1 (MQ=255)
                                                                                                       aaGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCgc                                                                                                                                  >  1:254648/1‑251 (MQ=255)
                                                                                                                    aaTGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCGGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGATCTACCCGCGGCATACCGGACATAGCTACCTTATCCAGCGCGCCTGTAGCCAGGGTGACGAACTGCATATCAGTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCAACCGt                                                                                                                    <  2:98429/251‑1 (MQ=255)
                                                                                                                       gcccaaaTCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGGTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCCCCCGTG‑cc                                                                                                                <  2:450667/247‑1 (MQ=255)
                                                                                                                        ccAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCCACCGTG‑CCg                                                                                                               <  1:25249/251‑1 (MQ=255)
                                                                                                                        ccAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCCACCGGG‑CCg                                                                                                               >  2:12750/1‑251 (MQ=255)
                                                                                                                             aTCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCCAGCGCGCCTATAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCA‑GCCGACCGTG‑CCGGATCg                                                                                                          <  1:101090/251‑1 (MQ=255)
                                                                                                                                   aGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCCGGTTGACGCACTGCATATCATTCTGGGGTTTGACGAAACCCGCGATCGCGCGGTGTTTGAAGACAGCGCTATGTCGCATCA‑GCC‑ACCCTGGCCGGATCGTcgggg                                                                                                    >  2:442423/1‑247 (MQ=255)
                                                                                                                                      cccAGCGCGCAAAAGCGGGAGGCGTTGCGCCGGTTTTTGGGGCTTGAGTTTCCCCGGACGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCCACCGTG‑CCGGATCGTCTGCgctg                                                                                                 <  2:136758/251‑4 (MQ=255)
                                                                                                                                              gcAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCCACCGTG‑CCGGATCGTCTGCGCTGGT‑TAttac                                                                                        <  1:215229/251‑3 (MQ=255)
                                                                                                                                               cAAAAGCTGGAGGCGTTGCACCGTTTGTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCCACCGTG‑CCGGATCGTCTGCGCTGGT‑TAttaca                                                                                       <  2:6650/251‑4 (MQ=255)
                                                                                                                                                  aaGCTGGAGCCGTTGCACCGCTTTTTGGGGCTTGAGTTTCCCCGGGAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTTCCCGCTGCATACCGGTCATATCTACCTTATCCAGCGCGCCTGTAGCCTGTTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCTCGATCGCGCGTTGTTTGAAGACATCGCTATGTCGCAGCA‑TCCCACCGTG‑CCGGATCGTCTGCGCTGGT‑TATTACAAAc                                                                                    <  1:486055/251‑1 (MQ=255)
                                                                                                                                                     cTGGCGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCA‑CACGCCCACCGTG‑CCGGATCGGCTGCGCTGGT‑TATTACAAACttt                                                                                 >  2:92958/1‑251 (MQ=255)
                                                                                                                                                          ggCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCCACCGTG‑CCGGATCGTCTGCGCTGGT‑TATTACAAACTTTTAAat                                                                            >  1:502069/1‑251 (MQ=255)
                                                                                                                                                             gTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCCACCGTG‑CCGGATCGTCTGCGCTGGG‑TATTACAAACTTTTAAATATc                                                                         >  1:475579/1‑251 (MQ=255)
                                                                                                                                                                gCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCCACCGTG‑CCGGATCGTCTGCGCTGGT‑TATTACAAACTTTTAAATATCAGa                                                                      <  2:90311/251‑1 (MQ=255)
                                                                                                                                                                     gTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTGTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCCACCGTG‑CCGGATCGTCTGCGCTGGT‑TATTACAAACTTTTAAATATCAGAAAAat                                                                 <  1:20338/251‑1 (MQ=255)
                                                                                                                                                                            ggggTTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCCTCCGTG‑CCGGATCGTCTGCGCTGGT‑TATTACAAACTTTTAAATATCAGAAAAATATTCGCa                                                          <  2:94924/251‑1 (MQ=255)
                                                                                                                                                                             gggCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCCACCGTG‑CCGGATCATCTGCGCTGGT‑TATTACAAACTTTTAAATATCAGAAAAATATTCGCAt                                                         <  1:244186/251‑1 (MQ=255)
                                                                                                                                                                              ggCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATAGCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCCACCGTG‑CCGGATCGTCTGCGCTGGT‑TATTACAAACTTTTAAATATCAGAAAAATATTCGCAtt                                                        <  2:254648/251‑1 (MQ=255)
                                                                                                                                                                               gCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCCACCGTG‑CCGGATCGTCTGCGCTGGT‑TATTACAAACTTTTAAATATCAGAAAAATATTCGCATTc                                                       <  1:235003/251‑1 (MQ=255)
                                                                                                                                                                                   gAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATTTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATAGCATTATGGGTTTTGACGATTCCCGCGATCGCTCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCCACCGTG‑CCGGATCGTCTGCGCTGGT‑TATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATg                                                    <  1:464361/250‑1 (MQ=255)
                                                                                                                                                                                              gCAGAAGAAAACCATTGGTGTGGTATTCGGTAATTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTATCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCCACCGTG‑CCGGATCGTCTGCGCTGGT‑TATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAag                                        <  1:506993/251‑1 (MQ=255)
                                                                                                                                                                                                      aaaccattggTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGCCATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCCACCGTG‑CCGGATCGTCTGCGCTGGT‑TATTACAAACTTTTAAATATCAGAAAAATATTCGCATTAATGCTTTCAACGAAGAGGGCATg                                >  2:354888/9‑251 (MQ=255)
                                                                                                                                                                                                         ccattggTGTGGTATTCTGTAAGTTCTACCCGCTTCATCCCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATAGCAGTATGTGTTTTGACGATCCCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCCACCGTG‑CCGGACCGTCTGCGCTGGT‑TATTACAAACTTTTAAATATCAGAAAAATATGCGCATTCATGCTTTCAACGAAGAGGGCATGGAg                             <  1:90778/246‑1 (MQ=255)
                                                                                                                                                                                                                    gTATTCGGTAAGTTCTACCCGTTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCA‑GCCCACCGTG‑CCGGATCGTCTGCGCTGGT‑TATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCa                  <  1:430250/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aCCCGCTGCATACCGGACATATCAACCTTAACCAGCGCGCCTGTACACAGGCTGACGACCTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGACCACAGCGCTATGTCGCAGCC‑GCCCACCGTG‑CCGGATCGTCTGCGGTGGTCTAT‑AAAAACTTTTAAAAATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGAATCCGCAGGGCTGGCGTGTGTgg  >  2:68995/1‑251 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         
AGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTAGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCA‑GCCGACCGTG‑CCGGATCGTCTGCGTTGGT‑TATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGG  >  NC_000913/4627307‑4627785

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑