Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,628,190 T→A A292A (GCT→GCA 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,1900TA94.4% 193.8 / NA 72A292A (GCT→GCAnadRDNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (34/34);  minor base C (3/1);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGG  >  NC_000913/4627943‑4628434
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
aGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAAGCGCGCTGGTCCACGCACaa                                                                                                                                                                                                                                                   <  1:543910/251‑1 (MQ=255)
     cgcgAAAACCCGTTCCGCTACGGGGAATATATTCCTACGGAAGGGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAAGACCACCAGTGCGTGGGAATATGGGCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACAt                                                                                                                                                                                                                                              <  1:598537/251‑1 (MQ=255)
       cgAAAACCCGTTTCGCTACTGGGAATATATTCCTGCGGAAGTGAAGCCGTTCTTTGTACGGACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGAGTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCAGTGCAGTAGTCCGATTTCGATAAAAGCGCGCTGGGCCACGCACAATACATTg                                                                                                                                                                                                                                            <  2:61642/251‑1 (MQ=255)
        gAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTgt                                                                                                                                                                                                                                           >  1:375723/1‑250 (MQ=255)
         aaaaCCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTTCGATTACGATAAAATCGCGCTGGGCCCCGCACAATACAtttat                                                                                                                                                                                                                                          >  1:429332/1‑248 (MQ=255)
           aaCCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGAttt                                                                                                                                                                                                                                        <  1:96748/251‑1 (MQ=255)
                             gAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCa                                                                                                                                                                                                                       >  1:184955/1‑250 (MQ=255)
                                 atatTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAAGACATTGATTTTGCAGTGAAATATGCCAATaaa                                                                                                                                                                                                                  <  2:601712/251‑1 (MQ=255)
                                     tCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTgg                                                                                                                                                                                                              >  1:569248/1‑251 (MQ=255)
                                      ccTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATTCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGAGTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGc                                                                                                                                                                                                             <  1:53132/251‑1 (MQ=255)
                                         acggAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACTCACAATACATTTATTTTGCAGTGAAATATGCCAATAAAGTGGcgtt                                                                                                                                                                                                          <  1:230320/248‑4 (MQ=255)
                                         acggAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAAACGCGCTGGGCCACGCACAATACATTTGTTTTGCAGGGAAATATGCCACTAAAggtgcggt                                                                                                                                                                                                          >  2:14936/4‑244 (MQ=255)
                                              aGTGACGCCGTTCTTTGTACGAACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCACTATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACAACTCGGCGGGGATGGAATCGCAATGCAGTATTCCGATTCCGATAAAATCGCGACGGGCCAAGACCAATACCTTGATTTTGCAATGAAATAAGCCacaacaggggcgtttattt                                                                                                                                                                                                     >  2:603970/1‑233 (MQ=255)
                                              aGTGAAGCCGTTCTTTGTACGTACCGTGGCGAGCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGGAAACAAACGTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTAttg                                                                                                                                                                                                     <  1:22767/251‑3 (MQ=255)
                                                     ccGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGCTTACGGTAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTAttgaaactt                                                                                                                                                                                               >  1:504760/1‑242 (MQ=255)
                                                       gttctttGTCCGTACCGAGGCGATACTTGGTGGCGAGTCGCGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATCCCACCAGTGCGTGGGAATATGGTCGCGCTTATGTCTTTTCACACCTCGGCGGTGATGAACTCGCCATGCAGTATTTCGATTACGATCAAAACGCGCTGGGCCACGAAACATACATTGATTTTGAAGTGAAAAATTCCAATAAAGTGGCGTTTAttgtcaatgatt                                                                                                                                                                                            >  2:308118/5‑240 (MQ=255)
                                                         tctttGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGAttt                                                                                                                                                                                           <  2:85747/248‑1 (MQ=255)
                                                              gtacgttccgaggCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACCATACATTGAGTTTGCAGTGAAATATTCCAATAAAGTGCCGTTTATTGACCCTGATTTTGTct                                                                                                                                                                                      <  2:368552/239‑2 (MQ=255)
                                                                 cGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTcactac                                                                                                                                                                                  >  1:198565/1‑248 (MQ=255)
                                                                    aCCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTcactaccc                                                                                                                                                                                <  1:475905/250‑6 (MQ=255)
                                                                      cgtcgCGAACCTTGGTGGCGAGTCGAACGGTCAATCCACCCTGGTAAACAAACTTTTCAATATCTTCAATACCAACATTTCGTGGAAATATGGTCGCTATTATGTCTTTTCACACCTCGGCGGTGATGATATCGCATTGAAGTATTCCGACTACGATAAAATCGCGCTGGGCCACGCACAATACACTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTctctaccctgg                                                                                                                                                                             <  2:485653/247‑11 (MQ=255)
                                                                               ccTTGGTGGCGAGTCGAGCGGTAGATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTGCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGACTACGATAAAATCTCGCTGGGCCCCGCCCAATTCATTTATTTTGCAGTGAAATATGCCATTAAAGTGGCTTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCa                                                                                                                                                                    <  2:381853/251‑1 (MQ=255)
                                                                                   ggtggcgAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCaaaa                                                                                                                                                                 <  1:134908/247‑1 (MQ=255)
                                                                                    gtggcgAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGGGGCGTTTATTTACACTGATTTTGTCACTACCCAGGCGTTCTGCaaaaa                                                                                                                                                                >  2:189226/3‑250 (MQ=255)
                                                                                           gTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAg                                                                                                                                                        <  2:161208/251‑1 (MQ=255)
                                                                                             cGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAggg                                                                                                                                                      <  2:429332/251‑1 (MQ=255)
                                                                                              gAGCGGTACATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGGGATGAGATCGCATTGCAGTATTTCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGAAAATGATTTTTTCACTACCCAGGCGCTCTGCAAAAAGTACGAAgggg                                                                                                                                                     >  2:157322/1‑250 (MQ=255)
                                                                                              gAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCCTTGCAGTATTTCGTTTTCGATAAAATCGCTCTGGGCCACGCCCAATACATTTATTTTGCAGTGAAATATGCCAATTAAGTGGCGTTTATTTACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAATTACGAAGGGc                                                                                                                                                     >  2:416604/1‑251 (MQ=255)
                                                                                               aGCGGTAAAGCCGCCCTGGTAAGCAAACTTGCCAATATCTTCAATCCCACCAGTGCGTTGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTTCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGAGTTTGTCACTACCCAGGCGTTCTGCAAAAAGGACTAAgg                                                                                                                                                       <  2:251242/248‑1 (MQ=255)
                                                                                                  ggTAAATCCACCCTGGTAAACAACCTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTAAGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGGGAAATATGCCAATAAAGTGGCGTTTATTGACACCGATTTTTTCACTAACCAGGCGTTCTGCAAAAAGAACGAAGGGCgtg                                                                                                                                                  >  1:171481/1‑248 (MQ=255)
                                                                                                  ggTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCgtg                                                                                                                                                  >  1:266671/1‑248 (MQ=255)
                                                                                                   gTAAATCCTCCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTTCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACTATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATGGACACTGATTTTGGCACTACCCAGGCGTTCTGCAAGAAGTACGAAGGTCGGGag                                                                                                                                                <  2:603462/251‑2 (MQ=255)
                                                                                                   gTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGAAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGAAAATGATTTTGTCACTACCCAGGGGTTCTGCAAAAAGTACGAagggaggg                                                                                                                                                  >  1:72499/1‑249 (MQ=255)
                                                                                                       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                                                                                                                        >  1:6614/1‑181 (MQ=255)
                                                                                                                                                                                                 tgatgaGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTtgatt                                                                                                                        <  2:6614/181‑1 (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)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
AGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGG  >  NC_000913/4627943‑4628434

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