Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,628,782 T→A intergenic (+235/+73) nadR → / ← ettA DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR/energy‑dependent translational throttle protein EttA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,628,7820TA95.7% 140.9 / ‑2.9 46intergenic (+235/+73)nadR/ettADNA‑binding transcriptional repressor/NMN adenylyltransferase NadR/energy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base A (20/24);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGA  >  NC_000913/4628544‑4629015
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         
aTAACCGCGATGAAACGGCACAAAGGCGAGGTTCAAAAAAATTTTTTTGATTGGTGGAAACCCTCTTGGCAAGAGGCCGTCGAGATAAAACAAAGTTTTAAAAAACCTGCACGGGTGCCCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCGGGAAAAGGTGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTAtt                                                                                                                                                                                                                              <  1:93043/251‑1 (MQ=255)
  aaCCGCGATGAAACGGCTCAAAGGCGAGTTATAAAATAATTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTg                                                                                                                                                                                                                             <  1:383935/250‑1 (MQ=255)
                  cTCAAAGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAAtttcttggtttgt                                                                                                                                                                                                             >  1:3273/1‑239 (MQ=255)
                                              gttGATGAGTTGAAATCCTCTCGGGAAGATGCTGTCAAGATTAAACAAGGTTTGACAAACCCTGCACGGGTGACCCGCTACGAAATTGCGGATCCTAAGGCTGTTGTAGGCCGGAAAAGTCGTTTACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                                                                                                  <  1:166798/249‑1 (MQ=255)
                                               tgtatgaGTTGAAACCCTATCTGGAAGAGGTTGTCGAGATTAAACAAGTTTTGAAAAACCCTGCACGGGTGACCCGCCACGAAAGTGCGGATATAAGGGCTTTTGTAGGCCGGAAAAGGCGTTCACTCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                                                                                                 <  1:196238/247‑1 (MQ=255)
                                                    gaGTTGAAAACTTCTCTGGAAGAGGCGGTCGATAGTAAACAAGGTTAGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGAAACTAAGTCTGTGGTAGGCGGGAAAAGTCGTCCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGc                                                                                                                                                                                                         <  1:294763/221‑1 (MQ=255)
                                                      gTTGAAACCCTCTCTGGAAGAGGCGGTGAAGATTAACCAAGGTTTAAAAAAACCTGCCCGGGTGCCCCGCTACGAAATTCCGGATCTTAGGGCTGTTTTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTTTCAGGCCTGGAGATTCACAACGTATTGAAATTGCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAttt                                                                                                                                                                         <  1:392895/251‑1 (MQ=255)
                                                                           agGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGATCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATAc                                                                                                                                                     >  2:407893/1‑251 (MQ=255)
                                                                             gtttacGAGATTAAAAAAGGTGTGAAAAAACCTGACCGGGGCCCCCGGTACGAAATTGCGGATACAAAGGCTGGGGTAGGCGGGAAAAGGCGGTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGGGGCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTCCACGCCTCATCCGGCATTTTACGCATTACTTCGCAATACGc                                                                                                                                                   <  1:175570/246‑1 (MQ=255)
                                                                                   gagaTTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTAc                                                                                                                                             >  1:347729/1‑251 (MQ=255)
                                                                                                        aaaCCCGCACGGGTGACCCCCTAGGAAATTGCGGATCCTAAGGCTTTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAgt                                                                                                                        <  1:193250/251‑2 (MQ=255)
                                                                                                                     tGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCCATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGcccc                                                                                                           >  1:332707/1‑250 (MQ=255)
                                                                                                                     gtaCCCGCTACGAAATGGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCa                                                                                                           <  1:87030/249‑1 (MQ=255)
                                                                                                                            ctaCGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAAtttcttggtttg                                                                                                                                                                                                              >  1:445579/3‑135 (MQ=255)
                                                                                                                            ctaCGAAATTGCGGATACTAAGGCTGTGGCAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAAtttcttggtttg                                                                                                                                                                                                              <  2:445579/143‑11 (MQ=255)
                                                                                                                               acaaaattgcggagtataagcgttttgtAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAgcgtgcgt                                                                                                   <  2:3273/226‑1 (MQ=255)
                                                                                                                                          atacTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCtt                                                                                      >  1:447652/4‑251 (MQ=255)
                                                                                                                                             cTAAGGCTGTTGTAGGCCGGACAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCTCATCCGGCATTTTACGCCTTACTTCACAATACGCTTGTACTTGATAAGCTTCGGCTCCAGTGCATCTGCCCACCGACTGCGTTTCTTGTACTCTTCGt                                                                                   >  2:426072/1‑251 (MQ=255)
                                                                                                                                                     tgttgtAGGCCGAAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTcag                                                                             <  1:414341/250‑3 (MQ=255)
                                                                                                                                                              ccGGAAAAGGCGTTCACTCCGTATCCGGCAGTCGGTGCTCAATGCTTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTAcc                                                                   <  2:8142/251‑1 (MQ=255)
                                                                                                                                                                 aaaaaggggtgtacgccgttcccGGCATTCGGTGCTCAATGCCTGATGCGTCGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTGCCTTc                                                                <  2:25622/230‑1 (MQ=255)
                                                                                                                                                                   aagggctttCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTGCTCAGTAAAGTTACCTTCga                                                              <  1:63140/244‑1 (MQ=255)
                                                                                                                                                                   aaaggcgtTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCg                                                               >  2:3333/7‑250 (MQ=255)
                                                                                                                                                                   aaaggcgtTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCg                                                               >  1:224243/7‑250 (MQ=255)
                                                                                                                                                                                 cgtaTCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACttt                                                 >  2:54351/4‑250 (MQ=255)
                                                                                                                                                                                   taTCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCt                                                                                               <  2:368743/201‑1 (MQ=255)
                                                                                                                                                                                   taTCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCt                                                                                               >  1:368743/2‑202 (MQ=255)
                                                                                                                                                                                      ccGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTc                                           >  1:104500/1‑251 (MQ=255)
                                                                                                                                                                                          cATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATcc                                       >  1:184969/1‑251 (MQ=255)
                                                                                                                                                                                           aTTCGGTGCTCAATGCCTTATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCt                                      <  1:3333/251‑1 (MQ=255)
                                                                                                                                                                                           aTTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAg                                                                                                                                                                                                           <  2:237326/86‑1 (MQ=21)
                                                                                                                                                                                           aTTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAg                                                                                                                                                                                                           >  1:237326/1‑86 (MQ=21)
                                                                                                                                                                                              cGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGt                                   >  2:79724/1‑251 (MQ=255)
                                                                                                                                                                                                     cAATGCCTGTTTCGACGCTTGCCGCGTCTTATCAGTCCTACAGTTTTACAACGTATTGAAATTGCTGGTGTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACGTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATcc                              <  2:474216/249‑1 (MQ=255)
                                                                                                                                                                                                        tGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACttt                                                 <  2:393366/227‑1 (MQ=255)
                                                                                                                                                                                                        tGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACttt                                                 >  1:393366/1‑227 (MQ=255)
                                                                                                                                                                                                                gCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTg                  <  1:54351/250‑1 (MQ=255)
                                                                                                                                                                                                                  gACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACATCGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCg               <  1:177602/251‑1 (MQ=255)
                                                                                                                                                                                                                    cGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGAt             >  1:113038/1‑251 (MQ=255)
                                                                                                                                                                                                                     gCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATa            <  1:383986/251‑1 (MQ=255)
                                                                                                                                                                                                                      cTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATAc           >  1:391172/1‑251 (MQ=255)
                                                                                                                                                                                                                       ttACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACg          >  2:365662/1‑251 (MQ=255)
                                                                                                                                                                                                                        tACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACg          <  2:104500/250‑1 (MQ=255)
                                                                                                                                                                                                                              gTCTTAGCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAg    <  1:482679/250‑1 (MQ=255)
                                                                                                                                                                                                                                cTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCCTTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGa  >  2:424237/1‑250 (MQ=255)
                                                                                                                                                                                                                                  tATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGa  <  2:272887/248‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         
ATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGA  >  NC_000913/4628544‑4629015

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