Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,362,045 T→C intergenic (‑111/+506) cadC ← / ← pheU DNA‑binding transcriptional activator CadC/tRNA‑Phe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,362,0450TC95.6% 160.8 / ‑2.8 45intergenic (‑111/+506)cadC/pheUDNA‑binding transcriptional activator CadC/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (23/20);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACAAGCCATTCGCCAACGCGAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTTTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTC  >  NC_000913/4361803‑4362284
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
aCTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACAAGCCATTCGCCAACGCGAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGc                                                                                                                                                                                                                                         >  2:312288/1‑251 (MQ=255)
         aGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACAAGCCATTCGCCAACGCGAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAActt                                                                                                                                                                                                                                >  1:104253/1‑251 (MQ=255)
            aaaGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACAAGCCATTCGCCAACGCGAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAActtctt                                                                                                                                                                                                                             <  1:336848/251‑1 (MQ=255)
             aaGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACAAGCCATTCGCCAACGCGAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTg                                                                                                                                                                                                                            <  1:281543/251‑1 (MQ=255)
               gAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACAAGCCATTCGCCAACGCGAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAAGCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGAt                                                                                                                                                                                                                          >  2:67465/1‑251 (MQ=255)
               gAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACAAGCCATTCGCCAACGCGAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCACAAACTTCTTGAt                                                                                                                                                                                                                          <  1:389283/251‑1 (MQ=255)
                        aagaTCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACAAGCCATTCGCCAACGCGAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTa                                                                                                                                                                                                                 <  2:320105/251‑1 (MQ=255)
                        aagaTCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACAAGCCATTCGCCAACGCGAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTa                                                                                                                                                                                                                 <  1:232657/251‑1 (MQ=255)
                                  aaTCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACAAGCCATTCGCCAACGCGAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAgg                                                                                                                                                                                                       >  2:288041/1‑251 (MQ=255)
                                       cGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACAAGCCATTCGCCAACGCGAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGAtttn                                                                                                                                                                                                  <  2:246944/251‑2 (MQ=255)
                                        ggCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACAAGCCATTCGCCAACGCGAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTa                                                                                                                                                                                                 <  1:126119/251‑1 (MQ=255)
                                                 gTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACAAGCCATTCGCCAACGCGAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAAcc                                                                                                                                                                                        <  2:234207/251‑1 (MQ=255)
                                                                  gCGGCTAATTTGGTTTATGGACGGAGTAACAAGCCATTCGCCAACGCGAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCtt                                                                                                                                                                        <  2:85700/250‑1 (MQ=255)
                                                                   cGGCTAATTTGGTTTATGGACGGAGTAACAAGCCATTCGCCAACGCGAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTAc                                                                                                                                                                      >  1:277492/1‑251 (MQ=255)
                                                                                      aCGGAGTAACAAGCCATTCGCCAACGCGAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCtt                                                                                                                                                   >  2:135275/1‑251 (MQ=255)
                                                                                       cGGAGTAAGAAGCCATTCGCCAACGCGAACTACAGGTTGTTGCAGAATAGAAACGCAGTCGAAAAGGGAATAATGCAATGATAATTGCCACAACCTATTTTTTCCAGCTATTGATGGGTTTATTTAATTATTATTGGTGCATGCAAGACGTTATTTTTCCCTGCCAGAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTaacgtcactattcagcatctgt                                                                                                                                                  <  2:557970/251‑22 (MQ=255)
                                                                                                    cATTCGCCAACGCGAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAAt                                                                                                                                     <  2:492894/251‑1 (MQ=255)
                                                                                                         gCCAACGCGAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCa                                                                                                                                <  1:288041/251‑1 (MQ=255)
                                                                                                          ccAACGCGAACTAAAGGTTGTTGAATAATAGAAACTTATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAg                                                                                                                               <  2:290318/251‑1 (MQ=255)
                                                                                                               gcgAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAaat                                                                                                                          >  1:471359/1‑251 (MQ=255)
                                                                                                                                gCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGt                                                                                                         <  1:372948/251‑1 (MQ=255)
                                                                                                                                    aataGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCgg                                                                                                     >  2:543589/1‑251 (MQ=255)
                                                                                                                                            cTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACtt                                                                                              >  1:517582/1‑250 (MQ=255)
                                                                                                                                                ttCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTAt                                                                                         >  1:57204/1‑251 (MQ=255)
                                                                                                                                                  cGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCTTGGAAAATTGTCCGGAGCACTTACTAttt                                                                                       >  2:541118/1‑251 (MQ=255)
                                                                                                                                                   gAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTAtttt                                                                                      >  1:174495/1‑251 (MQ=255)
                                                                                                                                                          gAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGAt                                                                               >  1:510684/1‑251 (MQ=255)
                                                                                                                                                          gAATGAGGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGAt                                                                               <  2:463276/251‑1 (MQ=255)
                                                                                                                                                           aaTGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGAtt                                                                              <  2:531978/251‑1 (MQ=255)
                                                                                                                                                               atgCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTa                                                                          <  2:395725/251‑1 (MQ=255)
                                                                                                                                                                              ccACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAGc                                                            >  2:223328/1‑250 (MQ=255)
                                                                                                                                                                                  aaCCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAGCAAGCg                                                       <  2:378215/251‑1 (MQ=255)
                                                                                                                                                                                   aCCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAGCAAGCg                                                       >  1:22327/1‑250 (MQ=255)
                                                                                                                                                                                       aTTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAGCAAGCGAACTg                                                  >  2:512696/1‑251 (MQ=255)
                                                                                                                                                                                                aTCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAGCAAGCGAACTGATACATTtc                                         >  2:342612/1‑251 (MQ=255)
                                                                                                                                                                                                    ataGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCtt                                     >  1:130449/1‑251 (MQ=255)
                                                                                                                                                                                                          gggTTTATTTACATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACCTAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTtct                               >  1:218364/1‑251 (MQ=255)
                                                                                                                                                                                                                atttaCATATTATTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCa                         <  1:135275/251‑1 (MQ=255)
                                                                                                                                                                                                                       tattatTGGTGAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTCGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAAtc                  >  1:87488/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gAATGCAAGACGTTATTTTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCg        >  2:529220/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gAATGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTc         <  1:512696/250‑1 (MQ=255)
                                                                                                                                                                                                                                  aaTGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGc                                                                                                  >  2:315147/1‑160 (MQ=255)
                                                                                                                                                                                                                                  aaTGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGc                                                                                                  <  1:315147/160‑1 (MQ=255)
                                                                                                                                                                                                                                  aaTGCAAGACGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGc       >  2:497986/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aCGTTATTCTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTc  >  1:272321/1‑248 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
ACTGTGTTGAGCAAAGAAAACCAGAAGATCGATTAATCTCGGCTCAAGGGTAAGTTGACGCCCATTGCGGCTAATTTGGTTTATGGACGGAGTAACAAGCCATTCGCCAACGCGAACTACAGGTTGTTGCATAATAGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTTTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTC  >  NC_000913/4361803‑4362284

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