Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,628,163 C→T D283D (GAC→GAT)  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,1630CT95.3% 118.7 / ‑3.8 43D283D (GAC→GAT) nadRDNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (19/22);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

TTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCA‑ATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCA‑GGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGG  >  NC_000913/4627921‑4628401
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
tttATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATTAGATCGCATTGCAGTATTCCGATTACGATa                                                                                                                                                                                                                                            >  1:148567/1‑250 (MQ=255)
   aTGAGTAGCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGTGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATaaaa                                                                                                                                                                                                                                         <  1:312347/250‑1 (MQ=255)
           cAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAAGAGATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACTCCTCGGCGGTGATGAGATCGCATTGGAGTATTCTGACTACGATAAAATCGCGCTgg                                                                                                                                                                                                                                <  1:137105/251‑1 (MQ=255)
               ggTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCa                                                                                                                                                                                                                            <  2:175525/251‑1 (MQ=255)
                              gAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACC‑ATACATTgt                                                                                                                                                                                                            >  1:21863/1‑250 (MQ=255)
                                                            ccTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCTAGCTGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTAGTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGc                                                                                                                                                                              <  1:642891/251‑1 (MQ=255)
                                                             cgacggaagggaatcggtgctttGTACGGACCGTGGCGATCCTTGGTGGCGAGTAGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGTAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGAGTTTGCAGTGAAATATGCCAATAAAGTGGcg                                                                                                                                                                             <  2:561247/231‑2 (MQ=255)
                                                                      tGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGGGATGAGATCGCATTGCAGGATTCCGATTACGATAAAATCGCGCTGGGCCAAGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGttttttgac                                                                                                                                                                    >  1:379039/1‑245 (MQ=255)
                                                                       gAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTAttgac                                                                                                                                                                    >  1:141713/1‑246 (MQ=255)
                                                                             gttctttGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTAttgacactgatt                                                                                                                                                             <  1:210179/247‑12 (MQ=255)
                                                                                    gtacGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAAGAAAGTGGCGTTTATTGACACTGATTTTGTca                                                                                                                                                       <  2:90310/247‑1 (MQ=255)
                                                                                             gTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGc                                                                                                                                             >  1:157008/1‑251 (MQ=255)
                                                                                                gCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTTGGCCACGC‑CACATACATTGTTTTTGCAGGGAAAAATGCCCATAAAGTGGGGTTTATTGACACGGATTTTGTCACTACCCACGCGtt                                                                                                                                          >  2:161335/1‑251 (MQ=255)
                                                                                                 cGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGGTTATTGACACTGATTTTGTCACTACCCAGGCggtt                                                                                                                                         >  1:140859/1‑248 (MQ=255)
                                                                                                  gATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGAAAATGATTTTGTCACTACCCAGGCGTTCt                                                                                                                                        >  1:451118/1‑251 (MQ=255)
                                                                                                              cgAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTGGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGGCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAAGCGCGCTGGGCCACGCACA‑ATACATTGATTTTGTAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTTCCCAGGCGTTCTGCAAAAAGCACg                                                                                                                            <  2:630423/251‑1 (MQ=255)
                                                                                                                aggcGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATACGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGaa                                                                                                                          <  1:247766/248‑1 (MQ=255)
                                                                                                                     aGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGGGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCg                                                                                                                     >  1:659709/1‑251 (MQ=255)
                                                                                                                                 aCCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGa                                                                                                                                                                                                            <  1:23969/152‑1 (MQ=255)
                                                                                                                                      ggTAAACAAACTTGCCAATATTTTCAATACCCCCAGTGCGGGGGAATATGGTCGCGATTATGTCTTGTCACACCGCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCTCGCACA‑ATACATTGATTTTGCAGTGAAAGATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTATCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGt                                                                                                      <  2:548739/249‑1 (MQ=255)
                                                                                                                                        tAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTTCGAAGGGCGTGAGCATCCGTTCGtac                                                                                                    >  1:513682/1‑247 (MQ=255)
                                                                                                                                         aaacaaacTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑g                                                                                                 <  1:355115/250‑1 (MQ=255)
                                                                                                                                           acaaacTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGTTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTTCCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑Ggcg                                                                                              <  2:451118/251‑1 (MQ=255)
                                                                                                                                           acaaacTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGGAGGGCGTGAGCATCCGTTCGTACA‑Ggcg                                                                                              >  1:361816/1‑251 (MQ=255)
                                                                                                                                              aacTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCA‑ATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCC‑GGCGCtg                                                                                           >  1:308665/1‑251 (MQ=255)
                                                                                                                                                  tGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTtgatt                                                                                        <  1:371793/250‑1 (MQ=255)
                                                                                                                                                  tGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTtgat                                                                                         <  1:336285/249‑1 (MQ=255)
                                                                                                                                                   gtcaataatctcaattccaccaGTGCGTGTGAATATGTTCGCGATTATGTCTTTGCACACCTCGGCGGTGATGAGATCGCATTGCAGTAGTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTtgattg                                                                                       <  1:595222/234‑1 (MQ=255)
                                                                                                                                                             ccaataccaccAGTGTGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCTGCGGTGATGAGATCGCATTGCAGTATTTCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATTCCGt                                                                            <  2:512062/246‑1 (MQ=255)
                                                                                                                                                                ataccaccAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATAAATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAAAACCGTTTc                                                                         >  1:506680/3‑251 (MQ=255)
                                                                                                                                                                 aaccaccAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACTATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCg                                                                        <  2:626983/250‑1 (MQ=255)
                                                                                                                                                                            cGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGGGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTTACACTGATTTTGTCACTAACCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCGGTCTGGGGATc                                                             >  1:479716/1‑251 (MQ=255)
                                                                                                                                                                               gggAATATGGGCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATTCCGTTTCGATCTGGTGATCc‑t                                                          <  1:671269/250‑1 (MQ=255)
                                                                                                                                                                                 gAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATcc                                                                                                            >  1:534478/1‑200 (MQ=255)
                                                                                                                                                                                 gAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATcc                                                                                                            <  2:534478/200‑1 (MQ=255)
                                                                                                                                                                                    tatGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGGATACCGTTTCg                                                                        >  1:334138/1‑232 (MQ=255)
                                                                                                                                                                                    tatGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCg                                                                        <  2:334138/232‑1 (MQ=255)
                                                                                                                                                                                    tatGGTCGCGATTATGTCTTTTCACACCTCGGAGGTGATGAGATCACATTGCCGTATTCCGATTACGATAAAATCGCGCTGGGCCACGAACA‑AAACATTGGTTTTGCAGGGAAATATGCCAATAAAGTGGCGATTATTGACACTGATTTTGTCACTAACCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAGCCTCCGTTCGTACA‑GGCGTTGATTGATTAATACCGTTTAGATCTGGTGATAC‑TGCTggg                                                    >  1:353187/1‑250 (MQ=255)
                                                                                                                                                                                                  tGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTTCAGTATTCCGATTAAGATAAAATCGCGCTGGGCCACGCACA‑ATAAATTGATTTTGCAGAGAAATATGCCAATAAAGAGGCGGTTATTGACAATGTTTTTGTTACTAACCAGGCGATCTGCAAAAAATACGAAGGGCGTGAGGATCCGTTC‑TGAAAGGCGTTGTTTGATGAAAATCGTTTAGATCTGGTGATCCCTG‑TGGGGAACAACACggcgg                                      >  2:529344/1‑247 (MQ=255)
                                                                                                                                                                                                   gTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCC‑GGCTGGAGAACAACACGCCgtg                                     <  1:527331/251‑1 (MQ=255)
                                                                                                                                                                                                    tCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCGATCTTGTGATCC‑TGCTGGAGAACAACACGCCgtgg                                    >  1:183137/1‑251 (MQ=255)
                                                                                                                                                                                                                            agaTCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTc            >  1:114686/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gCAGTATTCCGATTACGATAAAATCGCGCTTGTCCACGCACA‑ATACATTGATTTTGCAGTGAAAGATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCgg  <  2:502867/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
TTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCA‑ATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCA‑GGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGG  >  NC_000913/4627921‑4628401

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