Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,628,292 A→G E326E (GAA→GAG 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,628,2920AG98.5% 191.6 / NA 65E326E (GAA→GAGnadRDNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (36/28);  minor base T (0/1);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

ATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGAT  >  NC_000913/4628046‑4628531
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
aTCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATc                                                                                                                                                                                                                                             >  1:244427/1‑251 (MQ=255)
         ggTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGt                                                                                                                                                                                                                                      <  2:69856/249‑1 (MQ=255)
                  acTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCCGTATTCCGGTTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGCTTTTGTCACTACCCAGGCGTTCTGCCAAAAGTACGAAGGGCGTGAGCATCCGTTCGtcaaggcgttg                                                                                                                                                                                                                            >  2:195227/1‑240 (MQ=255)
                                 caataccaccAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGGTGAGATCGCATTGCAGTATTCCGAATACGATAAAATCACGCTTGGCCACGCACCATACCTTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGGTTTTTGACACTAATTTTTTAACAACCCAGGCGTTCTGCAAAAAGTACAAAGGGAGCGAGCATCCGTTCGTACAGGCGTTGATTGGTGAATACCGtt                                                                                                                                                                                                            >  1:112635/5‑251 (MQ=255)
                                  aataccaccAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGttt                                                                                                                                                                                                           <  1:379554/248‑1 (MQ=255)
                                        accaGTGCGTGGGAATCTGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCCGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGCTTATTGACACTGATTTTGTCACTAACCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCt                                                                                                                                                                                                     >  1:531537/1‑251 (MQ=255)
                                             tGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGGGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATCAAGTGGCGTTTATTGACACTGATTTTGTCACAAACCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAAAAACGTTTCGATCTGGTg                                                                                                                                                                                                 >  1:424677/1‑250 (MQ=255)
                                                gTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATcc                                                                                                                                                                                             >  2:133185/1‑251 (MQ=255)
                                                   ggAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAATATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGTTTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGTTTTTGGCACTACCCAGGCGTTCTGCAAAAAGTACCAAGGGCGGGAGCATCCGTTCGGACAGGCGTTGATTGATGAATACCGTTTCGGTCTGGGGATCctcg                                                                                                                                                                                          >  2:186740/1‑249 (MQ=255)
                                                   ggAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGACTTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGAGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACCGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCctgc                                                                                                                                                                                          >  1:394516/1‑251 (MQ=255)
                                                      atatGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGATCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTgg                                                                                                                                                                                       >  1:46073/1‑251 (MQ=255)
                                                                gATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCAGGGCCACGCACAATACATTGATTTTGCAGGGAAAAATGCCAATAACGGGGCGTTTAATGACACTGATTTTGTCACTACCCAGGCGTTCTGAAAAAAGTACGGAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGCGATCCTGCTGGAGAACAacac                                                                                                                                                                             >  1:54388/1‑251 (MQ=255)
                                                                  ttATGTCTTTTCACACCTCGGCGGTGATGCGATCGCATTTCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGAAGTGAAATATGCCAATAAAGTGGGGTTTATTGAAAATGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGAACGAAGGGCGGGAGCATCCGTACGGAAAGGCGGTGTTTGGTGAATACAGTTTCGATCTGGTGATCCTGCTGGAGAACCACACGc                                                                                                                                                                           >  2:225840/1‑251 (MQ=255)
                                                                          tttCACACCTCGGCGGTTATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCGGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCgtgg                                                                                                                                                                      <  1:514479/248‑1 (MQ=255)
                                                                                   tCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCCCAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGGGGGGGGCGGATggg                                                                                                                                                          >  1:373615/1‑250 (MQ=255)
                                                                                    cggcggTGAGGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACTCCGTGGGTGGCGGATGGtt                                                                                                                                                         <  2:500478/251‑1 (MQ=255)
                                                                                     ggcggTGATGAGATCGCATTGCATTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGttt                                                                                                                                                        <  2:266671/251‑1 (MQ=255)
                                                                                     ggcggTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACCTTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGttt                                                                                                                                                        >  1:589436/1‑251 (MQ=255)
                                                                                       cggTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATTCATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTGATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGGATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTTGGTGGCGTATGGTTTAc                                                                                                                                                      <  2:105205/251‑1 (MQ=255)
                                                                                         gTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGc                                                                                                                                                    >  1:474966/1‑251 (MQ=255)
                                                                                          tgatgaGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTtgatt                                                                                                                                                                                                                         <  2:6614/181‑1 (MQ=255)
                                                                                          tgatgaGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTtgatt                                                                                                                                                                                                                         >  1:6614/1‑181 (MQ=255)
                                                                                              gagaTCGCATTGCAGTATGCCGATTACGATAAAATCGCGCTGGGCCACGCACAATTCATTTATTTTGCAGTGAAATATGCCAATAAAGTTGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCTTGAGCATCCGTTCGTACAGGCGTTGATTGAGGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGTCGGATGGTTTACTCAGCCt                                                                                                                                               <  2:388023/251‑1 (MQ=255)
                                                                                                   cGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCa                                                                                                                                          <  2:184955/251‑1 (MQ=255)
                                                                                                    gCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGc                                                                                                                                           >  1:164844/1‑249 (MQ=255)
                                                                                                     cATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTAGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGCGGGTGGCGGCTGGTTTACGCCGCCTCGGCAGt                                                                                                                                        >  2:464843/1‑251 (MQ=255)
                                                                                                      aTTGCAGTATTCCGATGACGATAAAATCGCGCTGGGCCACGCACAAGACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGtt                                                                                                                                       <  2:404129/251‑1 (MQ=255)
                                                                                                      aTTGCAGTAATCCGATTACGATAAACTCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGtt                                                                                                                                       >  2:328956/1‑251 (MQ=255)
                                                                                                       ttGCAGTATTCCGATTACGATACAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTc                                                                                                                                      >  1:160325/1‑251 (MQ=255)
                                                                                                        tGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCAAAATACATTGATTTTGCAGTGAAATATGCCAATAAAGAGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTATGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGAGGGAGGCGGATGGTTTACGCAGCCTCGGCAGTTCg                                                                                                                                     >  2:163207/1‑251 (MQ=255)
                                                                                                            gTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTgg                                                                                                                                 <  1:133185/251‑1 (MQ=255)
                                                                                                             gattccgattACGATAAAATCGCGCTGGGCCACGCACAATACATGGATTTTGCAGTGAAATAGGCCAATAAAGTGGCGGTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGa                                                                                                                                <  2:26752/242‑1 (MQ=255)
                                                                                                              attccgattACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGa                                                                                                                                <  2:244427/242‑1 (MQ=255)
                                                                                                                   gattACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATTCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTTATTGATGAATACCGTTTCGATCTGGTGATCCTTCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCTTCGGCAGTTCGGTGGATCGCaa                                                                                                                          <  1:423587/248‑1 (MQ=255)
                                                                                                                    attACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCaa                                                                                                                          >  1:103129/3‑250 (MQ=255)
                                                                                                                          aTAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTc                                                                                                                   >  1:341776/1‑251 (MQ=255)
                                                                                                                                    gcTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGc                                                                                                         <  1:218600/251‑1 (MQ=255)
                                                                                                                                       gggCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCGGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCtg                                                                                                       <  1:416604/250‑1 (MQ=255)
                                                                                                                                          ccACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCtggtgg                                                                                                   <  1:611019/251‑1 (MQ=255)
                                                                                                                                               cacAATACGTGGATTGTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGGTTCGATCTGGTGAGCCTGCTGGAGAACAACACGCCGTGGGTGGTGGAGGGTTTACGCAGCCGCGGCAGTTCGTGGGAGCGAAAAGAGTTCCAGAACTTGCGGGGGGAGATg                                                                                              <  2:102153/249‑1 (MQ=255)
                                                                                                                                                     ccATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTTCCCAGTCTTTCTGCAAAAAGTACGAAGGGCGTGATCATCCGTTCGTACAGGCGTTGATTGAGGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGGTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCtgga                                                                                         <  1:219161/249‑4 (MQ=255)
                                                                                                                                                      cATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAag                                                                                       <  2:136359/251‑1 (MQ=255)
                                                                                                                                                            ttttGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAAc                                                                                  >  2:44631/1‑250 (MQ=255)
                                                                                                                                                                      aaaTATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCg                                                                           >  1:475409/1‑247 (MQ=255)
                                                                                                                                                                      aaaTATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCg                                                                           <  2:475409/247‑1 (MQ=255)
                                                                                                                                                                       aaTATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAAttc                                                                      <  1:520622/251‑2 (MQ=255)
                                                                                                                                                                        atatGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCTGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAAttcg                                                                     <  1:603970/251‑3 (MQ=255)
                                                                                                                                                                         tatGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAAttcgt                                                                    >  1:188362/1‑248 (MQ=255)
                                                                                                                                                                          atGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGGACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTTGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATAGCGAATTcgtg                                                                   <  2:589436/251‑1 (MQ=255)
                                                                                                                                                                           tGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGc                                                                  >  1:442930/1‑251 (MQ=255)
                                                                                                                                                                            gCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCCCTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGc                                                                  <  2:375723/250‑1 (MQ=255)
                                                                                                                                                                              cAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCggg                                                               <  2:239909/251‑1 (MQ=255)
                                                                                                                                                                                      tGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAag                                                         >  1:163842/1‑249 (MQ=255)
                                                                                                                                                                                      tGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTAGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAag                                                         >  2:264838/1‑249 (MQ=255)
                                                                                                                                                                                       ggcgattatttacactgcttttGTCACTACCCAGGCGTTCTGCAGAAATTACGAAGGGCGTGAGCATCCCTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCAGCTGGAGAACAACCCGCCGTGGGTGTCGGATGGTTTACGCAGCCTCGACATCTCGGTGGATCGAAAAGAGTTCCTGAACTTGCTTGTGGAGGTTCTGGAAGAGCACAATCTCGAATTCGTGCGGGTTAAAGAgg                                                      >  2:59369/19‑251 (MQ=255)
                                                                                                                                                                                                      tgATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAgcc                                       >  2:238162/2‑249 (MQ=255)
                                                                                                                                                                                                       gATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGCATTCGTGCGGGTTGAAGAGGCAGATTATGACAgccg                                      >  2:80080/1‑248 (MQ=255)
                                                                                                                                                                                                               cactacccAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAAGACCGTTTCGATCTGGTGATCCGGCTGGAGAACAACACGCCGTGGGTGGCGGATGGGTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGTAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGAGTATGACAGCCGTTTCctg                               <  2:415687/243‑1 (MQ=255)
                                                                                                                                                                                                                actccccAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCTTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCTGTTTCctgcg                             <  2:394516/245‑1 (MQ=255)
                                                                                                                                                                                                                   acccAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCctgcgctg                          >  2:180934/4‑251 (MQ=255)
                                                                                                                                                                                                                    cccAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCctgcgctgc                         >  2:460295/3‑251 (MQ=255)
                                                                                                                                                                                                                     ccAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCctgcgctgcg                        <  1:53901/250‑1 (MQ=255)
                                                                                                                                                                                                                           gTTCTGCAAAAAGTACGAAAGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGagc                  >  2:385894/1‑249 (MQ=255)
                                                                                                                                                                                                                                 cAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGc            >  1:199824/1‑251 (MQ=255)
                                                                                                                                                                                                                                           gAAGGGCGTGAGAATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGCACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAgatgat  >  2:525027/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
ATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGAT  >  NC_000913/4628046‑4628531

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