Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,626,042 T→C M44T (ATG→ACG)  radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,626,0420TC94.8% 162.9 / ‑0.3 59M44T (ATG→ACG) radADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (36/19);  total (39/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.14e-01

ATCCGTCACGCTGACCTGATGGGGGTATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAA  >  NC_000913/4625801‑4626282
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
aTCCGTCACGCTGACCTGATGGGGGTATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGcgcg                                                                                                                                                                                                                                         >  2:128097/1‑251 (MQ=255)
      cACGCTGACCTGATGGGGGTATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGGGGCGCGTAACg                                                                                                                                                                                                                                    >  2:386978/1‑250 (MQ=255)
               ctgctGGGGGGATTCTGCATCATCTCAGACAGCCTGCATCAGAAGTAATTGCTAACCCGCCATACTGCGGGCGGCACAGCATTAACGAGGTACACCGTGTCAGAAGCGCCAAACCGCGCCTTTGTTTGTAATGAATGCGGAGCCGAATATCCACGCTGCGAGGGGCCGTGCAATGCCTGACATGACTGGCATACAATCCCCGAGTGGCCACTTGTTGCGTAGCCAACGGTGGCGCGTAACAAGCGtctcc                                                                                                                                                                                                                           >  2:480184/5‑249 (MQ=255)
                  aTGGGGGTATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCgg                                                                                                                                                                                                                       >  1:61006/1‑251 (MQ=255)
                     ggggTATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTa                                                                                                                                                                                                                    <  1:491361/251‑1 (MQ=255)
                       ggTATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATg                                                                                                                                                                                                                  >  1:408585/1‑251 (MQ=255)
                         tATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGcc                                                                                                                                                                                                                >  1:224043/1‑251 (MQ=255)
                          aTTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGGGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCCATGCCg                                                                                                                                                                                                               >  2:373779/1‑251 (MQ=255)
                                       tcAGGCAGCCTGAATCAGAAGTAATTGCTCTCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCATGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGACTATGCCGGTAGCGCCGGGGt                                                                                                                                                                                                  >  1:399877/1‑251 (MQ=255)
                                       tcAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCggggg                                                                                                                                                                                                  >  2:8155/1‑250 (MQ=255)
                                            cAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGCGGTGCGTCTTGCTGCGTCGCCAACGGCGGCGCGTAACGAGAGGCTCAGCGGCCATGCCGGTCGCGCCGGGGGGGCa                                                                                                                                                                                              >  2:470635/1‑250 (MQ=255)
                                                    aTCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGt                                                                                                                                                                                                                               <  1:264111/209‑1 (MQ=255)
                                                    aTCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGAGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGt                                                                                                                                                                                                                               >  2:264111/1‑209 (MQ=255)
                                                     tCAGAAGTAATTGCTCGACCGCCATCCTGCGGGCGGCACAGCCTTACCGAGGTAAACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGAGGCAGTGCAGTGACTGTCATGACTGGAACACCATCCCCGAGGTGCGTCTTGCTGCGTCGCCCACGGTGGCGCGTAACGAGCGGCTCAACGGCTATGCCGGTATCGCCGGGGTGGCGAAAgccc                                                                                                                                                                                       >  2:483582/1‑245 (MQ=255)
                                                         aaGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAAct                                                                                                                                                                                <  1:257626/251‑1 (MQ=255)
                                                          aGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGt                                                                                                                                                                                                  <  1:209803/232‑1 (MQ=255)
                                                          aGTAATTGCTCCCCCGCAATCCTGCGGGCGGCACAGCATTAACGAGGTACCCCGTGGCAAACGCTCCAAAACGCGCCTTTGTTTGTAATGACTGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCCACGGTGGCGCGTAACGAGCGTCCCCGCGGCTATGCCGGTAGCGCCGGGGtaggtcgggagggcgtcgtg                                                                                                                                                                               >  2:209803/1‑232 (MQ=255)
                                                             aaTTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCCCAAGTACAGAAACTCTcc                                                                                                                                                                            >  2:221977/1‑251 (MQ=255)
                                                                  cTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGatat                                                                                                                                                                       >  1:360676/1‑251 (MQ=255)
                                                                    cgcccgccATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATc                                                                                                                                                                      >  1:145539/1‑250 (MQ=255)
                                                                      cccgccATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCACCAAACCGCGCTGTGGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCTGGGTGGCAAAAGTCCAGAAACGCTCCGAGATCAGc                                                                                                                                                                   <  2:245572/251‑1 (MQ=255)
                                                                          ccATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCACAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACCCCATCCCCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGACTCAGCGGCTCTGCCGGTTGCGCCGGGGTGCCAGAAGTCCCGAAACTCACCGATATCAGCCttt                                                                                                                                                               >  2:115357/1‑250 (MQ=255)
                                                                            aTCCTGCGGGCGGCACACCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCCACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTACGATATCAGCCTTgc                                                                                                                                                              >  1:101316/1‑249 (MQ=255)
                                                                               cTGCGGGCGGAACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAGACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCCGGCAGGGCCAATGCAGCGCCTGTCATGCCTGGAACACCATCACCGACGTGCTTCTTTCTGCGTCGCCAAAGGTGGCGCGTAAAGAGCGGCGCAACGGATAGGCCGGTGGCGCCGGGGGTGCAAAAGTCCAGCAACACACCGATATCATCATTGAagcg                                                                                                                                                          >  2:411433/1‑249 (MQ=255)
                                                                                  gggtcggcACAGCATTAACGAGGTACACCGTGGTAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGTGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCTCCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTg                                                                                                                                                       <  1:208338/247‑1 (MQ=255)
                                                                                         cacaGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGtt                                                                                                                                                <  1:126661/251‑1 (MQ=255)
                                                                                                           cacCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTaa                                                                                                                                >  1:75371/1‑249 (MQ=255)
                                                                                                            acCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAaga                                                                                                                             <  1:142048/251‑1 (MQ=255)
                                                                                                                  gCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCg                                                                                                                        >  1:54653/1‑250 (MQ=255)
                                                                                                                    aaaaGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGATATGCCGGGAGAGCCGGGGGGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGGTTTAAAGATTTCGAcc                                                                                                                     >  2:418521/1‑251 (MQ=255)
                                                                                                                        gCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGt                                                                                                                 <  2:145539/251‑1 (MQ=255)
                                                                                                                         cTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGCGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGACCAGAAACTCTCCGATATAAGCCTTGAATAACTGCCCCGTTTTTCCACCGGGTTTAAAGAGTTCGACagcgg                                                                                                                 >  2:284157/1‑249 (MQ=255)
                                                                                                                           ccAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGGCCAGCACCCCTCCGAGATCCGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGCTTTAAAGAGTTCGAcctcgcact                                                                                                              >  2:104495/1‑244 (MQ=255)
                                                                                                                            cAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTa                                                                                                             <  1:483582/251‑1 (MQ=255)
                                                                                                                             aaaaCGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAg                                                                                                            <  1:267524/251‑1 (MQ=255)
                                                                                                                                aCGCGCCTTGGTTTGTAATGAATGCGGGGCCGATTACCCGCGCAGGCAGGGGCAGTGCAGGGCCTGTCATGCCTGTAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTTGCCGATATCAGCCTCGAAGAGCGGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACAAggcg                                                                                                         <  2:22900/251‑1 (MQ=255)
                                                                                                                                 cgcgCCTTTGTTTGTAATGAATGCGGGGCCGATTAGCCGCGCTGGCAGGGGCAGTGCAGTGCCGGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcgg                                                                                                        <  2:103094/251‑1 (MQ=255)
                                                                                                                                   cgcCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcg                                                                                                      >  1:123644/1‑251 (MQ=255)
                                                                                                                                         tgtttgtAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGCTCGACCGCGTACTAGGCGGCGGCgtgg                                                                                                >  2:339437/1‑251 (MQ=255)
                                                                                                                                            ttgtAAAGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCAGCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATAGCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCgtggtg                                                                                              <  2:408585/250‑1 (MQ=255)
                                                                                                                                               tAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGGGCGTCTTGCTGCGTCGCCAACGGTGGAGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCAGCGGTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGAGGCGGCGTGGTGcc                                                                                            >  2:447006/1‑249 (MQ=255)
                                                                                                                                                   gaatgCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGa                                                                                        >  1:86818/1‑249 (MQ=255)
                                                                                                                                                    aatgCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCCTCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGGCTCAGCGGCTATGACGGTAGCGCCGGGGTGGCCAAAGTCCAGAACCTCTCCGATATCAGCCTTGAAGTGCTGCCGCGTTTTTCCACCGGTTTTAAAGAGTTCGACCGCGTCCTAGGCGGCGGCGTGGTGCCAGGag                                                                                       >  2:230845/1‑248 (MQ=255)
                                                                                                                                                      tgCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGCCCGCGTACTCGGCGGCGGCGTGGTGCCAGGAAGTGc                                                                                   >  1:164512/1‑251 (MQ=255)
                                                                                                                                                                      cgcgCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGt                                                                   >  1:418821/1‑251 (MQ=255)
                                                                                                                                                                               aGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAAACCTGGt                                                          >  1:240827/1‑251 (MQ=255)
                                                                                                                                                                                   gcagtgcagtGCCTGTCATGACTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCGATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCgg                                                      <  2:372213/251‑1 (MQ=255)
                                                                                                                                                                                     agtgcagtgcCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGCACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGAAACCCTggggggggg                                                    >  2:66870/1‑244 (MQ=255)
                                                                                                                                                                                        gcagtgcCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGGGCGGGGaa                                                  >  1:185869/1‑250 (MQ=255)
                                                                                                                                                                                                   cATGCCTGGACCCCCCTCGCCGCGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCTGCTATGCCGGTAGCGCCGGGGTGGCAAAAGGCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgctgt                                       <  2:176259/250‑2 (MQ=255)
                                                                                                                                                                                                   cATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgctgtt                                      >  2:393101/1‑249 (MQ=255)
                                                                                                                                                                                                                ccATCCCCGAGTTGCGTCTTGCTGCGTCTCCAACGGTGGCGCTTACCGTGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACGCTCCGAGATCCGCCGTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTGCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCTGGGAAATCCACGCTGTTGCTGCCAACGCtg                         <  2:369179/251‑1 (MQ=255)
                                                                                                                                                                                                                      ccGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCaaa                   <  1:373779/251‑1 (MQ=255)
                                                                                                                                                                                                                            tGCGTCTTGCTGCGTCGCCAACGGTCGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGcc                                                                                  >  1:141681/1‑182 (MQ=255)
                                                                                                                                                                                                                                   tgctgcGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGAGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCTGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCcagc         >  1:414208/1‑248 (MQ=255)
                                                                                                                                                                                                                                      tgcGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCcagcagctg    >  2:4909/1‑247 (MQ=255)
                                                                                                                                                                                                                                       gcGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCAGCAATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACGGGCCCAGCAGATGa   >  2:141425/1‑250 (MQ=255)
                                                                                                                                                                                                                                       gcGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTTTGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTTCCGCGTTTTTCCACCTTATTTAAAGATTTCGACCGCGTACTAGGCTGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGaa  <  1:225868/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gcGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTGGCCCAGCAGATGa   <  1:8155/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
ATCCGTCACGCTGACCTGATGGGGGTATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAA  >  NC_000913/4625801‑4626282

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