Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,959,564 C→T W955* (TGG→TGA recC ← exodeoxyribonuclease V subunit RecC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,959,5640CT80.7% 134.8 / 29.6 62W955* (TGG→TGArecCexodeoxyribonuclease V subunit RecC
Reads supporting (aligned to +/- strand):  ref base C (7/5);  new base T (25/25);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

ACCAGCAATGGCGCGGACATTCCTTCACGATACCCCTCAATCAGTTGTGAGAGGTAATGCAAAGCCTGTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAACCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCT  >  NC_000913/2959319‑2959803
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
acCAGCAATGGCGCGGACATTCCTTCACGATACCCCTCAATCAGTTGTGAGAGGTAATGCAAAGCCTGTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAAtcagcc                                                                                                                                                                                                                                            >  2:433592/1‑251 (MQ=255)
 acAGCAATTGCGCGGACATTCCTTCACGATACCCCTCAATCAGTTGTGAGAGGTAATGCAAAGCCTGTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCAGCCGGCTGCACCTGCGGCAAccagcca                                                                                                                                                                                                                                           <  1:170807/250‑1 (MQ=255)
      aaTGGCGCGGACATTCCTTCACGATACCCCTCAATCAGTTGTGAGAGGTAATGCAAAGCCTGTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAACCAGCCAGTTAt                                                                                                                                                                                                                                      <  2:77622/251‑1 (MQ=255)
         ggCGCGGACATTCCTTCACGATACCCCTCAATCAGTTGTGAGAGGTAATGCAAAGCCTTTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAACCAGCCAGTTATCTg                                                                                                                                                                                                                                   <  2:20249/251‑1 (MQ=255)
                 cATTCCTTCACGATACCCCTCAATCAGTTGTGAGAGGTAATGCAAAGCCTGTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGGCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGt                                                                                                                                                                                                                           <  1:150951/251‑1 (MQ=255)
                  aTTCCTTCACGATACCCCTCAATCAGTTGTGAGAGGTAATGCAAAGCCTGTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGtt                                                                                                                                                                                                                          >  1:456060/1‑251 (MQ=255)
                                cccnTCAATCAGTTGTGAGAGGTAATGCAAAGCCTGTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAGAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCaa                                                                                                                                                                                                            <  2:439599/251‑1 (MQ=255)
                                ccccTCAATCAGTTGTGAGAGGTAATGCAAAGCCTGTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCa                                                                                                                                                                                                             >  1:461859/1‑250 (MQ=255)
                                          aGTTGTGAGAGGTAATGCAAAGCCTGTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATAct                                                                                                                                                                                                  >  2:9236/1‑251 (MQ=255)
                                               tgAGAGGTAATGCAAAGCCTGTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGcc                                                                                                                                                                                             >  2:387631/1‑251 (MQ=255)
                                                gagaGGTAATGCAAAGCCTGTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAACCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGcc                                                                                                                                                                                             <  1:419111/250‑1 (MQ=255)
                                                   agGTAATGCAAAGCCTATTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCgg                                                                                                                                                                                         >  2:399485/1‑251 (MQ=255)
                                                       aaTGCAAAGCCTGTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTgg                                                                                                                                                                                     >  2:263251/1‑251 (MQ=255)
                                                       aaTGCAAAGCCTGTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTgg                                                                                                                                                                                     <  2:456060/251‑1 (MQ=255)
                                                           cAAAGCCTGTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGAGAGCTCAATTTCCATACTCTGCCCCGGCTGGCGAc                                                                                                                                                                                 >  2:73542/1‑251 (MQ=255)
                                                           cAAAGCCTGTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTCCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGGAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGAc                                                                                                                                                                                 <  1:378781/251‑1 (MQ=255)
                                                              aGCCTGTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAgg                                                                                                                                                                              <  2:167577/251‑1 (MQ=255)
                                                              aGCCTGTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTGTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAACCAGCCAGTTATCTGCACACCGTTGCAGGCGATATCAATTTCCATACTCTGCCCCGGCTGGCGACAgg                                                                                                                                                                              <  2:281292/251‑1 (MQ=255)
                                                               gCCTGTTCCGCTGCAAGCGACGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGc                                                                                                                                                                             <  2:137617/251‑1 (MQ=255)
                                                                   gTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATg                                                                                                                                                                         >  1:366209/1‑251 (MQ=255)
                                                                   gTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTACCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACATGCGGCGATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATg                                                                                                                                                                         >  1:29015/1‑251 (MQ=255)
                                                                      cGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGAct                                                                                                                                                                      >  1:382021/1‑251 (MQ=255)
                                                                       ggCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGActc                                                                                                                                                                     >  1:328292/1‑251 (MQ=255)
                                                                        gCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGActct                                                                                                                                                                    >  1:198828/1‑251 (MQ=255)
                                                                             aaGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAACCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCgg                                                                                                                                                               >  1:488627/1‑251 (MQ=255)
                                                                                      aaaaCGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAgctgct                                                                                                                                                      <  1:400110/251‑1 (MQ=255)
                                                                                          cGCCCCTCGCCGTCTTTGCGTAGAAAAAGGCGACGTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCAt                                                                                                                                                  <  2:126126/251‑1 (MQ=255)
                                                                                             cACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCAGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCAtctc                                                                                                                                               <  2:7352/251‑1 (MQ=255)
                                                                                                ccGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTg                                                                                                                                            <  2:424404/250‑1 (MQ=255)
                                                                                                    cGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCAc                                                                                                                                        <  2:328292/251‑1 (MQ=255)
                                                                                                          tGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAACCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTgtgt                                                                                                                                   <  2:38495/250‑1 (MQ=255)
                                                                                                            cGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTc                                                                                                                                >  2:214938/1‑251 (MQ=255)
                                                                                                              tAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTccc                                                                                                                              >  1:13104/1‑251 (MQ=255)
                                                                                                               aGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTccc                                                                                                                              <  1:387631/250‑1 (MQ=255)
                                                                                                                     aGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAAt                                                                                                                       <  2:29015/251‑1 (MQ=255)
                                                                                                                             ttCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCa                                                                                                                >  1:265574/1‑250 (MQ=255)
                                                                                                                              tCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAACCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCaa                                                                                                               >  1:463458/1‑250 (MQ=255)
                                                                                                                                              tGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTaa                                                                                               >  2:92617/1‑250 (MQ=255)
                                                                                                                                                       aGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGcc                                                                                     <  2:358277/251‑1 (MQ=255)
                                                                                                                                                            aGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCg                                                                                <  1:92617/251‑1 (MQ=255)
                                                                                                                                                                  cccAGCCAAAGTTGCATTCCCTGGGCCACAATTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAgc                                                                           <  2:84572/249‑1 (MQ=255)
                                                                                                                                                                        cAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCAGCTCGGAAGCGGCGGa                                                                     >  2:348814/1‑250 (MQ=255)
                                                                                                                                                                                 aTTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTc                                                           <  1:121942/251‑1 (MQ=255)
                                                                                                                                                                                     ccTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGca                                                       <  1:433592/251‑1 (MQ=255)
                                                                                                                                                                                      cTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCCAAAGAGCCATACGGTAAATCCCCTGCCACTCGAAAGCGCCGGAACCAGCGATTGGcat                                                      >  2:370905/1‑251 (MQ=255)
                                                                                                                                                                                         ggCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGcatcat                                                   <  2:13104/251‑1 (MQ=255)
                                                                                                                                                                                               acTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCt                                             <  1:416968/251‑1 (MQ=255)
                                                                                                                                                                                                 ttAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACCAGCGTTCGGCATCATCCTGCTc                                            >  2:109775/1‑250 (MQ=255)
                                                                                                                                                                                                     taaAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAACCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAAcc                                        >  1:449763/1‑250 (MQ=255)
                                                                                                                                                                                                     aaaaGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGTCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCa                                       <  2:321747/250‑1 (MQ=255)
                                                                                                                                                                                                              aCGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCa                              <  1:404800/251‑1 (MQ=255)
                                                                                                                                                                                                                gCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAAt                            >  1:319126/1‑251 (MQ=255)
                                                                                                                                                                                                                gCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAAt                            >  2:265043/1‑251 (MQ=255)
                                                                                                                                                                                                                 ccAGCGCAACAGGCCATCCGGCTGCACATGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCAGCATCATCCTGCTCAACCAGTGCATTctat                            >  1:260066/1‑247 (MQ=255)
                                                                                                                                                                                                                  cAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATa                           >  2:274483/1‑250 (MQ=255)
                                                                                                                                                                                                                   aGCGCAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAAc                         <  2:40890/251‑1 (MQ=255)
                                                                                                                                                                                                                      gcAACAGGCCATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAActgc                      >  2:291277/1‑251 (MQ=255)
                                                                                                                                                                                                                             gCCATCCGGCTGCACCTGCGGCAACCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCtgattg                >  1:351191/1‑250 (MQ=255)
                                                                                                                                                                                                                              ccATCCAGCTGCACATGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGATGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGAGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCtgattgat              >  2:211901/1‑251 (MQ=255)
                                                                                                                                                                                                                               cATCCGGCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTgattgatt             <  1:61323/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ggCTGCACCTGCGGCAATCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATa        <  2:366209/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ggCTGCACCTGCGGCAACCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGAttgatttgat         >  2:46719/1‑250 (MQ=255)
                                                                                                                                                                                                                                    ggCTGCACCTGCGGCAACCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGAACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGAttgatttgct         >  1:487527/1‑248 (MQ=255)
                                                                                                                                                                                                                                          aCCTGCGGCAACCAGCCAGTTATCTGCACACAGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAACAGCGCCATACGGAAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCt  >  1:118178/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
ACCAGCAATGGCGCGGACATTCCTTCACGATACCCCTCAATCAGTTGTGAGAGGTAATGCAAAGCCTGTTCGGCTGCAAGCGGCGGAAAACGCCACTCGCCGTCTTTGCGTAGAAAAAGGCGACTTTCACCATTACCACCGCTGGCACAGTAGACAAGGTGTTCCAGCCAAAGTTGCATTCCCTGGGCCACACTTAATAAAGAGGGACGCCAGCGCAACAGGCCATCCGGCTGCACCTGCGGCAACCAGCCAGTTATCTGCACACCGTTGCAGGCGAGATCAATTTCCATACTCTGCCCCGGCTGGCGACAGGCAATGACTCTGTCGGCAAGCTGCTGCATCTCCTGGCACTGTGTTTCCCAGAAAATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCT  >  NC_000913/2959319‑2959803

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