Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,627,159 C→T I416I (ATC→ATT)  radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,627,1590CT96.2% 160.0 / NA 53I416I (ATC→ATT) radADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (23/28);  minor base G (0/2);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTAGCTGATGCCAGCGGTATGACCAA  >  NC_000913/4626925‑4627400
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  
ggaaattctgctgTCTGTCTGGCACCGGCACGGGGGGCTGCAAATGGCCGATCAGGATGTGTTGGTGAACGTGGGCGGCGGCGTGAAGGTAACCGAACCCAGCGCCGACTTAGCGTTACTGCTGGCGTTGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCACAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAac                                                                                                                                                                                                                                    <  2:272108/244‑1 (MQ=255)
   aagccTGCTGGCTTTCTTGCCCCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑g                                                                                                                                                                                                                                <  2:386556/248‑1 (MQ=255)
          cgggctttcttGCACCGTCTCGGTGTTCTGCAAATGGCCGATCTGGATGTGGTTGTTAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTGCCAGCCGGCGTGACATACCGCTGCCATAGGAGCTGGTGGTGTTTGGTGAAGTCTGGCTGGCAGGGGAGATCCGCCGGTTGCCCAGCGTTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCg                                                                                                                                                                                                                         <  1:384155/242‑1 (MQ=255)
                gtctttcccgggagtgggggtCTGCAAATGGCCGACCTGGATTTGTTTGTGAACGTGGTCGGGGGCGGGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCGGCGTTAGAGACCGCTGCCTCAGGATCTGGTGGTGTTTGGTGAAGCGGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACTAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCggg                                                                                                                                                                                                                    <  2:505363/232‑1 (MQ=255)
                      cACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGTGTGAAGGTAACCGAAACCAGCGCCGTCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGt                                                                                                                                                                                                             <  1:273977/251‑1 (MQ=255)
                        ccGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTc                                                                                                                                                                                                           >  1:609796/1‑251 (MQ=255)
                            cacggtgggctgcgagcggcggagcAGGATGTGGTGGTGATCGGGGTCGGTGGCGTGAGGGCAGCCGAAACCAGGGCCGACTTAGCGTTATTGCTGGCGTTGGTTTCCAGCCTGCGTGAGAGCCCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCTGTGCCCAGCGGTCATGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGTCGGGCGATTGTTCcgg                                                                                                                                                                                                        <  2:265148/226‑1 (MQ=255)
                            cACGGTAGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGAAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGCCAGGAACGAATTTCAGAAGCGGCGGAACAA‑GGATTTCGCCGGGCGTTTGTTCcggc                                                                                                                                                                                                       >  1:632105/1‑251 (MQ=255)
                            aaCGGGGGTCTGAAAATGGCCTACCAGGATGTGTTGGGGAAAGTGGTCGGTGGCGTGAAGGTAACCGAATCCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCcgg                                                                                                                                                                                                        <  2:107367/249‑1 (MQ=255)
                              cGGTGGTCGGCAAATGACCGATCAGGATGTGTTTGGGAACGTGGTCGGCGGCGTGAAGGAAACCGACACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGAGCACAGCGGTCAGGCACGAATTTCAGAAGCGGCGAAA‑ACAGGTTTTCGACGGGAGCTTGTTCcggcgg                                                                                                                                                                                                     >  2:525456/1‑251 (MQ=255)
                                       gaaaggggCCGATCGGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGGCTTAGCGTTACTGCTGGCGATGGTTGCCAGCCGGCGGGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCTGGCGATTGTTCCGGCGGCTAACGtac                                                                                                                                                                                            <  2:81565/245‑3 (MQ=255)
                                              gCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGAGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGCACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTACAGTACCGaaaa                                                                                                                                                                                      >  2:491896/1‑250 (MQ=255)
                                                    cAGGATGTGTTGTTGAAAGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGTCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGGGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAgcgc                                                                                                                                                                                <  1:528094/250‑1 (MQ=255)
                                                          gtgtTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAgg                                                                                                                                                                         <  1:236611/251‑1 (MQ=255)
                                                               tgtgAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATg                                                                                                                                                                     <  1:491896/250‑1 (MQ=255)
                                                                gtgAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAACAAGCGCCGGAAGGGATGc                                                                                                                                                                    >  2:232563/1‑250 (MQ=255)
                                                                                  gTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgcgg                                                                                                                                                                                                                                          >  1:463124/1‑163 (MQ=255)
                                                                                    gAAGGGAACCGAAACCAGCCCCGCCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGTCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAAc                                                                                                                                               <  2:83922/251‑1 (MQ=255)
                                                                                     aaGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAAct                                                                                                                                              >  1:149109/1‑251 (MQ=255)
                                                                                      aGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCt                                                                                                                                                                                                   >  1:315603/1‑197 (MQ=255)
                                                                                          aaCCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTcc                                                                                                                                          <  1:285628/250‑1 (MQ=255)
                                                                                          aaCCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCAGGCGGCTAACGTAACGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTCAAAAACTCTCCg                                                                                                                                         >  1:73171/1‑251 (MQ=255)
                                                                                                  ccAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCtt                                                                                                                                 <  1:523349/251‑1 (MQ=255)
                                                                                                  ccAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCtt                                                                                                                                 >  1:576170/1‑251 (MQ=255)
                                                                                                   cagcgCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTa                                                                                                                                <  1:446095/247‑1 (MQ=255)
                                                                                                    agcgCCGACTTAGCGTTACTGCTGGAGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGGTCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGGTTGTTCCGGCGGCTAACGTACCGAAAAAAGCGACGGAAGGGATGAAGCATTTTGGCGTTAAAAAGctcctcaccccgctt                                                                                                                                 >  2:163077/4‑237 (MQ=255)
                                                                                                          gACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGGGtt                                                                                                                         >  1:166709/1‑251 (MQ=255)
                                                                                                          gACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAA‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAGAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGAAGCGCTTAGCGTGtt                                                                                                                         >  1:521929/1‑251 (MQ=255)
                                                                                                           agttAGCGTTACTGCTGGCGAGGGTTTCCAGCCTGCGTGACAGACCGCGGCCACAGGTGCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCGCAGCGGTCAGGGACGAAGTGCAGAAGCGTCGAAACAC‑GGTTTTCGCCGGGCGATGGTTCCGGTGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTGTGGCGTTAAAAAACTCTCCGACGCGGGTAGTGTGTTc                                                                                                                        <  2:79731/249‑1 (MQ=255)
                                                                                                                ggtgtACTGCTGGCGGGGGTTTCCAGCCTGCGTGACAGATCGCTGCCACAGTATCTGGTGGTGTTTGGTGAAGGCGGGCTGGCAGGGGAGATCCTCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTGTGCGTGTTcgacg                                                                                                                    <  2:290367/246‑1 (MQ=255)
                                                                                                                   ttACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACtt                                                                                                                >  1:303629/1‑251 (MQ=255)
                                                                                                                        gctgggtttGGTTTCCAGCCTTCGTGAAAGATCTCGGCCACAGGATCTGGGGGTGTTTGGGGAAGTCGGGCTGGCAGTGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATGTCAGAAGCGGCGAAACAC‑GGGTTTCGCCTGGCGATTGTTCCGGTGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGTCGCGCTTAGCGTGTTCGACGATTTATAAt                                                                                                           <  2:539755/243‑1 (MQ=255)
                                                                                                                        gctgGCGATGGTGTCCAGCCGGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCTGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATaa                                                                                                            <  1:609692/250‑1 (MQ=255)
                                                                                                                         ctgGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATg                                                                                                          >  1:636584/1‑251 (MQ=255)
                                                                                                                            gCGATGGTTTCAAGCCTGCTGGACAGACCGCTGCCACAGGAGCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTAGAATgagg                                                                                                       <  2:609796/251‑2 (MQ=255)
                                                                                                                                     tCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACgga                                                                                               <  1:525456/250‑1 (MQ=255)
                                                                                                                                                  aCAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGGTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATgtcgt                                                                                 <  2:166709/251‑1 (MQ=255)
                                                                                                                                                              cacAGGATCTGGTGGTGTTTGGAGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGGCTTATAATGAGGAATACGGGGGGAGATATGTCGTCTTTTGATTAcc                                                                     >  2:108868/3‑251 (MQ=255)
                                                                                                                                                               acAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATgag                                                                                                        <  1:159622/214‑1 (MQ=255)
                                                                                                                                                                        ggtggtGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAAc                                                              <  1:564928/248‑1 (MQ=255)
                                                                                                                                                                         gtggtgTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGcc                                                          <  1:123332/251‑1 (MQ=255)
                                                                                                                                                                              gtTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCa                                                      <  1:33711/250‑1 (MQ=255)
                                                                                                                                                                              gtTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCa                                                      >  1:673933/1‑250 (MQ=255)
                                                                                                                                                                                 tGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTCGCGTGTTCGACGACtttt                                                                                                              >  2:23174/1‑189 (MQ=255)
                                                                                                                                                                                 tGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTtat                                                                                                              <  1:23174/191‑1 (MQ=255)
                                                                                                                                                                                           gggCTGGCAGTGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGc                                        >  1:541912/1‑251 (MQ=255)
                                                                                                                                                                                                  cAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTAc                                 >  1:415020/1‑251 (MQ=255)
                                                                                                                                                                                                     gggAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTAcagc                              >  1:622946/1‑251 (MQ=255)
                                                                                                                                                                                                                         aGCGGTCAGGAACGAATGTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGt          <  1:499498/251‑1 (MQ=255)
                                                                                                                                                                                                                            ggTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACa                                               <  2:34286/211‑1 (MQ=255)
                                                                                                                                                                                                                            ggTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACa                                               >  1:34286/1‑211 (MQ=255)
                                                                                                                                                                                                                             gTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATg       >  1:227851/1‑250 (MQ=255)
                                                                                                                                                                                                                                 ggAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCaa  <  1:465289/251‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  
GGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTAGCTGATGCCAGCGGTATGACCAA  >  NC_000913/4626925‑4627400

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