Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,425,696 C→A N59K (AAC→AAA priB → primosomal replication protein N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,425,6960CA96.4% 182.5 / NA 57N59K (AAC→AAApriBprimosomal replication protein N
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (26/28);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TGATTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGACACGAAAACCAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCAC  >  NC_000913/4425459‑4425939
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
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)
                                                                                                                                                                                                                               agcgtaaacGAAAAACAGGTCATTACTCCCAGTATAACGGTCGTCAGTCGCATAACCTTTCAGGGGTTCATGTCATGCCACAAGGCAAAGAACTGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTTAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCCCCGCGTAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACt         <  1:398912/244‑1 (MQ=255)
                                                                                                                                                                                                                                gCGGAAACGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAGAAACTa        <  1:116587/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ggaaacGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTAGCATACACGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGCCTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTCTTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTAc       >  2:437039/5‑250 (MQ=255)
                                                                                                                                                                                                                                   gaaacGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGTGTTCATTTCATGCCACAAGGCAAAGAACGGACTGATCAAAATGGTTTTGCATGCCGAGCATATGGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATGTTCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGGTCAAGAGATCTACGATAAAGATATCGCTACGCTGAAAAACTACAt     <  1:183925/248‑1 (MQ=255)
                                                                                                                                                                                                                                   gaaacGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGCAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACAt     >  2:72591/4‑251 (MQ=255)
                                                                                                                                                                                                                                     aactaaaaACAGGCCATTACTCACTGTATAACTGTCGTCATTCTCATACCCTGTCAGGTGTTCATTTCATGCCACAAGGCAAAGAACGGACTGATCAAAATGTTTTTGCATGCCGAGCAGTTTGAATTGATAGATTCTGTAGACTAGCCAGTTGGCACGTTATTTCCGTCGTTGCAAGTTCTTCCGTTTCACCGCGGAAGTCTTTCATTAGATCTACTGTAAAGATATCGCTACGCTGAAAAACTACATCa   <  1:256548/247‑1 (MQ=255)
                                                                                                                                                                                                                                     aacGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCa   <  1:385227/250‑1 (MQ=255)
                                                                                                                                                                                                                                      acGAAAAACAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCAc  >  1:166939/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
TGATTTCGCAAACGAAACCGCTGATGATGCTGAAGCTGGGGATTCTGAAGAGTAATTTCTGATGACCAACCGTCTGGTGTTGTCCGGCACCGTGTGCAGGGCTCCCCTTCGAAAGGTCAGTCCATCAGGAATTCCTCACTGCCAGTTCGTGCTTGAGCATCGTTCTGTGCAGGAGGAAGCCGGCTTTCACCGGCAGGCGTGGTGTCAAATGCCCGTTATTGTTAGCGGACACGAAAACCAGGCCATTACTCACAGTATAACGGTCGGCAGTCGCATAACCGTTCAGGGGTTCATTTCATGCCACAAGGCAAAGAACGGACTGAGCAAAATGGTTTTGCATGCCGAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCAC  >  NC_000913/4425459‑4425939

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: