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

CTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCT  >  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)
                                                                                               tGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAAtgtg                                                                                                                                                <  1:23963/250‑1 (MQ=255)
                                                                                                  tGGCCGTTTTGTTCTTGCGAGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCAACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATAGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGATGATTGAGGTCATGCTCGTTGGCTGCGGTGGTGACCAGGCTGTGGGTCAGGCCACACTTGGGATAGAGACTAATGTGGGcg                                                                                                                                            >  2:208356/1‑250 (MQ=14)
                                                                                                   ggCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACTCCAATGTGGGCCt                                                                                                                                           <  1:94045/251‑1 (MQ=255)
                                                                                                    gCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCtt                                                                                                                                          <  1:89563/251‑1 (MQ=255)
                                                                                                        ttttGTTCTTGCGCGGATTCTGCTCCAAGGTTTTTACCTCGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCATCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATTCTCGTTGGCCGCTGTTGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCAt                                                                                                                                       <  1:290141/250‑1 (MQ=14)
                                                                                                           tGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCCTCGGCTGCGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTAAGGCCACTCTCGGCATCGCACCCAATGTGGGCATTTATGcc                                                                                                                                    >  2:231755/1‑250 (MQ=14)
                                                                                                              tCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCaa                                                                                                                                  >  2:124228/1‑249 (MQ=255)
                                                                                                                    gcggattcttgtgtaaGGTTTTTACCTTGCCGGGTCGCTCGGCGATCAGCCCGTCCACATCCACCTCGGCCAGCCCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACAAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCAc                                                                                                                          <  2:119037/237‑1 (MQ=11)
                                                                                                                    gcgGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCCGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCCATGGGGGCCTTCATGCCAAAGGGCCAc                                                                                                                          >  2:43188/1‑251 (MQ=255)
                                                                                                                          tCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTg                                                                                                                    >  1:104585/1‑251 (MQ=39)
                                                                                                                            tGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGCCACCAATGGGGGCCTTCATGCCAAAGTGCCACTGATTGcc                                                                                                                  >  2:195257/1‑251 (MQ=255)
                                                                                                                             gCTTCAAGGTTTTTACCTTGCCGTGACTGTCGGCGATCAGCCGGTCCACAGCCACCTCGGCCAGCTCATCGCGCTGTGGCGCTCCTTGGGAGCCGGCATCGGCTGAGATAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGc                                                                                                                   <  1:88003/249‑1 (MQ=11)
                                                                                                                                tCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTAGGCAAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGAATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCttt                                                                                                               <  2:358091/250‑1 (MQ=255)
                                                                                                                                  aagggtttcaccgcgCCGGTACGCTGCGCGACCAGCCGTTCCACATCCACCTCTGCCAGCTCCGCGCGCTGTGGCGCTCCGTGGTAGCCGGCATCGGCTGATACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCCCGTTGGCCGCGGTGGTGACCATGCTGTGGGTCAGGCCACTCTTGGCATCGTCACCAATGTGGGCCTTCATGCCAAAGTGCCTCTGATTGCCTTTCtt                                                                                                            <  2:67526/237‑1 (MQ=11)
                                                                                                                                  aaGGTTTTTACCTGGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCtt                                                                                                            <  2:212150/251‑1 (MQ=37)
                                                                                                                                    ggTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACATCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGATTCTCTCCATGAGGCAGATTACCCAGCTGCTTGAGGTAATGCTCGTTGGCCGAGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGAAACAAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTGCTTgg                                                                                                          >  1:68005/1‑251 (MQ=11)
                                                                                                                                      tttttACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTc                                                                                                        >  1:168295/1‑251 (MQ=25)
                                                                                                                                          tACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCTAGCTCCTCGCGCTGTGGCGCTCCTTGGAAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCAGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACCCCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGAt                                                                                                    >  1:325354/1‑251 (MQ=255)
                                                                                                                                          tACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGACTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGGGGGCCTTCATGCCACAGTGCCACTGATTGCCTTTCTTGGTCTGAt                                                                                                    >  1:360090/1‑251 (MQ=255)
                                                                                                                                             cTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTCGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCGTTCTTGGTCTGATGc                                                                                                  >  1:13305/1‑250 (MQ=25)
                                                                                                                                                        gCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATcg                                                                                      <  1:396595/251‑1 (MQ=18)
                                                                                                                                                          tCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACCAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGACTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATctcg                                                                                    >  2:35037/1‑248 (MQ=16)
                                                                                                                                                           tgtcgataagccctgccACACCCACCTCGGCCCGCTCCTCTCGCTGTTGCTCTCTTTGGTAGCCGGAATCGGCTGAGACAAATTGCGTCTCTCCATGAAGCAGATTACCCAGCTGACTGAGGCCATGCTCTTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCGTTCTTGGTCTGATGCATCTCCGGATcgc                                                                                     <  2:137741/234‑1 (MQ=11)
                                                                                                                                                                      cattCCACATCCACCTCGGCCAGCGCCTCCCGCTGGGGCCCTCCTTGGTAGCCGGCAGCGGCTGAGACAAATTGCTGCTCTCCATGAAGCAGATTACTCAGCTGATTGATGTCATGCGCGTTGGCCGCGGGGGTGACCAGGCTGTGGTTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGAGTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCttt                                                                        <  2:252027/248‑1 (MQ=11)
                                                                                                                                                                      cAGTCCACTTCCCCCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCttt                                                                        <  1:397851/251‑1 (MQ=18)
                                                                                                                                                                       aGTCCACAGCCTCCTCGGCCAGGTCCTCGCGCTGTGGCGCACCGTTGTAGCCGGCATCGTCTGAGCCAAATTGCTCCTCTCCCTGAAGCAGATTACCCAGCTGATTTAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTTATGCATCTCCGGATCGCGTTGCTGCTCTTTg                                                                       <  2:277612/251‑1 (MQ=11)
                                                                                                                                                                          tcataTCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCAACGGCTGAGACAAATTGATTCTCTCCATGAAGCAGATTAACCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGTCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCCTTGCTGCTCTTTGtt                                                                     <  2:189589/246‑1 (MQ=11)
                                                                                                                                                                            acaTCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCtt                                                                  <  1:308705/251‑1 (MQ=18)
                                                                                                                                                                             ctccccccgcGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCGGAGACAAATTGCTTCGCTCCAGGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATGGCCTTTCTTGGTCTGATGCATCTCCGGAGCGCGTTGCTGCTCTTTGTTCt                                                                   <  2:429328/240‑1 (MQ=11)
                                                                                                                                                                               gccACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGt                                                               <  1:163865/250‑1 (MQ=18)
                                                                                                                                                                                      cGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTg                                                        <  2:49586/251‑1 (MQ=21)
                                                                                                                                                                                                  cgcgCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAatg                                            >  2:102478/1‑251 (MQ=21)
                                                                                                                                                                                                    cgcTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAatgat                                          <  1:9038/251‑1 (MQ=21)
                                                                                                                                                                                                        gtgGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGAtggtg                                      >  1:36601/1‑251 (MQ=21)
                                                                                                                                                                                                              cTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAAGGATGGTGGCATcc                                >  2:98387/1‑251 (MQ=21)
                                                                                                                                                                                                                         ccGGCATCTGCTGAGACAAATTGCCTCTCTCCATGCAGCAGATTACCCCTCGGCTTGTGGTCATGCTCGTTGGCAGCGGTTGGGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGTCATCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTTCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGcc                     <  2:119809/251‑1 (MQ=11)
                                                                                                                                                                                                                          cGGCATCGGCTGAGACAAAGTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGACTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCt                    <  1:35037/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aTTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCt  >  1:461793/1‑251 (MQ=38)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
CTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCT  >  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: