Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,626,673 G→A L254L (TTG→TTA radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,626,6730GA98.3% 172.4 / NA 61L254L (TTG→TTAradADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (22/37);  total (23/38)
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.98e-01

GATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCT  >  NC_000913/4626430‑4626909
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             
gATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGcc                                                                                                                                                                                                                                        <  1:380788/251‑1 (MQ=255)
          gACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAAcc                                                                                                                                                                                                                               <  2:305367/250‑1 (MQ=255)
                      gTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGt                                                                                                                                                                                                                  <  1:298075/251‑1 (MQ=255)
                      gTGATGCATATGGCGGATGTACAGACATCGCCTGGCAGCGGGGAGCAGGTGCGGGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACCCGTAACCAAAGATGGCCCGCTGGCGGGCCCTTAAGTGCTGGGACACTGGTACGCATGTTCGGTGCCTTTGGAAGGGCGATCAGCCTCCCGTTTTCCCAACTTACGCCGCCATAAAAACCCCttaggggaggg                                                                                                                                                                                                                  >  2:396544/1‑242 (MQ=255)
                         aTGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGaa                                                                                                                                                                                                               <  1:35723/251‑1 (MQ=255)
                           gCATAGGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATg                                                                                                                                                                                                             <  2:372209/251‑1 (MQ=255)
                             atatGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACCCGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGGATCGACTGTTCGGGGCTATTGAATGGCGTTGCCGACTCCCGGTTTCGCACTTTACGCAGCCATAAAAACCGCTTCAGGGCGGTGAATg                                                                                                                                                                                                             >  1:457861/1‑249 (MQ=255)
                               atGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGGGTGGCGATTGTCATGGTTGGAAACGTAACCAAAGATGGCTCGCTGGGGGGCCCTAAAGAGGGGGAGCACTGTTTCGACTGGTCGGGGCGTTTGGATGGGGAGGCCGGCTCCCGTTTTTGCAATTTACGCCGCCCTAAAAACCCCTTCgggccgggggatgggat                                                                                                                                                                                                         >  2:292824/1‑236 (MQ=255)
                                   cGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGGGCTGGAACCCTGGATCGACTGTTCGGTGCTTTTGGGGGGCGGTGCCGACCCCCGTTTTCGCACTTTACGCAGCCATAAAAACCggctcgggcgggggaatgggctgggg                                                                                                                                                                                                     >  2:45882/1‑226 (MQ=255)
                                          accgTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTAACTGACTCGCTTCTCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGTTCGCTGGCGGGCCCTAAAGTTCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATGAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTcg                                                                                                                                                                                              <  2:176099/248‑1 (MQ=255)
                                                     cTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGa                                                                                                                                                                                   >  1:51192/1‑250 (MQ=255)
                                                            gggggCGCTGGGGCGTGAAACGGTGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTTGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGc                                                                                                                                                                           <  2:16368/248‑1 (MQ=255)
                                                              tGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATAGTCATGGATGGACACGTAACCAAAGATGGGTCGCTGGCGGGCCCTATAGTGCTGGACCACTGTATCGACTGGTCGGGGCTTTTGGATGGCGATGCCTACTCCCCGTTTCTCACTTTTCGCAGCCATAAAAAACGCTTCGGCGCGGTGGATGAGACGTGCGACTTCGCCTAGGACCGAACAgggctg                                                                                                                                                                          >  2:293936/1‑248 (MQ=255)
                                                                  gcAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTg                                                                                                                                                                     >  1:137365/1‑251 (MQ=255)
                                                                     ggTGCATGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGATTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAg                                                                                                                                                                  <  1:9771/251‑1 (MQ=255)
                                                                         cGTGAAACGGCGGCTTATCGTACGCGCTTCGCCAAATCGCGCGGTGTGTCTATTTTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTGCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCTGCCATAACAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAg                                                                                                                                                              <  1:38551/251‑1 (MQ=255)
                                                                           tGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGACCACTGTATCGCCTGTTCGGTGCTTTTGGGGGGCGATGCCGACTCCCGTTTTCGCCCTTTACGCAGCCATAAAAACCGCTTCGGCGCGGGGGCTGAGCTGGGCGTCTTCGCGATGCCCGAACAGGGGCTGCGTggaggcaac                                                                                                                                                             >  2:359403/1‑242 (MQ=255)
                                                                           tGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCa                                                                                                                                                               >  1:197324/1‑248 (MQ=255)
                                                                           tGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCa                                                                                                                                                               <  2:197324/248‑1 (MQ=255)
                                                                                 ggcggcTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGTAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCtt                                                                                                                                                      <  1:175006/251‑1 (MQ=255)
                                                                                  gtggcttacctGACGCGCCTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTTGACACGTAACCAAAGATGGCGCGCTGGCGTGCCCTAAAGTGCTGGAACACTGGATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTTAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTc                                                                                                                                                     <  2:293833/241‑1 (MQ=255)
                                                                                  gcggcTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTc                                                                                                                                                     <  1:200037/251‑1 (MQ=255)
                                                                                    ggtttagctGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTGGCGGGCGGGCCCTAAAGTGCTGGAACACGGTATCGACTGTGCGGTGCTTTTTGATGGCGATGCCGACTCCCGTTGTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTGCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCgg                                                                                                                                                   <  2:162350/243‑1 (MQ=255)
                                                                                             tacgcgattCGCCAAAACTCGCGGTGGGGCGATTGTCATGGGTGGATACGGAACCAAAGATGGTCCGCTGGTGGGCCGTAAAGTGCTGGACCACTGTATCGACTGTACGTTGTTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGTCCGAACAGGGTCTGCGTGAATACAGCAATATTTCGGCAATTTTc                                                                                                                                           <  2:394924/243‑1 (MQ=255)
                                                                                             gacgcGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAAtt                                                                                                                                              <  1:241998/243‑1 (MQ=255)
                                                                                                gcgcgtCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGGTGGGGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCGTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTaa                                                                                                                                       <  1:61136/246‑1 (MQ=255)
                                                                                                 cGCTTCGCAAAAACGCGGGGTGTGGCGATTGTCATGGTGGGGCGCGTACCCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACGCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAg                                                                                                                                      <  2:114656/251‑1 (MQ=255)
                                                                                                   cTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGCCTGCTCGGTGCCTTTGGATGGCGCTGCCGACTACCGTTTTCGCCATTTACGCACACATAAAAACCACTTCGGCGCGGGTAATGAGCAGGACTGCTTCGCGATAGCCCACCCGGGGCTGCGTGATGTAAGACACACTTAGGCAATTTTCTTTAGTc                                                                                                                                    >  2:405223/1‑251 (MQ=255)
                                                                                                     tCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTc                                                                                                                                    <  1:144475/249‑1 (MQ=255)
                                                                                                        ccAAAACGCGCGGTGTGGCGATTGTAATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCgcgga                                                                                                                               <  2:331747/251‑2 (MQ=255)
                                                                                                          aaaaCGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGCTGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCCCGATGGCCGAACAGGGGCTCCGTGAAGTCAGCAACCCTTCGGACATTTTTTTAAGTCgcgggg                                                                                                                              >  2:389398/1‑248 (MQ=255)
                                                                                                               ggggggTGTGGGGATTGTCATGGTTGGACACGGAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAg                                                                                                                        <  2:9683/247‑1 (MQ=255)
                                                                                                               gggcgGTGTGGCGATTGTCAGGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCTCTAAAGTGCTGGAACACTGTATCGACGGTTCGGTGCTTTTGGAGGGCGATGCCGACTCCCGTTTGCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAg                                                                                                                        <  2:469532/249‑1 (MQ=255)
                                                                                                                cgcgGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAgg                                                                                                                       >  1:413023/1‑250 (MQ=255)
                                                                                                                cgcgGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGAGGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGGACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGt                                                                                                                       <  1:32672/251‑1 (MQ=255)
                                                                                                                   ggTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGAc                                                                                                                    <  1:221702/251‑1 (MQ=255)
                                                                                                                            aTTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGGGCTGGAACCCTGGATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCCGCCATAAAAAACGCTTCGGCGCGGTGAATGGGCTGGGCGTCTTCGCGAAGACCGAAAAGGGGCTGCGGGACGTCAGCAACCCTTCGGCAATTTTCTTAAGGCgcggagtgggagggccctccggccg                                                                                                           >  2:414251/1‑230 (MQ=255)
                                                                                                                              tGTCATGGTTGGACACGTAACCAACGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTGTCTTAATTCGCGGAGATGAAGTGACCTCCGGCAGc                                                                                                          <  1:483537/250‑1 (MQ=255)
                                                                                                                              tGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGa                                                                                                                     >  1:102891/1‑239 (MQ=255)
                                                                                                                              tGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGa                                                                                                                     <  2:102891/239‑1 (MQ=255)
                                                                                                                                    ggttggacACGTAACCAAAGATGGCTCGCTGGCGGGCCCTCAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGGGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGGGATGAAGTGCCCTCCGGCAGCTCCGTg                                                                                                    >  2:473478/8‑250 (MQ=255)
                                                                                                                                    ggttggacACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTAGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGCCTTCGCGATGACCGAACCGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCCTAAGTCGCGGAAGTGAAGTGACCTCCGGCAGCTccgtt                                                                                                    >  2:404505/8‑246 (MQ=255)
                                                                                                                                        ggacACGTAACCAAAGATGGCTCGGTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTGAATTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGTtggt                                                                                               <  1:130828/248‑1 (MQ=255)
                                                                                                                                                        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)
                                                                                                                                                                                                          cGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGAtggtgg                                                                                             <  1:29904/189‑1 (MQ=255)
                                                                                                                                                                                                          cGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGAtggtgg                                                                                             >  2:29904/1‑189 (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)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             
GATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCT  >  NC_000913/4626430‑4626909

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