Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,428,087 C→T E214K (GAG→AAG)  insH5 ← IS5 family transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,428,0870CT98.2% 139.9 / NA 57E214K (GAG→AAG) insH5IS5 family transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (29/26);  minor base G (0/1);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCT  >  NC_000913/1427845‑1428331
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
cTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCAt                                                                                                                                                                                                                                               >  1:33239/1‑251 (MQ=255)
   gTATCTGGCTTTCACGAAGCCGAACAGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGCTGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGACTGGGACAACTTGCTTCCCTCCATGa                                                                                                                                                                                                                                             >  2:320034/1‑250 (MQ=14)
     aTCTGGCTGTCACGAAGCCGAACTGCCGCTTGATGACGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTGTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAg                                                                                                                                                                                                                                           <  2:389945/250‑1 (MQ=255)
      tCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGc                                                                                                                                                                                                                                          >  1:315050/1‑250 (MQ=255)
            tttCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTAc                                                                                                                                                                                                                                   >  1:265625/1‑251 (MQ=255)
                     gctgAACTGCCCCTTGATGATGCGAAACGGGTGCCCCACCCTGGCTCGGATGTTGGTTTTCATTTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTCCCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGAt                                                                                                                                                                                                                          <  2:394456/248‑1 (MQ=14)
                     gCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGAt                                                                                                                                                                                                                          <  1:397863/251‑1 (MQ=255)
                        gAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTg                                                                                                                                                                                                                        >  1:181876/1‑250 (MQ=255)
                                  ttGATGATGCGAAACGGGGGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGCTGTTGATGGCCGTTTTGTTCTTGCGCGGCTTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAAACCAATGGCTTATCTCCAAGAAGCAGATTACCCGGCTGATTAAGGCCATGCTc                                                                                                                                                                                                             >  2:350838/1‑251 (MQ=14)
                                   tgatgaTGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCt                                                                                                                                                                                                              <  1:9119/249‑1 (MQ=255)
                                   tgatgaTGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCg                                                                                                                                                                                                            >  1:212150/1‑251 (MQ=255)
                                        attcGAAACGGGTTCTCCTCCCTGGCACGGATGCGCGCTTTCATGTATTCGATGTTGATGTCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTGTTACCTTGCCGTGACGCTCGGCGATCAGCCTGGCCACATCCCCCTCGGCCAGCTCCTCGGGCTGTGGCGCTCCGTGGTAGCCGGCATCGGCTGATTCAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGc                                                                                                                                                                                                       <  1:200183/248‑1 (MQ=14)
                                            gAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCcgcg                                                                                                                                                                                                   >  1:212468/1‑251 (MQ=255)
                                                             cTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCtgtg                                                                                                                                                                                  >  1:455903/1‑251 (MQ=255)
                                                                gCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGt                                                                                                                                                                               >  1:80287/1‑251 (MQ=255)
                                                                    ggATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAgg                                                                                                                                                                           >  1:300225/1‑251 (MQ=255)
                                                                          tGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGCTTTTACCTTACCGGGACGCTCGGCGATCAGCCAGTCCACCTCCACCTCGGCAAGCTCCTCGCGATGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCCTGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCCGGCTGTGGGTCAGGCCActc                                                                                                                                                                     >  1:463101/1‑251 (MQ=14)
                                                                            gCTTTCATGTTTTTGGTGTTGAGTGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTGGCCGGGGCGCTCGGCGAGCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCTCGCTGTGGCGCGCCGTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCAGGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCtt                                                                                                                                                                   <  2:280743/251‑1 (MQ=11)
                                                                                          gATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGATGTGGGTCAGGCCACTCTTGGCATCGaaccca                                                                                                                                                      >  2:306790/1‑245 (MQ=255)
                                                                                               tGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AAtgtg                                                                                                                                                <  1:23963/250‑1 (MQ=255)
                                                                                                  tGGCCGTTTTGTTCTTGCGAGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCAACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATAGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGATGATTGAGGTCATGCTCGTTGGCTGCGGTGGTGACCAGGCTGTGGGTCAGGCCACACTTGGGATAGAGACT‑AATGTGGGcg                                                                                                                                            >  2:208356/1‑250 (MQ=14)
                                                                                                   ggCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACTCC‑AATGTGGGCCt                                                                                                                                           <  1:94045/251‑1 (MQ=255)
                                                                                                    gCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCtt                                                                                                                                          <  1:89563/251‑1 (MQ=255)
                                                                                                        ttttGTTCTTGCGCGGATTCTGCTCCAAGGTTTTTACCTCGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCATCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATTCTCGTTGGCCGCTGTTGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCAt                                                                                                                                       <  1:290141/250‑1 (MQ=14)
                                                                                                           tGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCCTCGGCTGCGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTAAGGCCACTCTCGGCATCG‑CACCCAATGTGGGCATTTATGcc                                                                                                                                    >  2:231755/1‑250 (MQ=14)
                                                                                                              tCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCaa                                                                                                                                  >  2:124228/1‑249 (MQ=255)
                                                                                                                    gcggattcttgtgtaaGGTTTTTACCTTGCCGGGTCGCTCGGCGATCAGCCCGTCCACATCCACCTCGGCCAGCCCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACAAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCAc                                                                                                                          <  2:119037/237‑1 (MQ=11)
                                                                                                                    gcgGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCCGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑CATGGGGGCCTTCATGCCAAAGGGCCAc                                                                                                                          >  2:43188/1‑251 (MQ=255)
                                                                                                                          tCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTg                                                                                                                    >  1:104585/1‑251 (MQ=39)
                                                                                                                            tGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGCCACC‑AATGGGGGCCTTCATGCCAAAGTGCCACTGATTGcc                                                                                                                  >  2:195257/1‑251 (MQ=255)
                                                                                                                             gCTTCAAGGTTTTTACCTTGCCGTGACTGTCGGCGATCAGCCGGTCCACAGCCACCTCGGCCAGCTCATCGCGCTGTGGCGCTCCTTGGGAGCCGGCATCGGCTGAGATAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGc                                                                                                                   <  1:88003/249‑1 (MQ=11)
                                                                                                                                tCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTAGGCAAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGAATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCttt                                                                                                               <  2:358091/250‑1 (MQ=255)
                                                                                                                                  aagggtttcaccgcgCCGGTACGCTGCGCGACCAGCCGTTCCACATCCACCTCTGCCAGCTCCGCGCGCTGTGGCGCTCCGTGGTAGCCGGCATCGGCTGATACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCCCGTTGGCCGCGGTGGTGACCATGCTGTGGGTCAGGCCACTCTTGGCATCGTCACC‑AATGTGGGCCTTCATGCCAAAGTGCCTCTGATTGCCTTTCtt                                                                                                            <  2:67526/237‑1 (MQ=11)
                                                                                                                                  aaGGTTTTTACCTGGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCtt                                                                                                            <  2:212150/251‑1 (MQ=37)
                                                                                                                                    ggTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACATCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGATTCTCTCCATGAGGCAGATTACCCAGCTGCTTGAGGTAATGCTCGTTGGCCGAGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGAAACA‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTGCTTgg                                                                                                          >  1:68005/1‑251 (MQ=11)
                                                                                                                                      tttttACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTc                                                                                                        >  1:168295/1‑251 (MQ=25)
                                                                                                                                          tACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCTAGCTCCTCGCGCTGTGGCGCTCCTTGGAAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCAGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACCCC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGAt                                                                                                    >  1:325354/1‑251 (MQ=255)
                                                                                                                                          tACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGACTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGGGGGCCTTCATGCCACAGTGCCACTGATTGCCTTTCTTGGTCTGAt                                                                                                    >  1:360090/1‑251 (MQ=255)
                                                                                                                                             cTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTCGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCGTTCTTGGTCTGATGc                                                                                                  >  1:13305/1‑250 (MQ=25)
                                                                                                                                                        gCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATcg                                                                                      <  1:396595/251‑1 (MQ=18)
                                                                                                                                                          tCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACCAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGACTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATctcg                                                                                    >  2:35037/1‑248 (MQ=16)
                                                                                                                                                           tgtcgataagccctgccACACCCACCTCGGCCCGCTCCTCTCGCTGTTGCTCTCTTTGGTAGCCGGAATCGGCTGAGACAAATTGCGTCTCTCCATGAAGCAGATTACCCAGCTGACTGAGGCCATGCTCTTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCGTTCTTGGTCTGATGCATCTCCGGATcgc                                                                                     <  2:137741/234‑1 (MQ=11)
                                                                                                                                                                      cattCCACATCCACCTCGGCCAGCGCCTCCCGCTGGGGCCCTCCTTGGTAGCCGGCAGCGGCTGAGACAAATTGCTGCTCTCCATGAAGCAGATTACTCAGCTGATTGATGTCATGCGCGTTGGCCGCGGGGGTGACCAGGCTGTGGTTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGAGTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCttt                                                                        <  2:252027/248‑1 (MQ=11)
                                                                                                                                                                      cAGTCCACTTCCCCCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCttt                                                                        <  1:397851/251‑1 (MQ=18)
                                                                                                                                                                       aGTCCACAGCCTCCTCGGCCAGGTCCTCGCGCTGTGGCGCACCGTTGTAGCCGGCATCGTCTGAGCCAAATTGCTCCTCTCCCTGAAGCAGATTACCCAGCTGATTTAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTTATGCATCTCCGGATCGCGTTGCTGCTCTTTg                                                                       <  2:277612/251‑1 (MQ=11)
                                                                                                                                                                          tcataTCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCAACGGCTGAGACAAATTGATTCTCTCCATGAAGCAGATTAACCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGTCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCCTTGCTGCTCTTTGtt                                                                     <  2:189589/246‑1 (MQ=11)
                                                                                                                                                                            acaTCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCtt                                                                  <  1:308705/251‑1 (MQ=18)
                                                                                                                                                                             ctccccccgcGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCGGAGACAAATTGCTTCGCTCCAGGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATGGCCTTTCTTGGTCTGATGCATCTCCGGAGCGCGTTGCTGCTCTTTGTTCt                                                                   <  2:429328/240‑1 (MQ=11)
                                                                                                                                                                               gccACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGt                                                               <  1:163865/250‑1 (MQ=18)
                                                                                                                                                                                      cGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTg                                                        <  2:49586/251‑1 (MQ=21)
                                                                                                                                                                                                  cgcgCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAatg                                            >  2:102478/1‑251 (MQ=21)
                                                                                                                                                                                                    cgcTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAatgat                                          <  1:9038/251‑1 (MQ=21)
                                                                                                                                                                                                        gtgGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGAtggtg                                      >  1:36601/1‑251 (MQ=21)
                                                                                                                                                                                                              cTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAAGGATGGTGGCATcc                                >  2:98387/1‑251 (MQ=21)
                                                                                                                                                                                                                         ccGGCATCTGCTGAGACAAATTGCCTCTCTCCATGCAGCAGATTACCCCTCGGCTTGTGGTCATGCTCGTTGGCAGCGGTTGGGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGTCATC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTTCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGcc                     <  2:119809/251‑1 (MQ=11)
                                                                                                                                                                                                                          cGGCATCGGCTGAGACAAAGTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGACTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCt                    <  1:35037/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aTTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCt  >  1:461793/1‑251 (MQ=38)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
CTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACC‑AATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCT  >  NC_000913/1427845‑1428331

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