Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,628,106 C→T G264G (GGC→GGT 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,1060CT97.2% 209.8 / NA 71G264G (GGC→GGTnadRDNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (37/32);  minor base G (0/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGG  >  NC_000913/4627862‑4628350
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
cgcCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGAt                                                                                                                                                                                                                                                <  1:168333/251‑1 (MQ=255)
  ccACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGAtt                                                                                                                                                                                                                                               >  1:161208/1‑250 (MQ=255)
   cacaGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTAt                                                                                                                                                                                                                                             <  1:619238/251‑1 (MQ=255)
      aGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCGTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTc                                                                                                                                                                                                                                          <  2:551242/251‑1 (MQ=255)
         atatGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGAGCCGCGAAAACCCATTCCGCTACGGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCGCCCTGGTAAACAAACTTGCCAATATCTTCAAGACCGCCAGTGCGTGGGAATATGGGCGCGATTATGTCtt                                                                                                                                                                                                                                        <  2:571692/250‑1 (MQ=255)
         atatGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTAGGAGTATCAGCGGTGCTCAGATCCGCGAAATCCCGTTCCGCTACGGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCTCCAGTGCGTGGGAATATGGTCGCGATTATGTCtt                                                                                                                                                                                                                                        <  2:555860/250‑1 (MQ=255)
           atGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTc                                                                                                                                                                                                                                     >  1:601712/1‑251 (MQ=255)
           atGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTc                                                                                                                                                                                                                                     >  1:26752/1‑251 (MQ=255)
              gAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGATTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGTGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTcaca                                                                                                                                                                                                                                  <  1:579423/251‑1 (MQ=255)
                      ggggATCGATACGGTGCTGGTCGATCCGAAACGTTCCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTc                                                                                                                                                                                                                              <  1:623519/247‑1 (MQ=255)
                        ggATCGATACGGTGCTGGTCGATCCGAAACGTTCCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACGGGGAATTTATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGAGCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGTTCGCGATTATGTCTTTTCACTCCTCGGCGGt                                                                                                                                                                                                                        <  2:397725/251‑1 (MQ=255)
                                     gCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCa                                                                                                                                                                                                            <  1:41476/250‑1 (MQ=255)
                                        ggTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGc                                                                                                                                                                                                        <  1:520433/251‑1 (MQ=255)
                                                     cGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTccgatt                                                                                                                                                                                            >  1:85747/1‑245 (MQ=255)
                                                       tACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGCGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTTTTccgatta                                                                                                                                                                                           >  1:538505/1‑243 (MQ=255)
                                                       tACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACg                                                                                                                                                                                         >  1:136359/1‑251 (MQ=255)
                                                            ttATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGGGATCGCATTGCAGTATTCCGATTACGATaa                                                                                                                                                                                     >  1:555636/1‑250 (MQ=255)
                                                                   tATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTTCGATTACGATAAAATCgcgc                                                                                                                                                                             >  2:475905/1‑251 (MQ=255)
                                                                      cAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATcaaatcat                                                                                                                                                                          <  1:328719/251‑8 (MQ=255)
                                                                           gTGCGCAGATCCGCGACAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTAACGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGACGTATTCCGTTTCAGATAAAATCGCGCTGGGCCAc                                                                                                                                                                     >  2:53901/1‑251 (MQ=255)
                                                                               gcAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTTCGATTACGATAAAATTGCGCTGGGCCAcccac                                                                                                                                                                 >  2:96748/1‑247 (MQ=255)
                                                                                 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)
                                                                                                  ttCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACAcc                                                                                                                                                                                                                                >  1:468340/1‑169 (MQ=255)
                                                                                                  ttCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACAcc                                                                                                                                                                                                                                <  2:468340/169‑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)
                                                                                                                                                                                           cACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCa                                                                                                                                                                                                            <  2:319499/100‑1 (MQ=255)
                                                                                                                                                                                           cACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCa                                                                                                                                                                                                            >  1:319499/1‑100 (MQ=255)
                                                                                                                                                                                                 ggTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGt                                                 <  2:69856/249‑1 (MQ=255)
                                                                                                                                                                                                       caaacTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCAtt                                                                                                                                                                                                          >  1:84187/1‑90 (MQ=255)
                                                                                                                                                                                                       caaacTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCAtt                                                                                                                                                                                                          <  2:84187/90‑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)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
CGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGG  >  NC_000913/4627862‑4628350

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