Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,626,820 A→G V303V (GTA→GTG radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,626,8200AG98.0% 147.7 / NA 49V303V (GTA→GTGradADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (20/28);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

ATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTT  >  NC_000913/4626582‑4627058
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                
aTGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACAATGTATCGCCTGTTCGGAGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCCTAAAAACCGCTTCGACGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGGTGGAGGTGTGGGAGGGAAcg                                                                                                                                                                                                                                     >  1:425064/1‑251 (MQ=255)
  ggctcggtggcgggcCCTAAAGGGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGATTCACGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGGCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAAcgcg                                                                                                                                                                                                                                   <  2:302943/244‑1 (MQ=255)
              ggcCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCtggt                                                                                                                                                                                                                       >  1:91866/1‑251 (MQ=255)
              ggcCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTGTTCGCACTTTACGCAGCCATAAAACCCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCtggt                                                                                                                                                                                                                       <  1:120539/251‑1 (MQ=255)
                ccctaaaGTGCTGGAACCCTGTATCGACTGTTCGGTGCTTTTGGATGTCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGTTCCGTGATGGTGGTGTGTGAAGGAACGCGTccgtttctggtgg                                                                                                                                                                                                                     <  1:171858/247‑12 (MQ=255)
                     aaGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGAtt                                                                                                                                                                                                                >  1:343875/1‑251 (MQ=255)
                            ggAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGgc                                                                                                                                                                                                           <  1:404747/249‑1 (MQ=255)
                             gAACACTGTATCGACTGGTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCt                                                                                                                                                                                                        <  1:161784/251‑1 (MQ=255)
                                          aCTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACTGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCg                                                                                                                                                                                           <  1:39303/251‑1 (MQ=255)
                                                   tGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAACAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCCGCTCCGTGCTGGTGGTGTGGGAAGGACCGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCCCTCGGTGCTGGCg                                                                                                                                                                                  >  2:212386/1‑251 (MQ=255)
                                                           gATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGaatccgcg                                                                                                                                                                          >  1:446654/1‑245 (MQ=255)
                                                              gttgttGCCGACTCCCGTTTTGTCACGTTACGCGTCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATTACCGAACAGGGGCGGCGTGAAGTTAGCAACCCTTCGGCAATTTTCTTAAGGCGCGGAGATGAAGTTACCGCCGGCAGCTCCGTGATGGTGGTGTGAGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATggttaatccgcgcc                                                                                                                                                                        <  2:392566/245‑13 (MQ=255)
                                                                     ccGACTCCCGTTTTCGCACTTTACGCAGCCATAAAACCCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGCTCACTCGATGATGGCGAATCCGCGCCGCGTGGcc                                                                                                                                                                >  2:280454/1‑250 (MQ=255)
                                                                     ccGACCCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCGCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCa                                                                                                                                                                <  2:76635/251‑1 (MQ=255)
                                                                       gcctccctgtttCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGGGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGt                                                                                                                                                              <  2:406707/242‑1 (MQ=255)
                                                                             cGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCt                                                                                                                                                        <  1:9160/251‑1 (MQ=255)
                                                                                tctccaccgttcctcagcaaaaaaaaCCGCTTCGGCGCGGTGAAGGAGCTGGGTGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCACCCGGTCGGCAAGTTTCGTAAGGCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGGCCGTTGCTGGTGGAGATTCAGGCGCTGGGAGATCACTTGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTgg                                                                                                                                                      <  2:345242/229‑1 (MQ=255)
                                                                                     cactttacgcggcaatagaaaCCGCTTCGGCGCGGTGCATGAGCTGGTCGTCATCGCGATGACCGAACAGGGGCTGCGTGAAGTAAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGAGATGATGGAGGTGGAAAGCACGCGTCCGTTGCTGGAGGAGCTTCAGGCGCTGGACGATCACTCTATAATGGCGCATCAGCGCCGCATGGCAGTGGGGCTGGGACac                                                                                                                                                 >  2:170939/19‑249 (MQ=255)
                                                                                            cgcAGCCATAAAAACCGCTGCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCt                                                                                                                                         <  2:129521/251‑1 (MQ=255)
                                                                                               aGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGGTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGGGGCAGTGGGTCTGGCACAAAACCGTATGGc                                                                                                                                      >  2:764/1‑251 (MQ=255)
                                                                                                        aaCCGCTGCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCctgct                                                                                                                             <  1:357774/251‑1 (MQ=255)
                                                                                                         aCCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCctgctg                                                                                                                            >  1:59988/1‑251 (MQ=255)
                                                                                                                cGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCtgtc                                                                                                                      <  1:292824/250‑2 (MQ=255)
                                                                                                                      ggTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCAcc                                                                                                               <  1:359403/251‑1 (MQ=255)
                                                                                                                           aTGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGCTTCAGGCGCTGGTCGATCACTCGATGATGGCGTATCCGCGCCGCGTGGCAGTGGGGCTGGAACACAACCGTATGGCAATCCTGCTGGCTGCCTTGCACCGtacc                                                                                                          >  2:413595/1‑248 (MQ=255)
                                                                                                                                  gggCGTCGGCGCGATGAACGAACAGGGGATGCGTGAAGTCAGCAACCCTTCGGCAATTTTCGTAAGTCGCGGAGATGAAGTGACCGCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACggtggt                                                                                                    <  2:51192/250‑1 (MQ=255)
                                                                                                                                   ggCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCTCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTTCACCGTACCGGTGGTCt                                                                                                  >  2:214503/1‑251 (MQ=255)
                                                                                                                                    gggccTTCGCGATGACCGAACAGGGGCTTCGGGAAGTCAGCAACCCTGCGGCAATTTTCGTAAGTCGCGGATATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTg                                                                                                 <  2:423422/247‑1 (MQ=255)
                                                                                                                                              gATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCg                                                                                       >  1:72453/1‑251 (MQ=255)
                                                                                                                                                    cGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAgg                                                                                 <  1:96014/251‑1 (MQ=255)
                                                                                                                                                      aaCAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGAt                                                                               >  1:233169/1‑251 (MQ=255)
                                                                                                                                                      aaCAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTAGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGGGGAGATTCAGGCACTGGTCGGTCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGCCTGGCAATCCTGCTGGCTGTCTTTCACCGGCAGGGTGGTCTGCAAATGGCCGATCAggct                                                                               >  1:9517/1‑249 (MQ=255)
                                                                                                                                                           gggCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGa                                                                                      >  1:213920/1‑239 (MQ=255)
                                                                                                                                                           gggCTGCGTGAAGTCAGCAACCCTGCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGa                                                                                      <  2:213920/239‑1 (MQ=255)
                                                                                                                                                                gCGTGAAGTCAGCATCCCTGCGGTAATTTTCTTAATGCGCGGAGATGAAGTGACCTCCGGCAGCTCCGGGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTtgt                                                                       <  2:347776/249‑1 (MQ=255)
                                                                                                                                                                   ggAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGGTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACg                                                                  <  2:413023/250‑1 (MQ=255)
                                                                                                                                                                       gTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGAGGAGATTCAGGAGCTGGTCGATCACTCGATGATGGCGAATCCGCGACGAGTGGCAGTGGGGCTGGAAAAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTATGCAAATGGCCGTTAAGGATGTGTTGTTGAACGTGGt                                                             >  1:334216/1‑251 (MQ=255)
                                                                                                                                                                         cAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGaa                                                                                                                                                    <  1:26416/163‑1 (MQ=255)
                                                                                                                                                                          aGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTcgg                                                          >  1:92991/1‑251 (MQ=255)
                                                                                                                                                                                    cGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGCTCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGa                                                  >  1:49711/1‑249 (MQ=255)
                                                                                                                                                                                      gCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGCTGATGGCGACTCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACCGTGGTCTGCCAATGGCCGATCCGGCTGTGTTTATGAACGCTGGCGGCGGCGTGACGGt                                              >  2:17055/1‑251 (MQ=255)
                                                                                                                                                                                       cAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGt                                              >  1:341481/1‑250 (MQ=255)
                                                                                                                                                                                                  aGTCGGGGAGATGAAGTGACGTCCGGCAGGCGCGTGATGGTGGTGTGGGAAGGAACGCGTTCGTTGCTGGTGGAGAGTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAg                                  <  2:137365/251‑1 (MQ=255)
                                                                                                                                                                                                        ggagATGAAGTGACCTCCGGCAGCCTCGTGATGGGGGGGGGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAAGCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGa                            <  2:204339/248‑1 (MQ=255)
                                                                                                                                                                                                           gATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACtt                         <  1:45882/251‑1 (MQ=255)
                                                                                                                                                                                                            aTGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTa                        <  1:389398/251‑1 (MQ=255)
                                                                                                                                                                                                                aGTGACCTCCGGCAGCGCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGt                    <  1:396544/251‑1 (MQ=255)
                                                                                                                                                                                                                     ccTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTActg               <  1:414251/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cgTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGCCTTAGCGTTACTGCTGGCGATGGttt  >  1:326597/2‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                
ATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTT  >  NC_000913/4626582‑4627058

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