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.2% 156.4 / NA 53A298A (GCG→GCAnadRDNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (20/30);  minor base C (2/0);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAA  >  NC_000913/4627967‑4628453
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
aCTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTggaaaat                                                                                                                                                                                                                                              >  1:83890/1‑245 (MQ=255)
     gAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGGGAGAGCATTGCCAAATTTCGATTACGATAAAAACGCGCTTGGCCACGCCCCATACCTTGATTTTGCAGTGGAataatcca                                                                                                                                                                                                                                          >  2:108040/1‑245 (MQ=255)
          tatccctagggaggTGAAGCCTTTCTTTGTACGTACCGTGGCGTCCCTTGGGGGGGATTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCCCCAGTGCGTGGGAATAGGGTCGCGATTATGTCTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAg                                                                                                                                                                                                                                    <  2:165132/237‑1 (MQ=255)
           atTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCTATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAAGATGCCAATAAAg                                                                                                                                                                                                                                    <  1:352336/250‑1 (MQ=255)
           atTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGGGGTGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTTATTTTGGCGTGAAATATGCCAATAAAgg                                                                                                                                                                                                                                   >  2:66016/1‑250 (MQ=255)
             tCCTACGGAAGTGAAGCCGGTCTTGTTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACTCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTg                                                                                                                                                                                                                                  <  2:466039/250‑1 (MQ=255)
                         gAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATACACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGACGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGAt                                                                                                                                                                                                                      <  1:58279/250‑1 (MQ=255)
                            gCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCCCAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTAttgacactg                                                                                                                                                                                                                  >  1:206426/1‑243 (MQ=255)
                              cGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACCTTGATTTTGCCGTGAAATATGCCAATAAAGTGGCGTTTAtttacactgat                                                                                                                                                                                                                >  1:85353/1‑241 (MQ=255)
                                            aCCGTGGCGATCCTTGGTGGCGTGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATTCCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTGTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACAGTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATGGACACTGATTTTGTcactaccca                                                                                                                                                                                                  <  2:279857/251‑7 (MQ=255)
                                             ccGTGGCGATCCTTGGTTGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAATCTTGCCAATATATTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTcactacccag                                                                                                                                                                                                 <  1:14540/251‑8 (MQ=255)
                                                 ggCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCTCTACCCAGGCGt                                                                                                                                                                                             <  1:116136/251‑1 (MQ=255)
                                                     aTCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTTAAACTGATTTTTTCACTACCCAGGCGTTCTg                                                                                                                                                                                         >  1:243321/1‑251 (MQ=255)
                                                            gtggcgAGTCGAGCGGTAAATCCACCCTGTTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAg                                                                                                                                                                                  <  1:122377/249‑1 (MQ=255)
                                                                  aGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTTTTGACACTGATTTTGTCACTACCCCGGCGTTCTGCAAAAAGTACGa                                                                                                                                                                             >  1:314166/1‑250 (MQ=255)
                                                                     cGAGCGGTAAATCCACCCTGGGAAACAAACTTGCCAATATCTTCAATACCTCCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAggg                                                                                                                                                                         <  1:45040/251‑1 (MQ=255)
                                                                           gTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCgtg                                                                                                                                                                     >  1:187666/1‑247 (MQ=255)
                                                                                                                 aatacccccaGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGAGTTTGTCACTTCCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGttt                                                                                                                             <  1:208931/244‑1 (MQ=255)
                                                                                                                  ataccaccAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTc                                                                                                                            <  2:375063/249‑1 (MQ=255)
                                                                                                                    accaccAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCg                                                                                                                           <  2:343130/250‑1 (MQ=255)
                                                                                                                          aGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTAttgac                                                                                                                                                                                                                      <  2:414581/153‑5 (MQ=255)
                                                                                                                          aGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTAttgac                                                                                                                                                                                                                      >  1:414581/1‑149 (MQ=255)
                                                                                                                              cgggggAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAAGCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGa                                                                                                                  <  2:58264/246‑1 (MQ=255)
                                                                                                                              cGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCCCAATACATTGATTTTGCAGTGAAATATGCCCATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTTCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGa                                                                                                                  >  1:58264/1‑249 (MQ=255)
                                                                                                                              cGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATc                                                                                                                <  1:73849/251‑1 (MQ=255)
                                                                                                                                 gggAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCct                                                                                                              <  1:66016/250‑1 (MQ=255)
                                                                                                                                                aTTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACg                                                                                              <  1:322976/251‑1 (MQ=255)
                                                                                                                                                     gTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATCCCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCgt                                                                                          <  1:108040/250‑1 (MQ=255)
                                                                                                                                                      tCTTTTCACACCTCGGCGGTGATGAGATCGCATTTCAGTATTTCGATTACGAAAAAATCGCGCTGGGCCACGCACAATACCTTTATTTTGCAGTGAAATATGCCAATAAAGGGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTATGCAAAAAAGACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGGATACCGGTTCGGTCTGGTGGTCCTGCTGGGGAAAAACACGCCgtgg                                                                                        >  2:153365/1‑251 (MQ=255)
                                                                                                                                                            cacaCCACGGCGGTGATGAGTTCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTCCCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGAATGAGGAATTCCGTGTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCg                                                                                  <  1:46028/251‑1 (MQ=255)
                                                                                                                                                              cacCGCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGa                                                                                <  1:121861/251‑1 (MQ=255)
                                                                                                                                                               cccTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGAt                                                                               <  2:85353/250‑1 (MQ=255)
                                                                                                                                                                                aTCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACCTTGATTTTGCAGTGAAATATGCCAATCAAGTGGCGTTTATTGACACTGATTTTGTCACTAACCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATAACGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCgg                                                              >  1:227141/1‑251 (MQ=255)
                                                                                                                                                                                  cGCATTGCAGTATTCCGATTACGATCAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGc                                                             >  1:269327/1‑250 (MQ=255)
                                                                                                                                                                                         cagtactccgataacgctaacctCGCGCTGGGCCACTCACAATACATTTATTTTGCAGAGAAATATGCCAATAAGGAGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGCAGGGCGTGAGCATCCTTTCGCACAGGCGTTGATTGATGAATACAGTTTCGATCTGGTGATCCTGCTGGAGAACAAAACGCCTTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTccgg                                                     >  2:298891/23‑248 (MQ=255)
                                                                                                                                                                                         cAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGCTGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCggg                                                     >  2:322169/1‑250 (MQ=255)
                                                                                                                                                                                                ccgattACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGAAACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCg                                               <  2:366605/245‑1 (MQ=255)
                                                                                                                                                                                                   attACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCaaa                                           <  2:243321/249‑1 (MQ=255)
                                                                                                                                                                                                              aTCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGaa                                <  1:124694/251‑1 (MQ=255)
                                                                                                                                                                                                                cgcgCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGGGGATCGCAAAGAGTTCCAGAACt                              >  2:12951/1‑251 (MQ=255)
                                                                                                                                                                                                                     tGGGCCACGCACAGTACATTGATTTTGGAGTGAAATATGCCATTAAAGAGGCGTTTATTAACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGCAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGGTGAATACCGTTTCGATCTGGTGATCCTGCTGGTGAACAACACGCCGTGGGTGGCGGTTGGTTTACGCAGCCTCGGCAGTTCGGTGGAACGCAAAGCGTTCCAGAACTTGCtg                         >  2:63140/1‑251 (MQ=255)
                                                                                                                                                                                                                       ggCCACGCACAATGCATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGTCACTGATTTTGTCACTCCCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGTGGATGGTTTACGCAGCCTCGGCAGTTCGGGGGATCGCAAAGAGTTCCAGAACTTGCtggt                       <  1:95368/251‑1 (MQ=255)
                                                                                                                                                                                                                       ggCCACGCACAATACATGGAGTTTGCAGTGAAATATGCCAATAAAGTGGCTTTTATTGACACTGATTTTGTCACTACCCAGTCGTTCTGCAAAAAGTACGAAGGGCTGGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCTGCCTCGGCAGTTCTGTGGATCGCAAAGAGGTCCAGAACTTGCtggt                       <  2:189985/251‑1 (MQ=255)
                                                                                                                                                                                                                           ccgcacAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGgag                   <  1:480659/246‑1 (MQ=255)
                                                                                                                                                                                                                             gcacAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCtggtgg                     >  2:422154/4‑247 (MQ=255)
                                                                                                                                                                                                                             gcacAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACAATGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCACGGCCGTTCGGTGGGACTCAACGAGTTACAGAAAGTGCtggtgg                     >  2:253817/4‑247 (MQ=255)
                                                                                                                                                                                                                              cacAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTGTATTGACACTGATTTTGTCACGACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGGTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATg                <  2:269327/249‑1 (MQ=255)
                                                                                                                                                                                                                                cAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCt              >  1:205444/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aaTACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCtg             >  1:403702/1‑250 (MQ=255)
                                                                                                                                                                                                                                    aCATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCTTTCTGCAAAAAGTACGAAGTGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGaa          <  1:339407/251‑1 (MQ=255)
                                                                                                                                                                                                                                        tGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCGAAGAGTTCCAGATCTTGCTGGTGGAGATGCtggacg         >  2:196238/1‑243 (MQ=255)
                                                                                                                                                                                                                                        tGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTGCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAagag       <  2:128399/250‑1 (MQ=255)
                                                                                                                                                                                                                                             ttGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACaa  <  1:427413/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
ACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAA  >  NC_000913/4627967‑4628453

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