Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,626,601 G→T P230P (CCG→CCT)  radA → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,626,6010GT93.9% 199.9 / ‑2.1 68P230P (CCG→CCT) radADNA recombination protein
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base T (26/36);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGG‑CACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGG  >  NC_000913/4626356‑4626845
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
gACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGt                                                                                                                                                                                                                                                 <  1:66870/251‑1 (MQ=255)
gACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGG‑CACGTAACCAAAGATGGTTCGCTGGCTGGCCCGaaa                                                                                                                                                                                                                                                   <  1:141425/249‑1 (MQ=255)
 aCAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGG‑CACGTAACCAAAACTGGCTCGCTGGCGGGCCCCAAAggg                                                                                                                                                                                                                                                >  1:81246/1‑249 (MQ=255)
        cAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGa                                                                                                                                                                                                                                           <  1:411433/249‑1 (MQ=255)
        cAATATGTAGTCGGAAACCAGCATCGAGCAGAGCTGCCTGATTGCCGAAGAAGAGCATCCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTTCAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGTCGATTGTCATGGTTGGA‑CTCGTAACCAAAGATGGCTCGCTGGCGGGTCCTAAAGTGCTg                                                                                                                                                                                                                                             <  1:56120/247‑1 (MQ=255)
        cAATATGGTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTg                                                                                                                                                                                                                                             <  1:418521/247‑1 (MQ=255)
            atGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGCTCCAGGTGATGCATATGGCGGATGTACAGTCCTCGCCTTGCAGCGTGGCGCAGGTGCGTGAAACTGCGGCTTATATGACGCGCTTCGCCCAAACCCGCGGTGTGGCGCTTTTCATTGTTTGG‑CCCCTAACCCAAGATGGCTCCCTGGCGGGCCCTCAAGTGTTGGTACACtt                                                                                                                                                                                                                                     >  2:38037/1‑250 (MQ=255)
                tgtCGGAAACCGGCATCGATCAGATCTGCCTCGTTACCGAAGAAGAGCAACCGAAGCTGATGGTAAGTGACTCTGTCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAATGTGCTGGAACACTGTAt                                                                                                                                                                                                                                  <  2:217997/250‑1 (MQ=255)
                              aTCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGAAGCGCTTCGCCCAAACGCGCGGTGTGGCGGTTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGTCGGGCCCTAAAGGGCTTGAAAACTTTATAGGCTGTTCGGTGCt                                                                                                                                                                                                                   >  1:204899/1‑251 (MQ=255)
                                   gcAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGGGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTgg                                                                                                                                                                                                              <  1:284157/250‑1 (MQ=255)
                                           gccggATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATTCATATGGCGGAGGTACATTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCg                                                                                                                                                                                                        <  2:283149/245‑1 (MQ=255)
                                                    ccGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTcc                                                                                                                                                                                             <  1:206336/251‑1 (MQ=255)
                                                            gagCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCg                                                                                                                                                                                     <  1:38037/251‑1 (MQ=255)
                                                            gagCAACCGAAGCTGATGGGAATTGACTCGAGCCAGGTGATGCATATGGCGGATGTACAGTCAGCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGTCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCGCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGGTCGGTGCGTTTGGATGGCGATGCCGACTCCCGTTTTCg                                                                                                                                                                                     <  2:75371/251‑1 (MQ=255)
                                                                 aCCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGGTTTTGGATGGCGATGCCGACTCCCGTTTTAGCActt                                                                                                                                                                                >  1:304387/1‑249 (MQ=255)
                                                                     aagtggagggTAAGTGACTCGATACAGGTGATGCATATGGCGGAGTTACAGTCAGCGACGGGCAGCGTGGCGGAGGTGCGTGAAGAGGCGGCTTAGCTGACGCGCTTCGCCAAAACGCGAGGTGTGGCGATGGTCATGGGTGGA‑CACGTAACCAAAGATGGCTCGCGGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCTGTTTTCGCAcgttacc                                                                                                                                                                            <  2:101316/243‑7 (MQ=255)
                                                                       gCTGATGGTAATTGACTAGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTTTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTTGAACACTGTATCGCCTGTTCGGGGCTTTTTGGTGGCGATGCCGACTCCCGTTTTCGCActtttcgc                                                                                                                                                                           >  1:356192/1‑243 (MQ=255)
                                                                             gtgaaTTGACTCGTTCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGTA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACTCTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGTCTCCCGTTTTCGCACTTTACGCAGCCAt                                                                                                                                                                     <  1:278031/247‑1 (MQ=255)
                                                                                    gacccGATCCAGGTGAGGCATATGGCGGATGTACCTGCACCGCCTGGCAGCGTGGCGCAGGGGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAAcc                                                                                                                                                              <  2:12537/246‑1 (MQ=255)
                                                                                       ccGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCtt                                                                                                                                                          <  1:415845/250‑1 (MQ=255)
                                                                                        cGATCCAGGTGATGCATAGGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTc                                                                                                                                                         <  2:86818/251‑1 (MQ=255)
                                                                                            ccAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCGGACGCGCTTCTCCAAAACGCGCGGTGTGTCGATTGTCATGGTTGGA‑CCCGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGgcg                                                                                                                                                     <  1:281096/251‑1 (MQ=255)
                                                                                            ccAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGGGAAACGGCGGCTTATCTGACGCGCTTCGCCCAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGGGCTGGGACACTTTATAGACTGTTTGGTGCTTTTTGGTGGGCAGGCCGACCCCCGTTTTCGCACCTTTCGCAGCCATAAAAACCCCTTCgggg                                                                                                                                                     >  2:263853/1‑249 (MQ=255)
                                                                                                gTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCCAAAGTGCTTGAACACTGGATCGACTGTTCGGTGCTTTTGGGTGGCGATGCCGCCCCCCGTTTTCGGCCTTTACGCAGCCATAAAAACCGCTTCGGCGCggg                                                                                                                                                 >  1:70731/1‑250 (MQ=255)
                                                                                                              ggATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCg                                                                                                                                   >  1:395397/1‑251 (MQ=255)
                                                                                                              ggATGTACAGGCATCGCCTGGCCGCGTGGCGCAGGTGCGTGAAACGGCTGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATGGTCATGGTTGGA‑CACGTAACCAAAGATGGCGCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATTCCGACTCCCGTTTTCTCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCg                                                                                                                                   <  2:414208/251‑1 (MQ=255)
                                                                                                                         cATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATg                                                                                                                        >  1:501662/1‑251 (MQ=255)
                                                                                                                            cGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGGGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCTACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGAcc                                                                                                                     >  1:176752/1‑251 (MQ=255)
                                                                                                                                  gTAGCGTGGCGCAGGTACGTGAAACGGCGGCTTATTTGACACGCTTCGCCACAACGCGCGGTGTGGCGATTGTCATGGTGGGG‑CACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGAGGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAAc                                                                                                                 >  1:242143/1‑249 (MQ=255)
                                                                                                                                   caGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGAt                                                                                                                                                                                                      >  1:127043/2‑163 (MQ=255)
                                                                                                                                   caGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGAt                                                                                                                                                                                                      <  2:127043/162‑1 (MQ=255)
                                                                                                                                    aGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAggg                                                                                                             <  1:192747/251‑1 (MQ=255)
                                                                                                                                     gCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAgggg                                                                                                            >  1:227028/1‑251 (MQ=255)
                                                                                                                                          gcgcAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCg                                                                                                       <  1:73599/251‑1 (MQ=255)
                                                                                                                                             cAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTg                                                                                                     >  1:357240/1‑250 (MQ=255)
                                                                                                                                                  gCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTc                                                                                                <  1:14916/250‑1 (MQ=255)
                                                                                                                                                     tGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGGAGGCAGc                                                                                             >  1:353413/1‑250 (MQ=255)
                                                                                                                                                     tGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGc                                                                                             >  1:396639/1‑250 (MQ=255)
                                                                                                                                                         aCGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTTGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAAccc                                                                                        >  1:128835/1‑251 (MQ=255)
                                                                                                                                                         aCGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTGTACGCAGCCATAAAATCCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAAccc                                                                                        <  1:485927/251‑1 (MQ=255)
                                                                                                                                                             gggctgttttGACCCGCTTCGCAAAATCGGGCGTGGTGGGTATTGTCA‑GGTGGGGGCACGTAACCAAAGATGGTCCGCGGGCTGGCCCGAAAGTGCTGGATCCCTGTATCGACTTTTCGGTGCTTTTGGATGGCGATGCCGACCCCCGTTTTCGCCCCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTAGGCGTCTTCGCGATGACCGACCAGGGGCTTCGTGAAGTCAGCAACCCTTCg                                                                                    <  2:242143/244‑1 (MQ=255)
                                                                                                                                                                    actGACGCGCTCCGCCAAAACTCGCGGTGTGGCGATTGTCATGGTTGGA‑CACTTATCCAAAGATGGCTCGCTGGCGGGCCCTAAAGGTCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGGGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCATCGGCAAtt                                                                              <  2:219252/248‑1 (MQ=255)
                                                                                                                                                                         cgcGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGGGCTTTTGGATTGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGACTGAACTGGGCGTCTTCGCGATGACCGAACAGGGGGTGCGGGAAGTCAGCCAACCTTCGGCAATTTTCtt                                                                         >  2:424515/3‑250 (MQ=255)
                                                                                                                                                                         cgcGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTAGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCCTAAAAACCCATTCGGCGCGGGGAATGAGCTGGGCGTCTTCGCGATGGCCGACCAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCGATTTTCtt                                                                         >  2:22313/3‑250 (MQ=255)
                                                                                                                                                                            gCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGt                                                                     <  1:870/251‑1 (MQ=255)
                                                                                                                                                                               tCGCAAAAACGCGCGGTGAGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTcg                                                                   <  2:185869/250‑1 (MQ=255)
                                                                                                                                                                                cGCCAAAACGCGCTGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTcgcg                                                                 <  1:482654/251‑1 (MQ=255)
                                                                                                                                                                                    aaaaCGCGCGGTGTGGCGATTGTCATGGGTGGA‑CACGTAACCAAAGATGGCTCGCTGGAGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAAACATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCCCGATGAACGACCAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTACGTcgaggagg                                                             >  1:202048/1‑245 (MQ=255)
                                                                                                                                                                                     aaaCGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGAt                                                            <  1:37313/251‑1 (MQ=255)
                                                                                                                                                                                      aaCGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGTGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATg                                                           <  1:220887/251‑1 (MQ=255)
                                                                                                                                                                                                  ggCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTTGGCGTCTTCGCGGTGACCGCACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGGGATGAAGTGACCTCCg                                               >  2:146629/1‑251 (MQ=255)
                                                                                                                                                                                                  ggCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTCCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGtcctcc                                                <  1:115357/250‑1 (MQ=255)
                                                                                                                                                                                                          tCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGACGTGACCTCCGGCAGCTcc                                       >  1:249461/1‑250 (MQ=255)
                                                                                                                                                                                                          tCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTcc                                       >  1:373548/1‑250 (MQ=255)
                                                                                                                                                                                                            aTGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGACGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGCTGAAGTGCCCTCCGGCAGCTcggg                                     >  2:8641/1‑248 (MQ=255)
                                                                                                                                                                                                               ggttggacACGTCACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGATGGAACCCTGTAACGAATGTTCGGTGCTTTTGGATGGCGATGCCGACTACCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGAACGTACAGGGGCTGCGAGAAGTCAGCAACCCTTCGGCCATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCggccgctccgtgg                                   >  2:240669/8‑240 (MQ=255)
                                                                                                                                                                                                                       aCGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCt                                                                          <  1:137280/204‑1 (MQ=255)
                                                                                                                                                                                                                       aCGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGGGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCt                                                                          >  2:137280/1‑204 (MQ=255)
                                                                                                                                                                                                                           aaCCAAAGATGGCTCTCTGGTGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTTCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGGCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTggg                       <  2:70731/251‑1 (MQ=255)
                                                                                                                                                                                                                           aaCCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGCTGGTGGTGTggg                       >  1:321089/1‑251 (MQ=255)
                                                                                                                                                                                                                            aCCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGgga                      <  2:396639/251‑1 (MQ=255)
                                                                                                                                                                                                                                aaGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGggaagga                  <  1:263853/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGAGGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAAc                <  2:124382/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ggctCGCTGGCGGGGCCTAAAGGGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACGTGACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGGGGTGTGGGAAGGAAcgcg             <  2:81246/248‑1 (MQ=255)
                                                                                                                                                                                                                                             ggcgggcCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTccgttgc     >  1:348109/5‑246 (MQ=255)
                                                                                                                                                                                                                                             ggcgggcCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTccgttg      <  1:326786/246‑5 (MQ=255)
                                                                                                                                                                                                                                              gcgggcCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTccgttgct    <  2:249461/248‑7 (MQ=255)
                                                                                                                                                                                                                                                gggcCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTCCGTGATGGGGGTGTGGGAAGGAAAGCGTCCGTTGCtgg  >  2:476498/2‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
GACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGG‑CACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGG  >  NC_000913/4626356‑4626845

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: ‑