Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,425,688 C→A H57N (CAC→AAC)  priB → primosomal replication protein N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,425,6880CA98.1% 160.2 / NA 52H57N (CAC→AAC) priBprimosomal replication protein N
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (26/25);  minor base G (0/1);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CCGTGAGCGTCGCGATGATTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGACACGAAAACCAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTA  >  NC_000913/4425444‑4425933
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
ccGTGAGCGTCGCGATGATTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAGACGaaa                                                                                                                                                                                                                                                 <  1:120067/251‑1 (MQ=255)
 cGTGAGCGTCGCGAAGATTAAGAAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGaaa                                                                                                                                                                                                                                                 >  1:21165/1‑250 (MQ=255)
               tgatTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACAGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGGGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCAATTACTc                                                                                                                                                                                                                                  >  2:181259/1‑251 (MQ=255)
                      gCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGtat                                                                                                                                                                                                                           >  1:291399/1‑251 (MQ=255)
                       cAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGCCCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGTTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGtata                                                                                                                                                                                                                          <  2:364661/251‑1 (MQ=255)
                       aaaaagAAAGCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGtata                                                                                                                                                                                                                          <  2:221629/246‑1 (MQ=255)
                              aaCCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCg                                                                                                                                                                                                                   >  1:157330/1‑251 (MQ=255)
                                    tgatgaTGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCCCAGTATAACGGTCGGCAGTc                                                                                                                                                                                                             >  2:197301/1‑251 (MQ=255)
                                      atgatgCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGc                                                                                                                                                                                                           >  1:330957/1‑251 (MQ=255)
                                         atgCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGTAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATa                                                                                                                                                                                                        <  1:250523/251‑1 (MQ=255)
                                         atgCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGAAGACCAACCGTCTGGTGTTGTCCGGCACCGTGAGCAGGGCTCCCCTTCGAACGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGACGCAt                                                                                                                                                                                                         >  1:24534/1‑250 (MQ=255)
                                              gAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAAcc                                                                                                                                                                                                     <  1:205207/249‑1 (MQ=255)
                                              gAAGCTGGGGATTATGAAGAGTAATTTCTGATGAACAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGGTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACCGTATAACGGTCGGCCGTCGCATAACCGt                                                                                                                                                                                                   >  2:355840/1‑251 (MQ=255)
                                                    ggggATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAgg                                                                                                                                                                                              <  2:261894/250‑1 (MQ=255)
                                                                      tttCTGATGGCCAAACGTCTGGTGTTGTCCGGCACCGTGGGGAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGACTTCCTAACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTGGCGGAAACGTAAACCAGGCCATTAATCACAGTCTAACGGTCGGGCGACGCCTAACCGTTCAGGGGttcatttcatcacaca                                                                                                                                                                            >  2:398912/1‑244 (MQ=255)
                                                                          tgatgaCCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGc                                                                                                                                                                        >  1:276948/1‑250 (MQ=255)
                                                                           gatgaCCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTa                                                                                                                                                                                                                                     >  1:64861/1‑188 (MQ=255)
                                                                           gatgaCCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTGTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAATCGAAAAACAGGCCATTa                                                                                                                                                                                                                                     <  2:64861/188‑1 (MQ=255)
                                                                              gaCCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCaa                                                                                                                                                                      >  2:110746/1‑248 (MQ=255)
                                                                               aCCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAAAGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGa                                                                                                                                                                   >  1:197324/1‑250 (MQ=255)
                                                                                ccAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGCAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAAc                                                                                                                                                                 <  1:61468/251‑1 (MQ=255)
                                                                                ccAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGGCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAAc                                                                                                                                                                 <  1:425525/251‑1 (MQ=255)
                                                                                               tgtCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGGTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAGGAACGGACTGAGCAAAAt                                                                                                                                                   <  2:276948/250‑1 (MQ=255)
                                                                                                    ggcgccGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGtttt                                                                                                                                             <  2:318698/247‑1 (MQ=255)
                                                                                                    ggCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGtttt                                                                                                                                             >  1:318698/1‑251 (MQ=255)
                                                                                                      cACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTTTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGc                                                                                                                                           <  2:263091/251‑1 (MQ=255)
                                                                                                      cACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGc                                                                                                                                           <  1:315775/251‑1 (MQ=255)
                                                                                                                      ccccTTCGAAAGGTCAGTCCATAAGGAATTCCTCACGGCCAGTTCGTGCTGGAGCAGCGTTCTGGGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGGGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGGCGCATAACCGTTCAGGGGTGCATTTCAGGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCAGGCCGAGTAGattg                                                                                                                            <  2:404856/250‑1 (MQ=255)
                                                                                                                      ccccTTCAAAAGGTCAGGCCATCAGGACTTCCCCACTCCCAGTTCGTGCGTGAGCATCGTACTGTGCAGGAGGAAGCCGCTTTTCCCCGGCAGTCGTGGTGTCAAATTCCCGTTATTGTTAGCGGAACCTACAAACAGTCCATTACTCACAGTATAACGGTCGGCAGTCGCATCGCCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACTGACTGAGCAAAATGGTTTTTCATGCCGAGCAGattg                                                                                                                            <  2:105034/250‑1 (MQ=255)
                                                                                                                         cTTCGCAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTGTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGATCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCAGAGCAGattggatt                                                                                                                        >  2:125568/1‑251 (MQ=255)
                                                                                                                               aaGGTCAGTCCAGCAGTAATTCCTCACTGCCCGTTCGTGCTTGAGCATCGGTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGTCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGa                                                                                                                  <  1:292989/251‑1 (MQ=255)
                                                                                                                                       ccccgcAGGAATTCCTTCCTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCTTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACTGACTGAGCAAAATGTTGTTGCATTCCGAGCAGATTGAATTGATAGATTCTGgag                                                                                                          <  1:390896/246‑1 (MQ=255)
                                                                                                                                           tCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACt                                                                                                       >  1:7816/1‑250 (MQ=255)
                                                                                                                                           tCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACt                                                                                                       >  1:244992/1‑250 (MQ=255)
                                                                                                                                               gAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGt                                                                                                                                                                                                                             <  1:63857/128‑1 (MQ=255)
                                                                                                                                               gAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGt                                                                                                                                                                                                                             >  2:63857/1‑128 (MQ=255)
                                                                                                                                                  ttCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAACCAGGCCATTACTCACAGTATAACGATCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCatat                                                                                               >  1:404308/1‑251 (MQ=255)
                                                                                                                                                                     gCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCcgtcg                                                                             <  1:197301/250‑1 (MQ=255)
                                                                                                                                                                       ttGAGAATCGTTCTGTGCAGGTGGAAGCCGGATTTCACCGGCAGGCGTGGCGTCAAATGCCCGATATTGTTAGCGGAACCGCACAACAGGACATTACTAACAGTCTAACGGTCGGCAGTCGCATAAACGTTCAGGGGTTCATTTCATACCACAAGGCCAAGACCGGACTGAGCAAAATGGTTTTGCATGACGAGCAGATTGAATTGATAGATTCTGGAGACACGCCACATGGCACGTTATTTCCGtcgtcg                                                                          >  2:64683/1‑251 (MQ=255)
                                                                                                                                                                                     gtgCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGt                                                             <  1:110746/250‑1 (MQ=255)
                                                                                                                                                                                       gCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTTGTGTCAAATGCCCGTTATTGTTAGCGGAAACTAAAATCAGGCCATTACTCACAGTATAACGTTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGATCAGATTGAATTGATAGATTCTGTAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTc                                                          <  2:286626/251‑1 (MQ=255)
                                                                                                                                                                                           gtgtaaGCCGGCGTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCCAAGGATACCGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGTTTTTGCATGCCGAGCAGATTTAATTCATAGATTCGGGAGCCTAGCCAGAGGGCACGTTAATTCCGTCTTTGCAAGATCTTCCGGTGCACcg                                                      <  2:398851/247‑1 (MQ=255)
                                                                                                                                                                                                      cTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGtt                                           >  1:373183/1‑251 (MQ=255)
                                                                                                                                                                                                        ttCACGGGCAGGCGTGGTGGAAAATGCCCTTTATTGGTAGCGGAAACGAAAAACAGGCCATTGCTCACATTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTCCATGCCACAAGGCAAAGAACGGACTGAGCAAAGTGGTTGTGCATGCCGAGCAGATGGAATTGAGAGAGTCTGGAGACTATCCATATGGCACGTTATTTCCGTCGTCTCAAGGTCTGcc                                                               <  2:399614/229‑1 (MQ=255)
                                                                                                                                                                                                        ttCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGCAACGAAAAACACGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTAAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATCGATTCTGGCGACTAGACATATGGCACGTTATTTCCGTCGTCGCACGTTCTGcc                                                               >  1:399614/1‑229 (MQ=255)
                                                                                                                                                                                                                cAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGACGGCGTTCAag                                       >  1:302949/1‑245 (MQ=255)
                                                                                                                                                                                                                 aGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCg                                 >  1:289590/1‑250 (MQ=255)
                                                                                                                                                                                                                           tCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCCTTACTCACAGTATAAAGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCa                                         >  1:135027/1‑232 (MQ=255)
                                                                                                                                                                                                                           tCAAATGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTGCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAAGGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATAGGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCa                                         <  2:135027/232‑1 (MQ=255)
                                                                                                                                                                                                                               aTGCCCGTTATTGTTAGCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGGCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCg                  <  2:363364/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ttATTGTTAGCGGAAACGAAAAACAGGCCATTACTAACAGTATAACGGTCGGCAGCCGCATAACCGTTCAGGGGTTCATTTAATGCCACAAGGCAAAGAACGGACTGAGCAACATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGAACGTTATTTCCGTCGTCGCAAGATCTGCCGTTTAACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCt           >  1:221197/1‑251 (MQ=255)
                                                                                                                                                                                                                                               gCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAGAAACTa  <  1:116587/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
CCGTGAGCGTCGCGATGATTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGACACGAAAACCAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTA  >  NC_000913/4425444‑4425933

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: