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

ATCGCGGCAACTTTATGTACAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACC  >  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)
                             aaTGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTa                                                                                                                                                                                                  <  2:176958/251‑1 (MQ=255)
                                ggCGGTGAAGGCCAGTGCGGTCAATGCGCCATGAAGGTCGATTGCAAAGAAAAAGCTGCCGAACCACCAGCGCGGGCTCGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACACTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACGGGAAGTTCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAAGCATAACC‑CGTAtgt                                                                                                                                                                                               <  2:443332/251‑1 (MQ=255)
                                   ggTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAAc                                                                                                                                                                                                        <  2:226737/240‑1 (MQ=255)
                                   ggTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAAc                                                                                                                                                                                                        >  1:226737/1‑240 (MQ=255)
                                    gTGAAGGCCAGTACGGTCAATTCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCaa                                                                                                                                                                                           <  1:387450/251‑1 (MQ=255)
                                      gAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAAcg                                                                                                                                                                                         >  1:160474/1‑251 (MQ=255)
                                          gCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGGGCAACGCGCt                                                                                                                                                                                     >  1:373737/1‑251 (MQ=255)
                                           ccAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCGCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCtt                                                                                                                                                                                    <  1:284951/251‑1 (MQ=255)
                                            cAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CTTATGTGCAACGCGCTTg                                                                                                                                                                                   <  2:207955/251‑1 (MQ=255)
                                                               tGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATg                                                                                                                                                                                                                                        <  1:198024/180‑1 (MQ=255)
                                                               tGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATg                                                                                                                                                                                                                                        >  2:198024/1‑180 (MQ=255)
                                                               ggcaGGTCGATTTCAAAGAAAAAGCGGCCGAACCGCCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCGAGATGCGCGACCAATACAACCTGACGCTCTCCCGCAAGCAATCGCAGCAGTGCAACGCATGGAACAAGATGTAACCGGTTACCGACTAGAATTACGAGCACTATGATCGCACCGCGAAGGTGCAGGGGAAACA‑AACCTCGTTTTTGCAACTCGCGCGCCAGGCGCGAAagggcgaa                                                                                                                                                                <  2:436191/248‑7 (MQ=255)
                                                                gAAGGTCGATGTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAAc                                                                                                                                                               <  2:352094/251‑1 (MQ=255)
                                                                 aaGGTCGATTTCAAAGAAAAAGCTGCCGAACCCCCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGCAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGACAGAGCTAAcc                                                                                                                                                              >  2:121589/1‑251 (MQ=255)
                                                                      cGATTTCAAAGCAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGCAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTacac                                                                                                                                                         >  2:395436/1‑251 (MQ=255)
                                                                          ttCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAg                                                                                                                                                      <  2:208206/250‑1 (MQ=255)
                                                                               aGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCAACTACGTCTAGATGCGCGACCAATACAACCTTACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATTAACGCATCGCGTAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTGGCGGGAt                                                                                                                                                <  2:411861/251‑1 (MQ=255)
                                                                                  aaaaGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACA‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCt                                                                                                                                             >  1:93229/1‑251 (MQ=255)
                                                                                  aaaaGCTGACGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGCTGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTACCCTACACTAGCGGGATTCt                                                                                                                                             >  2:170248/1‑251 (MQ=255)
                                                                                    aaGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCttt                                                                                                                                           >  1:141225/1‑251 (MQ=255)
                                                                                         gCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAAc                                                                                                                                                                                                                                     >  1:162579/1‑157 (MQ=255)
                                                                                         gCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAAc                                                                                                                                                                                                                                     <  2:162579/157‑1 (MQ=255)
                                                                                         gCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGtt                                                                                                                                      <  1:378806/251‑1 (MQ=255)
                                                                                           cGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGCAGTGCGAgcgc                                                                                                                                                                                                                                            >  1:443476/1‑148 (MQ=255)
                                                                                                 accaGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTa                                                                                                                              <  1:121589/251‑1 (MQ=255)
                                                                                                    aGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTAccc                                                                                                                           <  1:127777/251‑1 (MQ=255)
                                                                                                        cgTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTGTTGTTAACCCCTACCCCAcg                                                                                                                       <  2:230112/251‑1 (MQ=255)
                                                                                                                    gCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTGTTTGTTAACCCCTACCCCACGCGTACAACcgcg                                                                                                           <  1:324536/251‑1 (MQ=255)
                                                                                                                      cATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTg                                                                                                         >  2:290758/1‑251 (MQ=255)
                                                                                                                          gcgtgcACCTACTTCTATATGCGCGAACAATACAACCTGACACTCTCTCGCCAGCAAACGAAGCAGTTCAAAGCATGGAACAAGATGTATCCTGATACCGACAGGATGTGGGAGCGCGATGACCGCATCGCGAAGGTGCAGGGCTATCATAACC‑CGTATGTGCAACGCGCGTGCCAAGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTTTTAACACCTTCCCCACACTGACAACCGGGTggggg                                                                                                     >  2:423557/5‑250 (MQ=255)
                                                                                                                           cgcgcACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTgggg                                                                                                      >  1:109155/1‑249 (MQ=255)
                                                                                                                             cgcACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAgac                                                                                                  <  1:381790/251‑1 (MQ=255)
                                                                                                                                aCCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACATCCGCGTGGGGAgacgac                                                                                               <  1:448195/251‑1 (MQ=255)
                                                                                                                                 ccTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGacgacg                                                                                              <  1:233326/251‑1 (MQ=255)
                                                                                                                                 ccTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGacgacg                                                                                              >  1:128002/1‑251 (MQ=255)
                                                                                                                                                                  cacTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTa                                                             <  1:187183/251‑1 (MQ=255)
                                                                                                                                                                   acTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGGATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCGACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGGCGCGGATTTTTAGCTATGCTTATCCCCCGCATTTat                                                            <  2:25757/251‑1 (MQ=255)
                                                                                                                                                                            ccAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGa                                                    >  1:104299/1‑250 (MQ=255)
                                                                                                                                                                            ccAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGa                                                    >  1:348571/1‑250 (MQ=255)
                                                                                                                                                                                      cAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCa                                         <  1:146279/251‑1 (MQ=255)
                                                                                                                                                                                          tGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGcc                                      >  2:172013/1‑250 (MQ=255)
                                                                                                                                                                                            ttCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATACGTATCCCCCGCATTTATCATCCTGAACCACTGCCCAGCCAtt                                   >  2:38940/1‑251 (MQ=255)
                                                                                                                                                                                                aCGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATtctc                                >  1:180104/1‑250 (MQ=255)
                                                                                                                                                                                                aCGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTca                               <  1:331564/251‑1 (MQ=255)
                                                                                                                                                                                                 cGCATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTcac                              >  1:81281/1‑251 (MQ=255)
                                                                                                                                                                                                   cATGGAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACAt                            >  1:42333/1‑251 (MQ=255)
                                                                                                                                                                                                       gAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGa                                            <  1:582/231‑1 (MQ=255)
                                                                                                                                                                                                       gAACAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGa                                            >  2:582/1‑231 (MQ=255)
                                                                                                                                                                                                         aCAAGATGTATCCGGTTACCGACTGGAAGGGCGAGCGCGATGAACGCATCGCGAAGGGGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCACGCt                      >  2:171802/1‑251 (MQ=255)
                                                                                                                                                                                                          cAAGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACGAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCgcgcgt                     <  2:360462/251‑3 (MQ=255)
                                                                                                                                                                                                           aaGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCttt                    >  1:407079/1‑251 (MQ=255)
                                                                                                                                                                                                           aaGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCttt                    <  2:160474/251‑1 (MQ=255)
                                                                                                                                                                                                            aGATGTATCCGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGATCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTg                   <  2:92698/251‑1 (MQ=255)
                                                                                                                                                                                                                    ccGGTTACCGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATg           <  2:104299/251‑1 (MQ=255)
                                                                                                                                                                                                                            cGACTGGAAGTGCGAGCGCGATGAACGCATCGGGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAAc   <  2:45230/251‑1 (MQ=255)
                                                                                                                                                                                                                            cGACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAAc   <  2:407079/251‑1 (MQ=255)
                                                                                                                                                                                                                             gACTGGAAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAAcc  <  1:139993/251‑1 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
ATCGCGGCAACTTTATGTACAGCCAGTGGAATGGCGGTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACC‑CGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATGCCGCCAACC  >  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: ‑