Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,628,208 G→A A298A (GCG→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,2080GA96.0% 147.4 / ‑4.3 50A298A (GCG→GCA) nadRDNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (24/24);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCA‑ATACATTGATTTTGC‑GGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCA‑GGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCT  >  NC_000913/4627981‑4628441
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           
ccTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCTAGCTGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTAGTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGc                                                                                                                                                                                                                      <  1:642891/251‑1 (MQ=255)
 cgacggaagggaatcggtgctttGTACGGACCGTGGCGATCCTTGGTGGCGAGTAGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGTAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGAGTTTGC‑AGTGAAATATGCCAATAAAGTGGcg                                                                                                                                                                                                                     <  2:561247/231‑2 (MQ=255)
          tGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGGGATGAGATCGCATTGCAGGATTCCGATTACGATAAAATCGCGCTGGGCCAAGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGttttttgac                                                                                                                                                                                                            >  1:379039/1‑245 (MQ=255)
           gAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTAttgac                                                                                                                                                                                                            >  1:141713/1‑246 (MQ=255)
                 gttctttGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTAttgacactgatt                                                                                                                                                                                                     <  1:210179/247‑12 (MQ=255)
                        gtacGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAAGAAAGTGGCGTTTATTGACACTGATTTTGTca                                                                                                                                                                                               <  2:90310/247‑1 (MQ=255)
                                 gTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGc                                                                                                                                                                                     >  1:157008/1‑251 (MQ=255)
                                    gCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTTGGCCACGC‑CACATACATTGTTTTTGC‑AGGGAAAAATGCCCATAAAGTGGGGTTTATTGACACGGATTTTGTCACTACCCACGCGtt                                                                                                                                                                                  >  2:161335/1‑251 (MQ=255)
                                     cGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGGTTATTGACACTGATTTTGTCACTACCCAGGCggtt                                                                                                                                                                                 >  1:140859/1‑248 (MQ=255)
                                      gATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGAAAATGATTTTGTCACTACCCAGGCGTTCt                                                                                                                                                                                >  1:451118/1‑251 (MQ=255)
                                                  cgAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTGGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGGCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAAGCGCGCTGGGCCACGCACA‑ATACATTGATTTTGT‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTTCCCAGGCGTTCTGCAAAAAGCACg                                                                                                                                                                    <  2:630423/251‑1 (MQ=255)
                                                    aggcGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATACGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGaa                                                                                                                                                                  <  1:247766/248‑1 (MQ=255)
                                                         aGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGGGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCg                                                                                                                                                             >  1:659709/1‑251 (MQ=255)
                                                                          ggTAAACAAACTTGCCAATATTTTCAATACCCCCAGTGCGGGGGAATATGGTCGCGATTATGTCTTGTCACACCGCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCTCGCACA‑ATACATTGATTTTGC‑AGTGAAAGATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTATCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGt                                                                                                                                              <  2:548739/249‑1 (MQ=255)
                                                                            tAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTTCGAAGGGCGTGAGCATCCGTTCGtac                                                                                                                                            >  1:513682/1‑247 (MQ=255)
                                                                             aaacaaacTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑g                                                                                                                                         <  1:355115/250‑1 (MQ=255)
                                                                               acaaacTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGTTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTTCCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑Ggcg                                                                                                                                      <  2:451118/251‑1 (MQ=255)
                                                                               acaaacTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGGAGGGCGTGAGCATCCGTTCGTACA‑Ggcg                                                                                                                                      >  1:361816/1‑251 (MQ=255)
                                                                                  aacTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCA‑ATACATTGATTTTGC‑GGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCC‑GGCGCtg                                                                                                                                   >  1:308665/1‑251 (MQ=255)
                                                                                      tGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTtgatt                                                                                                                                <  1:371793/250‑1 (MQ=255)
                                                                                      tGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTtgat                                                                                                                                 <  1:336285/249‑1 (MQ=255)
                                                                                       gtcaataatctcaattccaccaGTGCGTGTGAATATGTTCGCGATTATGTCTTTGCACACCTCGGCGGTGATGAGATCGCATTGCAGTAGTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTtgattg                                                                                                                               <  1:595222/234‑1 (MQ=255)
                                                                                                 ccaataccaccAGTGTGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCTGCGGTGATGAGATCGCATTGCAGTATTTCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATTCCGt                                                                                                                    <  2:512062/246‑1 (MQ=255)
                                                                                                    ataccaccAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATAAATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAAAACCGTTTc                                                                                                                 >  1:506680/3‑251 (MQ=255)
                                                                                                     aaccaccAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACTATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCg                                                                                                                <  2:626983/250‑1 (MQ=255)
                                                                                                                cGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGGGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTTACACTGATTTTGTCACTAACCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCGGTCTGGGGATc                                                                                                     >  1:479716/1‑251 (MQ=255)
                                                                                                                   gggAATATGGGCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATTCCGTTTCGATCTGGTGATCc‑t                                                                                                  <  1:671269/250‑1 (MQ=255)
                                                                                                                     gAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATcc                                                                                                                                                    >  1:534478/1‑200 (MQ=255)
                                                                                                                     gAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATcc                                                                                                                                                    <  2:534478/200‑1 (MQ=255)
                                                                                                                        tatGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGGATACCGTTTCg                                                                                                                >  1:334138/1‑232 (MQ=255)
                                                                                                                        tatGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCg                                                                                                                <  2:334138/232‑1 (MQ=255)
                                                                                                                        tatGGTCGCGATTATGTCTTTTCACACCTCGGAGGTGATGAGATCACATTGCCGTATTCCGATTACGATAAAATCGCGCTGGGCCACGAACA‑AAACATTGGTTTTGC‑AGGGAAATATGCCAATAAAGTGGCGATTATTGACACTGATTTTGTCACTAACCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAGCCTCCGTTCGTACA‑GGCGTTGATTGATTAATACCGTTTAGATCTGGTGATAC‑TGCTggg                                                                                            >  1:353187/1‑250 (MQ=255)
                                                                                                                                      tGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTTCAGTATTCCGATTAAGATAAAATCGCGCTGGGCCACGCACA‑ATAAATTGATTTTGC‑AGAGAAATATGCCAATAAAGAGGCGGTTATTGACAATGTTTTTGTTACTAACCAGGCGATCTGCAAAAAATACGAAGGGCGTGAGGATCCGTTC‑TGAAAGGCGTTGTTTGATGAAAATCGTTTAGATCTGGTGATCCCTG‑TGGGGAACAACACggcgg                                                                              >  2:529344/1‑247 (MQ=255)
                                                                                                                                       gTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCC‑GGCTGGAGAACAACACGCCgtg                                                                             <  1:527331/251‑1 (MQ=255)
                                                                                                                                        tCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCGATCTTGTGATCC‑TGCTGGAGAACAACACGCCgtgg                                                                            >  1:183137/1‑251 (MQ=255)
                                                                                                                                                                agaTCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTc                                                    >  1:114686/1‑251 (MQ=255)
                                                                                                                                                                          gCAGTATTCCGATTACGATAAAATCGCGCTTGTCCACGCACA‑ATACATTGATTTTGC‑AGTGAAAGATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCgg                                          <  2:502867/251‑1 (MQ=255)
                                                                                                                                                                                        ccGATAAAATCGCGCTGGGCATCGCAAA‑ATCCATCGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTAGTGACACTGAGTTTGTCACTACCCAGGCGGTCTGCAAAAAGTACGATGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCGAGCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTGACGCAGCCTCGGCAGTTCGGTGGATCGCAAagag                            <  2:19668/250‑1 (MQ=255)
                                                                                                                                                                                          gATAAAAACGGGCTGGGCCACGTACA‑ATACATTGATTTTGCAGGT‑AAATATGCCAATAAAGTGGCGATTATGGACACTGAGTTTGTTACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGAATGCTTTTGTACA‑GGCGTTGATTGATGAATACCGTTTCGGTCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGAAAAGAGtt                          <  2:475808/251‑1 (MQ=255)
                                                                                                                                                                                            aaaaaTCGCGCAGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCTTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTcc                        <  2:588537/250‑1 (MQ=255)
                                                                                                                                                                                              aaaTCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTca                        >  2:119534/1‑248 (MQ=255)
                                                                                                                                                                                                aTCGCGCTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGa                     >  1:150004/1‑250 (MQ=255)
                                                                                                                                                                                                   gcgcTGGGGCACGCTCA‑ATACATTGATTTGGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTTCCCAGGCGTTCTGCAAAAAGTACGAATGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTGCGATCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAAcgt                 <  2:148567/251‑3 (MQ=255)
                                                                                                                                                                                                      cTGGGCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGACGAATACCGTTTCGATCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCt              >  1:264780/1‑251 (MQ=255)
                                                                                                                                                                                                        gggCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTTGCGTTTATTGACACTGATTTTGTCCCTACCCCGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCtg             <  2:21863/250‑1 (MQ=255)
                                                                                                                                                                                                         ggCCACGCACA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCtggt           >  1:418302/1‑251 (MQ=255)
                                                                                                                                                                                                                cacA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCAGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATg    >  1:44696/3‑251 (MQ=255)
                                                                                                                                                                                                                 acA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGc   >  1:613283/2‑251 (MQ=255)
                                                                                                                                                                                                                 acA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTAGATCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGc   >  1:667978/2‑251 (MQ=255)
                                                                                                                                                                                                                  cA‑ATACATTGATTTTGC‑AGTGAAATATGCCAATAAAGTGTCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACA‑GGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCt  <  1:536773/251‑1 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           
CCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCA‑ATACATTGATTTTGC‑GGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCA‑GGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCC‑TGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCT  >  NC_000913/4627981‑4628441

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: ‑