Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,630,022 T→C E167E (GAA→GAG ettA ← energy‑dependent translational throttle protein EttA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,630,0220TC90.2% 94.9 / 2.9 41E167E (GAA→GAGettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base C (14/23);  total (17/24)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.59e-01

GCTCCAGCCAGGAGGAGTAGTTACCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACAACGGTGCCTTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGC  >  NC_000913/4629775‑4630265
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
gCTCCAGCCAGGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCTCAACGCTACGCGACGACGCTca                                                                                                                                                                                                                                                   <  2:299106/251‑1 (MQ=255)
   ccagccagGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTcacca                                                                                                                                                                                                                                                <  1:240826/251‑1 (MQ=255)
     agccagGAGGAGTAGTTGCCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTTGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCaccacc                                                                                                                                                                                                                                              <  1:216035/251‑1 (MQ=255)
                      gccTTCCCACGGAATTCCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGAttt                                                                                                                                                                                                                             <  1:329029/250‑1 (MQ=255)
                                       cTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCCCACCGCTACGAGAAGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGACCCAGTCCGGc                                                                                                                                                                                                            >  1:71182/1‑251 (MQ=255)
                                       attCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGGTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCACTCGTAGACGTGCAGGAAGCGTTCCAGCCAGGCCACGGTTTCGGCATCCAGGTGGTTGGTCGTGTCGTCGAGCAGAAGCATGTCTGGTTTTTCCAGACGCAGGCGGCACACCGCTACGCGACGACGCTCACCACCGAAGCGGTTAGAGATTTTCGCGTCCCAGTCCGGc                                                                                                                                                                                                            >  1:71188/2‑251 (MQ=255)
                                                    aaGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGACGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACCCTACGCGACGACGCTCACCCCCGGAGAGGTTAGCGTTTTTCGCGTCCCCGTCCGGCGACCGTAGCGCAt                                                                                                                                                                                              >  1:339263/1‑251 (MQ=255)
                                                                  aGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTcc                                                                                                                                                                                >  1:298332/1‑251 (MQ=255)
                                                                   gCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCa                                                                                                                                                                               <  1:82496/251‑1 (MQ=255)
                                                                                 cGAGGAAGTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAg                                                                                                                                                                  <  1:150098/250‑1 (MQ=255)
                                                                                        gTAACGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGAc                                                                                                                                                          <  1:188476/251‑1 (MQ=255)
                                                                                          aaCGGTCGTGGGTAATCGCCACCACGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACcg                                                                                                                                                        >  1:384866/1‑251 (MQ=255)
                                                                                                        accgcccccacgttgacctCGAAGTCGGGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGGTGGTCGGTGCGTCGGTCAGCTGCATGTCTGTTTTTGCCAGCAGCAGGCGGCAAAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCTTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTAAGTTCAGGTTATGACCGTCGTGAGCCTGAAt                                                                                                                                          <  2:45584/233‑1 (MQ=255)
                                                                                                                caCGGTGCCCTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCtt                                                                                                                                  <  1:179145/250‑1 (MQ=255)
                                                                                                                caCGGTGCCCTCGAAGTCGGGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCCCAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCtt                                                                                                                                  <  1:168481/250‑1 (MQ=255)
                                                                                                                        cctCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACgg                                                                                                                          <  1:110781/249‑1 (MQ=255)
                                                                                                                            gAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCtt                                                                                                                      >  1:368110/1‑251 (MQ=255)
                                                                                                                              aGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGCTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGt                                                                                                                    >  1:9463/1‑251 (MQ=255)
                                                                                                                                  gTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGAATCCAGGCGGTTGGTCGGTTCGTCGAGCAGCAGAATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGCCGACGCTCACCACCAGAGAGGATAGCGATTTTAGCGTCACAATCCGGCAGAAGTAGCGCATCCGCCGCACGCTCCAGCTGAACGTTCAGCAGATGACCGACGTGAGCCTGACAGACCTCTTCCAAACCGCCTTGTTc                                                                                                                  >  2:370394/1‑249 (MQ=255)
                                                                                                                                           gCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCt                                                                                                       <  1:349988/251‑1 (MQ=255)
                                                                                                                                                  agccagGCCACGGATTCGGCATCCCGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTCTGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCCGCCTGTCaa                                                                                                 >  2:58972/1‑250 (MQ=255)
                                                                                                                                                   gccagGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCaaa                                                                                                >  1:229797/1‑250 (MQ=255)
                                                                                                                                                     cagGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATc                                                                                             >  1:191384/1‑251 (MQ=255)
                                                                                                                                                                      aTCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGCCGTGAGCCTGCATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAACTCGGCATCCGGCTCGGCGt                                                                            >  2:21168/1‑251 (MQ=255)
                                                                                                                                                                            gtggtTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAgcg                                                                      <  1:275728/251‑1 (MQ=255)
                                                                                                                                                                                  ggtCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATaca                                                                <  1:93909/251‑1 (MQ=255)
                                                                                                                                                                                  ggtCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATaca                                                                <  1:83874/251‑1 (MQ=255)
                                                                                                                                                                                  ggtCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATac                                                                 >  1:110205/1‑250 (MQ=255)
                                                                                                                                                                                             cgAGCAGCCGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACACCGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGCCGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAgg                                                     <  2:9463/251‑1 (MQ=255)
                                                                                                                                                                                             cgAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAgg                                                     <  1:237776/251‑1 (MQ=255)
                                                                                                                                                                                              gcgctgcagcaTGTCTGGTTTTTCCAGCAGCAGGCGGCTCAACGTTACGCTCCGACGCTCCCCACCGGAGAGGTTAGCGATTTTCGCGTCCCATGCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCTTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGc                                                    <  2:110205/246‑1 (MQ=255)
                                                                                                                                                                                              gAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGc                                                    <  1:377042/251‑1 (MQ=255)
                                                                                                                                                                                               ggaggcagccTGTGTGGTTTGGCCAGCAGCGGGGGGGAGAACGGTACGCGACGAGGCTCGCCACCGGAGAGGGTAGGGATGTTCGCGTCCCAGTCGGGCAGACGTAGGGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGAGCTCTGCCAGACGGGCTTGTTCAGCGGCCAGCTTGGCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCg                                                   <  2:111651/247‑1 (MQ=255)
                                                                                                                                                                                                      gcaatgcTGGTTTTTCCAGCAGCAGGCGGCACAACCCTACGCGACGGCGCTCACCACCGGAAAGGCTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCACCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCACGTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATTCACTTCATCCAGGCGTTTCAgg                                            <  1:249593/245‑1 (MQ=255)
                                                                                                                                                                                                              ggTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCAGCCGCCGCACGCTCCAGCTGTACGTTCAGATTTTGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAAc                                    <  1:55967/251‑1 (MQ=255)
                                                                                                                                                                                                                      cagcagcagGCGGCACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAgc                                                                       <  1:260804/208‑1 (MQ=255)
                                                                                                                                                                                                                            ccggcggcAAAACGCTACGAGACGACGCTCACCACCGGAGAGGCTAGCGATTTTCGCGTCCCAGTCCGGCAGACGGAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCCGGTTATGACCGTCGAGAGACTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGCAGGACAGCGCATACACTTCATCCAGGCGTTTCAGGACGTTAACCACGTCAGAAAc                        >  2:170772/3‑249 (MQ=255)
                                                                                                                                                                                                                                 ggcACAACGCTACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCAGAAACCGcttct                 <  1:21168/251‑1 (MQ=255)
                                                                                                                                                                                                                                           tACGCGACGACGCTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGGTTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGTCGTTAACCACTTCAGAAACCGCTTCTTCAATGGACt       >  1:267630/1‑251 (MQ=255)
                                                                                                                                                                                                                                                gacgacGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTATGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTcacgc  >  1:90104/1‑251 (MQ=255)
                                                                                                                                                                                                                                                gacgacGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCATTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTcacgc  >  1:25917/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
GCTCCAGCCAGGAGGAGTAGTTACCTTCCCACGGAATACCTTCACCGCGGTCAAGTTCGAGGATCCAGCCTGCAACGTTATCGAGGAAGTAACGGTCGTGGGTAATCGCCACAACGGTGCCTTCGAAGTCGTGCAGGAAGCGTTCCAGCCAGGCCACGGATTCGGCATCCAGGTGGTTGGTCGGTTCGTCGAGCAGCAGCATGTCTGGTTTTTCCAGCAGCAGGCGGCACAACGCTACGCGACGACGTTCACCACCGGAGAGGTTAGCGATTTTCGCGTCCCAGTCCGGCAGACGTAGCGCATCCGCCGCACGCTCCAGCTGTACGTTCAGATTATGACCGTCGTGAGCCTGAATGATCTCTTCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGC  >  NC_000913/4629775‑4630265

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: