Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,626,595 T→G A228A (GCT→GCG)  radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,626,5950TG98.6% 199.6 / NA 71A228A (GCT→GCG) radADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (31/39);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

CCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGC‑TGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCAC  >  NC_000913/4626350‑4626839
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
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)
                                                                                                                      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                                                                                                                                                                >  1:371176/4‑163 (MQ=255)
                                                                                                                                                                         atctGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCa                                                                                                                                                                <  2:371176/160‑1 (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)
                                                                                                                                                                                                                                         gatgtctCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAAc          <  2:453332/245‑1 (MQ=255)
                                                                                                                                                                                                                                              ctCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGGCTCCCGTTTTCGCACTTTACGCAGCAATAACAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGTTGGTGGTGTGGGAAGGAACGCGTc     >  1:288140/2‑251 (MQ=255)
                                                                                                                                                                                                                                                cGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTcc    <  1:560386/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 gttggc‑gggcCCTAAAGTGCTGGAACCCTGTATCGACGGTTCGTGGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTGCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTccgt  <  2:286672/249‑3 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
CCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGC‑TGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCAC  >  NC_000913/4626350‑4626839

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