Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,626,586 T→C G225G (GGT→GGC)  radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,626,5860TC97.2% 204.4 / NA 73G225G (GGT→GGC) radADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (31/38);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGC‑TGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAA  >  NC_000913/4626345‑4626830
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
gCCTTCCGACTGATAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCGGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGCTCCAGGTGATGCATATGGCGGATTTACAGTCATCGCCTGGCAGCGTGGCGGAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAATATGGCTCgctggca                                                                                                                                                                                                                                              <  2:308959/251‑2 (MQ=255)
 ccTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCgctggcgg                                                                                                                                                                                                                                             <  1:360040/251‑3 (MQ=255)
   tgcagACTAACAAGCTCAATATGTTGTCGGATACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGACGAGCAAACGAAGATGATGGAAATTGACTCGATCCAGGTGATGCATATGGGGGATGTACAGTCATCGCCTGGCAGCGTGGCGCTGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGATGGGG‑CGATTGT‑CCTAGTTGGACCCGTAACCAACGATGGGtggctggcgggg                                                                                                                                                                                                                                           >  2:355585/5‑247 (MQ=255)
     ccGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGccc                                                                                                                                                                                                                                        >  1:558923/1‑251 (MQ=255)
           gACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAAACAAAGATGGCTCGCTGGC‑GGGCCCTAAAGt                                                                                                                                                                                                                                  >  1:476342/1‑251 (MQ=255)
                   cAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTg                                                                                                                                                                                                                              <  1:492655/247‑1 (MQ=255)
                       atGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACt                                                                                                                                                                                                                       <  1:431334/250‑1 (MQ=255)
                       atGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGGGCTGGAACCCTg                                                                                                                                                                                                                      >  1:578348/1‑251 (MQ=255)
                         gttgtCTGAAACCAGCATCGAGCAGATCTGCCGGATGGCCGAAGAAGAGCATCCGAAGCGGATGTGAATTGACTCGAGCCAGGTGATGCATATGGCGGATGTACAGGTATCGCCTGGCATCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGATGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTa                                                                                                                                                                                                                    <  2:328378/251‑1 (MQ=255)
                               ggAAACCAGCATCGAGCAGATATGCATGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGGGGCGCAAGGGCGTGAAACAGCGGCATATCTGACGCGCTTCGCCAAAACGCGCGGGGGGG‑CGATTGT‑CAAGGGTGGGCACGTAAACAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGAc                                                                                                                                                                                                               >  1:365655/1‑250 (MQ=255)
                                       gCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTg                                                                                                                                                                                                      <  1:105339/251‑1 (MQ=255)
                                            gagcAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCAGATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCtttt                                                                                                                                                                                                 <  1:26883/248‑1 (MQ=255)
                                                      gCCTGATTGCCGAAGAAGAGCAACCGATGCTGATAGTATTTGTCTAGATTCTGTTGCTGCATATGGCGTATGGTCAGTCTGCGCCTGGCAGCGTTGCGCCGGTGCGTGTAACGGCGTCTTATCTGACTCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTATCCAAAGATGGTTCGCTGGC‑GGGGTATGAAGTGCTGGAACACTGTATCTACTGTTCGGTGCttttttctga                                                                                                                                                                                           <  2:95665/247‑7 (MQ=255)
                                                      gCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGGGATGCATATGGCGGATGTACAGGCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGc                                                                                                                                                                                          <  2:371001/248‑1 (MQ=255)
                                                         tGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGAGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCCGGAACACTGTATCGCCTGTTCCGTGCTTTTGGATGGCGATGcc                                                                                                                                                                                    >  1:285309/1‑251 (MQ=255)
                                                                 gaagaaGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCg                                                                                                                                                                            >  1:286672/1‑251 (MQ=255)
                                                                   agaagaTCAGCCGAAGTTGGTGGTAATTGACTCGTGGCAGGTGATGCATATGGCGTATGTACTGGCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAACCGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGta                                                                                                                                                                          <  2:223582/251‑2 (MQ=255)
                                                                   agaagaGCAGCCGAAGCTGTTGGTAATTGACGCGTTCCAGGTTAGGCATATGGCGGAGGTCTAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGTCGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCGCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGGTCGGTGCTTTTGGATGGCGATGCCGTCTCCCGtt                                                                                                                                                                          <  2:6351/251‑1 (MQ=255)
                                                                      agagTAACCGAAGCGGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATTTTCAGTCAGCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGTCGCGCGTCGCCAAAACGCGCGGTGTGG‑CTATTGG‑CATGGTTGGACACGTAACCAAATATGGCTCGCTGGC‑GGGCACTAAAGTGCGGGAACACTGTACCGACTGTGCGGTGCTGTTGGATGGCGACGCTGACTCCAGTTTTc                                                                                                                                                                       <  2:388471/251‑1 (MQ=255)
                                                                           aaCCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCAct                                                                                                                                                                  <  1:214730/251‑2 (MQ=255)
                                                                              cGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAATATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCActtta                                                                                                                                                               <  1:394882/251‑5 (MQ=255)
                                                                               gAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCActttac                                                                                                                                                              <  1:633276/251‑6 (MQ=255)
                                                                                     gATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGc                                                                                                                                                         <  1:505656/250‑1 (MQ=255)
                                                                                       gggtaatggacgtgATCCAGGTGATGCATATGGCGGAGGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGGGTGAAACGGTGGCTTATCTGACGCGCTACGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCGGTC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCAt                                                                                                                                                      <  2:155283/238‑1 (MQ=255)
                                                                                           aatggatgcgagccAGGTGATGCAAATGGCGGATGGACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCGGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGTCTCCCGTTTTCGCACGTTAGGCAGCCATaaaa                                                                                                                                                  <  2:32440/239‑1 (MQ=255)
                                                                                            aTTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATaaaaa                                                                                                                                                 <  1:509558/251‑1 (MQ=255)
                                                                                            aTTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAAAACTGTATCGACTGTTCGGGGCTTTTGGATGGCGGTGCCGACTCCCGTTTTCGCACTTTACGCAGCCATaaaa                                                                                                                                                  >  1:508403/1‑250 (MQ=255)
                                                                                              tGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAAc                                                                                                                                                <  1:277807/250‑1 (MQ=255)
                                                                                                  ccGATCCGGGTGATGCAAATGGCTGATTTACAGTCATCGCCGGGCGGCGTGGCGCAGGTGCGTGAAACGGCGGTTTATCTGACGCGTGTCGTCAAAACGCGCGGTGTGG‑CGATTTG‑CATGGTTGGACACGTAACCAAAGCTGGTCGGCGGGCGTGGTC‑TAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCtt                                                                                                                                           <  2:376233/250‑1 (MQ=255)
                                                                                                        cAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACCCTGTATCGACTGTTCGGGGCTTTTGGATGGCGGTGCCGACTCCCGTTTTCGCACTTTACGCCGCCCTAAAAACCGCTTCGgcgg                                                                                                                                     >  1:10935/1‑250 (MQ=255)
                                                                                                           gggaggCATATGGTGTATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAACCGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGG‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGGGCGGGAACACTGTATCGACTGTTCGGGGCTTTTGGATGGCGATGCCGACTCCCGTGTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGt                                                                                                                                  <  2:654811/246‑1 (MQ=255)
                                                                                                            tGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTTTCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CCTGGTTGGACACGTAACCCAAGATGGCTCGCTTGC‑GGGCCCTAAAGTGCTGGAACACAGTATTCACTGTTCGGGGCTTTTGGATGGCGATGCCGCCTCCCGTTTTCGCACTTTTCGCCGCAATTAAAACCGCTTctgcacgggg                                                                                                                                 >  2:99834/1‑242 (MQ=255)
                                                                                                             gATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTg                                                                                                                                 <  1:202487/250‑1 (MQ=255)
                                                                                                             gATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGa                                                                                                                                <  1:6086/251‑1 (MQ=255)
                                                                                                                  atatGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGAAACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGTCTGTTCGGTGCTTTTGGATGGCGATGCCGAATCCCGTTTTTGCACTTTAAGCAGCCCTAAAAAACGCTTCGGCGAGGTGAATg                                                                                                                             >  1:300109/1‑249 (MQ=255)
                                                                                                                       gCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGCTGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGCAAACTGTATCGGCTGTTCGGTGCTTTTGGATGGCGCTGCCGACTCCCGTTTTCGCACTTTACGCAGGCATAAAAAACGCTTCGGCGCGGTGAATGAGCTggg                                                                                                                      >  1:410471/1‑251 (MQ=255)
                                                                                                                        cGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACTGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATTGCTCgctggc                                                                                                                                                                                                                                               <  2:446261/130‑1 (MQ=255)
                                                                                                                        cGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCCAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCgctggc                                                                                                                                                                                                                                               >  1:446261/1‑130 (MQ=255)
                                                                                                                           aTGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTc                                                                                                                  >  1:72568/1‑251 (MQ=255)
                                                                                                                           aTGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGGGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTc                                                                                                                  >  1:527580/1‑251 (MQ=255)
                                                                                                                              tACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTc                                                                                                               >  1:492594/1‑251 (MQ=255)
                                                                                                                               aCAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTc                                                                                                               <  1:119060/250‑1 (MQ=255)
                                                                                                                                       cGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGCTGCCGACTCCCGTTTTCGCACTTTACGCCGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGAcc                                                                                                      >  1:221082/1‑251 (MQ=255)
                                                                                                                                        gCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTTGACACGTAACCAAAGATAGCTCGCTGGC‑GGGACCTAAAGTGCTGGAACACTGTATCGGATGTTCGGTGCTTTTGGAAGGCGCGGCCGCCTCCCGTTTTCGCCCTTTAAGCCGCCATAAAAACCGCTTCGGCGCGGGGCAAGGGCTGGGCGTCTTcactagtaccc                                                                                                     >  2:458050/1‑241 (MQ=255)
                                                                                                                                        gCCTGGCAGCGGGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGA‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGGACACTGGATCGACTGTTCGGTGGTTTTGGATGGCGATGCCGACCCCCGTTTTCGCCCTTTACCCCGCAAAAAAAACACCTTCGGCGCGGTGAATGAGCAGGGCGTCTTccgcatggccc                                                                                                     >  2:348741/1‑241 (MQ=255)
                                                                                                                                         ccTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGGGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGGGCTGGGACACTGTATCGACTGGTCGGGGCTTTTGGATGGCAATGCCGACTCCCGTTTTCGCACTTTACGCAACCATAAAAACCGCTTCGGCGCGGGGGATGAGGTGGGCGTCTTCGCGAttgccg                                                                                                     >  2:545849/1‑245 (MQ=255)
                                                                                                                                            ggcagggggtggagggggCGTGAAACGGCGGCTTATCTGACGAGATTCCGCAAAACGCGTGGTGTGG‑GGAT‑GTGAATGGTTGGACACGTAACCAAAGATGTCTCGCTGGT‑GGGACCTAAAGTGCTGGAAAAATGTATCGACTGTTCTGTGCTTTTTGAAGGCGATGGCGACTCCCGTTTGCGTACTTTACGCAGCAATAAAACCCGCTGCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGCCCGaaaa                                                                                                 <  2:265528/234‑3 (MQ=255)
                                                                                                                                             gcaGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAg                                                                                                <  1:120830/249‑1 (MQ=255)
                                                                                                                                                           gTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTGGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGt                                                                                  <  1:296828/251‑1 (MQ=255)
                                                                                                                                                             gCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCa                                                                                <  1:545849/251‑1 (MQ=255)
                                                                                                                                                                          gctattctGACGCGCTTCGCCAAAACGCGCGGGGTGG‑CGATAGT‑CATGGGTGGACAAGTAACCAAAGATGGCTCGCTGTC‑GGGCCCTAAGGTGCTGGAACACTGTATCGACTGGTCGGTGCTTTTGGAGGGCGATGCCGACTCCCGTTTTCGCACTGTACGCAGCCATAAAAACCGCTGCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGc                                                                   <  2:6321/244‑1 (MQ=255)
                                                                                                                                                                             tatctGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCa                                                                  <  1:99834/245‑1 (MQ=255)
                                                                                                                                                                              atctGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCa                                                                                                                                                       <  2:371176/160‑1 (MQ=255)
                                                                                                                                                                              atctGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCa                                                                                                                                                       >  1:371176/4‑163 (MQ=255)
                                                                                                                                                                              atctGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAAtt                                                               >  1:338213/4‑251 (MQ=255)
                                                                                                                                                                               tctGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAAttt                                                              >  1:162822/3‑251 (MQ=255)
                                                                                                                                                                                 tGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTc                                                            >  1:648762/1‑251 (MQ=255)
                                                                                                                                                                                       gCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGt                                                      <  1:257533/251‑1 (MQ=255)
                                                                                                                                                                                         gtggCCAAAACGCGCGG‑GTGGGCGATTGG‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCGTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTGCTTAAGTcg                                                    <  2:545962/248‑1 (MQ=255)
                                                                                                                                                                                             ccAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCgcgga                                                <  1:498873/251‑2 (MQ=255)
                                                                                                                                                                                             ccAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGGGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCgcggg                                                >  1:453332/1‑250 (MQ=255)
                                                                                                                                                                                                 aactcgtggTGTGG‑AGATTGG‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCATAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGAGGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATg                                            <  2:241219/244‑1 (MQ=255)
                                                                                                                                                                                                           gtgG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCt                                   <  1:140892/250‑1 (MQ=255)
                                                                                                                                                                                                           gtgG‑CGATTGG‑CATGGTTGGACACGTAACCAAAGATGGCTCGCGGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCAACTGTTCGGGTCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTTCGCAGCCATAAAAACCGCTTCGGCGCGGGGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGGGAAGACAGCAAACCTTCGGCAATTTTCTTAAGTGGCGGGGCTGAGGTGACCTc                                  >  2:87298/1‑251 (MQ=255)
                                                                                                                                                                                                              ggtgATTTT‑CATGGTTGGATACGTACCCAAAGACGCATCGCTGGC‑GGGCCCTAATGTGCTGGAACACTGTAGCGACGGTTCGGTGATTGTGGATGGCGATGCAGTCTCCCTTCTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGTGCTGTGAATGAGCTTGGCTTCTTCGCGATGACCGAACCGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAAGTTTCTTAAGGCGCGTAGATTAAGTGACCTcc                                 <  2:406871/247‑1 (MQ=255)
                                                                                                                                                                                                                 gATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGGCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCggcc                             >  1:363401/1‑249 (MQ=255)
                                                                                                                                                                                                                   ttGT‑CATGGTTGGTCACGTAACCAAAGCTGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCGTTTGGATGGCGATGCCGACTCCCGTATTCGCACTTTACGCAGCCATACAAACCACTTCGGCGCGGTGAATGAGCTCGGCTTATTCGCGTTGACAGCAAAGGGGCTGCGGGAGGGCAGCGACCATCAGGCAGTTTTCTGAAGTGGGGGAGATGAAGAGAACTCCggagac                           >  2:50588/1‑247 (MQ=255)
                                                                                                                                                                                                                     gT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTc                         >  1:265148/1‑251 (MQ=255)
                                                                                                                                                                                                                           ggttggacACGTAACCAAAGATGGCTCGCTGGC‑GGGACCTAAAGTGCGGGAACACTGTATCGCCTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGa                    <  2:179031/244‑1 (MQ=255)
                                                                                                                                                                                                                               ggacACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAACTGACCTCCGGCAGCTCCGTGAtggt                >  1:107367/4‑251 (MQ=255)
                                                                                                                                                                                                                                     gTAACAAAAGATGACTCGCTGGC‑GGGCCCTCAACTGCTGGCACACTGTATCGACTGCTCGGCGCTTTTGGATGGCGATGCCGACTCCCGTTTTAGCACCTTACGAAGCCATTAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTTGCGATGACCGAACAGGAGCTGCGTGAAGTCAGCAACCCTTCGGCCATTTTATTAAGTCGCGGGGATGGAGTGATCTCCggacgccccgcagggcgggtggg          >  2:323335/1‑230 (MQ=255)
                                                                                                                                                                                                                                      tAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATACAAACCTCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGGCACCCTTCGGCCATTTTCTTAAGTCGGGGAGAGGAAGTGGCCTTCGGAAGCTCCGTGATGGTGGTGTgg         >  2:519503/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aGATGGCTCGTTGGC‑GGGCCCTAAAGTGCGGGAACACGGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAATCCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTGCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAATCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGggaaggaa  <  2:492594/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
GCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGC‑TGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAA  >  NC_000913/4626345‑4626830

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑