Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,090,968 G→A E208K (GAG→AAG)  endA → DNA‑specific endonuclease I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,090,9680GA98.4% 219.5 / NA 64E208K (GAG→AAG) endADNA‑specific endonuclease I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (26/37);  minor base C (1/0);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGCGGCAACTTTATGTACAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACC  >  NC_000913/3090741‑3091212
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
gtcgcggcatcggtATGTACAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGGCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTTCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACGCTCTCGCCAGCAAACGCAGCTTTTCAACGCAGGGAACAATATGTATCCGGTTACCTACTGGAAGTGCGAGCGCGATgttctcatc                                                                                                                                                                                                                                <  2:67415/238‑9 (MQ=255)
    cgGCAACTTTATGTACAGCCAGTGGAATGGCGGTGCAGGCCAGTACGGTCCATTCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGCACCACCAGCGCGTGCCCGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACCCTCTCCCCAGCAAACGCAGCTGTTCAACGCATGGAACCAGATGTATCCGGTTACCGACTGGAAGTGCGCGCGCGATGAACGCATcgcg                                                                                                                                                                                                                             >  2:147229/1‑250 (MQ=255)
     ggCAACTTTATGTACAGCCAGTGGAATGGCGGTGAAGGCAAGTGCGGACACAGCGACATGAAGGGCGAATTCAAAGAAAAAGCGGGCGAAACACAAGCGCGGGCACGCGGTGCCATTGCGCGCAACTACTTCTATAAGAGCGACCAATACAAACAGACACACTCTCGCCAGCAAACGCAACTGATCAAAGAATGGAAAAAGATGGATACGGTTACAGAATGGAAGAGCGAGCGCGAGGGACGCATCGCGa                                                                                                                                                                                                                            >  2:369284/1‑250 (MQ=255)
              aTGTACAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGGGCAgg                                                                                                                                                                                                                   <  2:379209/250‑1 (MQ=255)
                gTACAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCa                                                                                                                                                                                                                <  2:333260/251‑1 (MQ=255)
                 tACAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCa                                                                                                                                                                                                                <  2:185892/250‑1 (MQ=255)
                       cAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATaa                                                                                                                                                                                                         <  1:338225/251‑1 (MQ=255)
                             aaTGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTa                                                                                                                                                                                                  <  2:176958/251‑1 (MQ=255)
                                ggCGGTGAAGGCCAGTGCGGTCAATGCGCCATGAAGGTCGATTGCAAAGAAAAAGCTGCCGAACCACCAGCGCGGGCTCGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACACTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACGGGAAGTTCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAAGCATAACCCGTAtgt                                                                                                                                                                                               <  2:443332/251‑1 (MQ=255)
                                   ggTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAAc                                                                                                                                                                                                        <  2:226737/240‑1 (MQ=255)
                                   ggTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAAc                                                                                                                                                                                                        >  1:226737/1‑240 (MQ=255)
                                    gTGAAGGCCAGTACGGTCAATTCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCaa                                                                                                                                                                                           <  1:387450/251‑1 (MQ=255)
                                      gAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAAcg                                                                                                                                                                                         >  1:160474/1‑251 (MQ=255)
                                          gCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGGGCAACGCGCt                                                                                                                                                                                     >  1:373737/1‑251 (MQ=255)
                                           ccAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCGCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCtt                                                                                                                                                                                    <  1:284951/251‑1 (MQ=255)
                                            cAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCTTATGTGCAACGCGCTTg                                                                                                                                                                                   <  2:207955/251‑1 (MQ=255)
                                                               tGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATg                                                                                                                                                                                                                                        <  1:198024/180‑1 (MQ=255)
                                                               tGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATg                                                                                                                                                                                                                                        >  2:198024/1‑180 (MQ=255)
                                                               ggcaGGTCGATTTCAAAGAAAAAGCGGCCGAACCGCCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCGAGATGCGCGACCAATACAACCTGACGCTCTCCCGCAAGCAATCGCAGCAGTGCAACGCATGGAACAAGATGTAACCGGTTACCGACTAGAATTACGAGCACTATGATCGCACCGCGAAGGTGCAGGGGAAACAAACCTCGTTTTTGCAACTCGCGCGCCAGGCGCGAAagggcgaa                                                                                                                                                                <  2:436191/248‑7 (MQ=255)
                                                                gAAGGTCGATGTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAAc                                                                                                                                                               <  2:352094/251‑1 (MQ=255)
                                                                 aaGGTCGATTTCAAAGAAAAAGCTGCCGAACCCCCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGCAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGACAGAGCTAAcc                                                                                                                                                              >  2:121589/1‑251 (MQ=255)
                                                                      cGATTTCAAAGCAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGCAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTacac                                                                                                                                                         >  2:395436/1‑251 (MQ=255)
                                                                          ttCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAg                                                                                                                                                      <  2:208206/250‑1 (MQ=255)
                                                                               aGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCAACTACGTCTAGATGCGCGACCAATACAACCTTACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATTAACGCATCGCGTAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTGGCGGGAt                                                                                                                                                <  2:411861/251‑1 (MQ=255)
                                                                                  aaaaGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACACGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCt                                                                                                                                             >  1:93229/1‑251 (MQ=255)
                                                                                  aaaaGCTGACGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGCTGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTACCCTACACTAGCGGGATTCt                                                                                                                                             >  2:170248/1‑251 (MQ=255)
                                                                                    aaGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCttt                                                                                                                                           >  1:141225/1‑251 (MQ=255)
                                                                                         gCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAAc                                                                                                                                                                                                                                     >  1:162579/1‑157 (MQ=255)
                                                                                         gCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAAc                                                                                                                                                                                                                                     <  2:162579/157‑1 (MQ=255)
                                                                                         gCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGtt                                                                                                                                      <  1:378806/251‑1 (MQ=255)
                                                                                           cGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGCAGTGCGAgcgc                                                                                                                                                                                                                                            >  1:443476/1‑148 (MQ=255)
                                                                                                 accaGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTa                                                                                                                              <  1:121589/251‑1 (MQ=255)
                                                                                                    aGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTAccc                                                                                                                           <  1:127777/251‑1 (MQ=255)
                                                                                                        cgTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTGTTGTTAACCCCTACCCCAcg                                                                                                                       <  2:230112/251‑1 (MQ=255)
                                                                                                                    gCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTGTTTGTTAACCCCTACCCCACGCGTACAACcgcg                                                                                                           <  1:324536/251‑1 (MQ=255)
                                                                                                                      cATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTg                                                                                                         >  2:290758/1‑251 (MQ=255)
                                                                                                                          gcgtgcACCTACTTCTATATGCGCGAACAATACAACCTGACACTCTCTCGCCAGCAAACGAAGCAGTTCAAAGCATGGAACAAGATGTATCCTGATACCGACAGGATGTGGGAGCGCGATGACCGCATCGCGAAGGTGCAGGGCTATCATAACCCGTATGTGCAACGCGCGTGCCAAGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTTTTAACACCTTCCCCACACTGACAACCGGGTggggg                                                                                                     >  2:423557/5‑250 (MQ=255)
                                                                                                                           cgcgcACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTgggg                                                                                                      >  1:109155/1‑249 (MQ=255)
                                                                                                                             cgcACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAgac                                                                                                  <  1:381790/251‑1 (MQ=255)
                                                                                                                                aCCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACATCCGCGTGGGGAgacgac                                                                                               <  1:448195/251‑1 (MQ=255)
                                                                                                                                 ccTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGacgacg                                                                                              <  1:233326/251‑1 (MQ=255)
                                                                                                                                 ccTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGacgacg                                                                                              >  1:128002/1‑251 (MQ=255)
                                                                                                                                                                  cacTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTa                                                             <  1:187183/251‑1 (MQ=255)
                                                                                                                                                                   acTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGGATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCGACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGGCGCGGATTTTTAGCTATGCTTATCCCCCGCATTTat                                                            <  2:25757/251‑1 (MQ=255)
                                                                                                                                                                            ccAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGa                                                    >  1:104299/1‑250 (MQ=255)
                                                                                                                                                                            ccAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGa                                                    >  1:348571/1‑250 (MQ=255)
                                                                                                                                                                                      cAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCa                                         <  1:146279/251‑1 (MQ=255)
                                                                                                                                                                                          tGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGcc                                      >  2:172013/1‑250 (MQ=255)
                                                                                                                                                                                            ttCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATACGTATCCCCCGCATTTATCATCCTGAACCACTGCCCAGCCAtt                                   >  2:38940/1‑251 (MQ=255)
                                                                                                                                                                                                aCGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATtctc                                >  1:180104/1‑250 (MQ=255)
                                                                                                                                                                                                aCGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTca                               <  1:331564/251‑1 (MQ=255)
                                                                                                                                                                                                 cGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTcac                              >  1:81281/1‑251 (MQ=255)
                                                                                                                                                                                                   cATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACAt                            >  1:42333/1‑251 (MQ=255)
                                                                                                                                                                                                       gAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGa                                            <  1:582/231‑1 (MQ=255)
                                                                                                                                                                                                       gAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGa                                            >  2:582/1‑231 (MQ=255)
                                                                                                                                                                                                         aCAAGATGTATCCGGTTACCGACTGGAAGGGCGAGCGCGATGAACGCATCGCGAAGGGGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCACGCt                      >  2:171802/1‑251 (MQ=255)
                                                                                                                                                                                                          cAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACGAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCgcgcgt                     <  2:360462/251‑3 (MQ=255)
                                                                                                                                                                                                           aaGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCttt                    >  1:407079/1‑251 (MQ=255)
                                                                                                                                                                                                           aaGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCttt                    <  2:160474/251‑1 (MQ=255)
                                                                                                                                                                                                            aGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGATCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTg                   <  2:92698/251‑1 (MQ=255)
                                                                                                                                                                                                                    ccGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATg           <  2:104299/251‑1 (MQ=255)
                                                                                                                                                                                                                            cGACTGGAAGTGCGAGCGCGATGAACGCATCGGGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAAc   <  2:45230/251‑1 (MQ=255)
                                                                                                                                                                                                                            cGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAAc   <  2:407079/251‑1 (MQ=255)
                                                                                                                                                                                                                             gACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAAcc  <  1:139993/251‑1 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
ATCGCGGCAACTTTATGTACAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACC  >  NC_000913/3090741‑3091212

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