Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,428,107 C→T G207D (GGT→GAT)  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,1070CT96.0% 125.3 / NA 51G207D (GGT→GAT) insH5IS5 family transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (27/21);  minor base G (0/2);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTG  >  NC_000913/1427869‑1428352
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     
gAACTGCCGCTTGATGATGCGCAACGGGTGCTCCACCCTGGCACGGATGCTGGCCTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTCCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGAGACGCTCGGCGATCAGCCAGTCCACCTCAACCTCGGCCAGCTACTCGCGCTGTGGCGCTCCTTGGTAGCCTGCATCGGCTGAGACCAATTGCTCCTCTCCATGAAGAAGATTACCCCGCTGGTTgg                                                                                                                                                                                                                                           >  2:538755/1‑250 (MQ=17)
       cggttGATGATGCGAATCGGGTGCTCCACCCTGGCACGGATGCTGGCGTTCATGTAGTCGATGGTGATGGCCGTTTTGTTCTTGCGCTGATTCTGCTCCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCAGGAAGCAGATTATCCAGCTGATTGTGGTCATg                                                                                                                                                                                                                                    <  1:306035/248‑1 (MQ=14)
             atgatgCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATGCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATGATCCAGCTGATTGAGGTCATGCTCGtt                                                                                                                                                                                                                              <  1:140742/251‑1 (MQ=255)
               gatgCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGGTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTgg                                                                                                                                                                                                                            >  1:226611/1‑251 (MQ=255)
                    gAAACGGGTGCTCCACCCTGGCTCGCATGCTGGCTTTCATGTATTCGATGTTAATGGCTGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCcgcg                                                                                                                                                                                                                       <  2:33404/251‑1 (MQ=255)
                     aaaCGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCgg                                                                                                                                                                                                                      >  1:371252/1‑251 (MQ=255)
                       aCGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCggtg                                                                                                                                                                                                                    <  1:368395/251‑1 (MQ=255)
                         gggTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCggtggt                                                                                                                                                                                                                  <  1:405876/251‑1 (MQ=255)
                                       ggCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCTGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTggg                                                                                                                                                                                                    <  1:538755/251‑1 (MQ=255)
                                         cACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGTTTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGCCAAATTGCTCCTCTCCATGACGCAGATTATCCAGCTGATTGCGGTCAAGCTCGTTGGCCGCGGTGGTGACCAGGCTGGGGGTc                                                                                                                                                                                                  >  2:25182/1‑251 (MQ=14)
                                                    gCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGTGTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCtt                                                                                                                                                                                       >  2:139271/1‑251 (MQ=255)
                                                    gCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCtt                                                                                                                                                                                       <  1:217749/251‑1 (MQ=255)
                                                       ttCATGTAGCCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAAGTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGCCATGCTCGTTGGCCGCGGGGGTGACCAGGCTGTGGGTCAGGCCACTCTTgg                                                                                                                                                                                     <  1:52293/250‑1 (MQ=255)
                                                            gTATTCGATGTTGATGGCCGTTTTGATCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCCTGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCCTCg                                                                                                                                                                                >  2:389723/1‑250 (MQ=255)
                                                                tCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTCCCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGacac                                                                                                                                                                            <  1:570273/250‑1 (MQ=255)
                                                                 cGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACAcc                                                                                                                                                                           >  1:496003/1‑250 (MQ=255)
                                                                    tGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTAGTTGGACGCGGTGGCGACCAGGATGTGGGTCCGGCAACGCTTGGAATCGACACCAAtg                                                                                                                                                                       >  2:102405/1‑251 (MQ=14)
                                                                            gCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCtt                                                                                                                                                               >  2:25281/1‑251 (MQ=255)
                                                                                 tttGTTCTTGCGCGGATTCTGCTTCAAGGATTTTACCTTGCCGGGACGCTCGGCGATCAGCCAATCCACATCCACCTCGGCCAGCCCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAACAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGCCCAGGCTGTGGGTCAGGCCACTCTTGGCAACGACAACAATGTAGGCCTTCATGc                                                                                                                                                          >  2:294727/1‑251 (MQ=14)
                                                                                    gTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCCTCGGCTGAGACCAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCCTGccc                                                                                                                                                        >  2:92669/1‑249 (MQ=255)
                                                                                    gTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGGGGGCCTTCATGcccc                                                                                                                                                       >  2:550161/1‑249 (MQ=255)
                                                                                                       ttCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGaca                                                                                                                                                                             <  1:169949/210‑1 (MQ=38)
                                                                                                       ttCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGaca                                                                                                                                                                             >  2:169949/1‑210 (MQ=38)
                                                                                                         cAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCCTGCCAAAGTGCCACTGCTTGCCTTTCt                                                                                                                                  >  2:268771/1‑251 (MQ=255)
                                                                                                         cAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGAGGGTCAGGCCACTCTTGGCATCGCCACCACTGTGGGCCTTCATGCCAAAGGACCACTGATTGCCTTTCt                                                                                                                                  >  2:529620/1‑251 (MQ=255)
                                                                                                           aGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAGGTGCCACTGATTGCCTTTCTTg                                                                                                                                >  2:52235/1‑251 (MQ=37)
                                                                                                                    ccTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGc                                                                                                                       >  1:133575/1‑251 (MQ=21)
                                                                                                                      ttGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGCGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCACAGTGCCACTGATTGCCTTTCTTGGTCTGATGCAt                                                                                                                     >  1:256012/1‑251 (MQ=21)
                                                                                                                       tGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCAt                                                                                                                     <  1:92669/250‑1 (MQ=25)
                                                                                                                          cGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGGACGGCCACTCGTGGCATCGGCCCCCATGTGGGCCTTCATGCCAAAGTGCCCCTGATTGCCTTTCTTGGTCTGATGCCTCTcc                                                                                                                 >  2:417956/1‑251 (MQ=11)
                                                                                                                            ggACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCgg                                                                                                               <  1:461136/251‑1 (MQ=21)
                                                                                                                                   cGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATcgcg                                                                                                         <  1:212962/250‑1 (MQ=18)
                                                                                                                                    ggCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGtt                                                                                                       <  1:529620/251‑1 (MQ=18)
                                                                                                                                     gcgctgaGCCAGTCCACGTCCCCCTCGGCCAGCTCCCCGCGCTGTGGCTCTCCTTGGTAGCCGGCATCGGTAGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGTTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACTCCAATGTGGGCCTTCTTGCCAAAGGGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGtta                                                                                                      <  2:419540/245‑2 (MQ=11)
                                                                                                                                     gcgatcaccccggctacaTCCACCCCGTCCAGTGCTTCGCGCTGAGGCTGTCATTGGGAGCCGGCATACGCTGAGACAAATTTCTCCACTGCTAGAAGCAGATTAGCCAACGGATTGGGGTAATTCTCGTTGGCCGCGGTAGTTACCAGGCTGTGTGTCAGGCCACTCTTGTCATCGTCTCCAATGTGGTCCTGCATTCCAAAGTGCCAATGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCTCGTtg                                                                                                      <  2:525494/236‑1 (MQ=2)
                                                                                                                                      cGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTtgc                                                                                                     >  1:14363/1‑251 (MQ=18)
                                                                                                                                           aGCCAGTCCACATCCCCCTCAGCCAGCTGCTCGCGCTGTAGCGCGACTTGGCAGCCGGCGTCGGATGAAACAAATGGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGAGGTGACCAGGCTGTGGGTCAGACCACTCTTGGCATCGACACCAATGTGAGCCTTCATGCCAACGGGACACTGATTGCCTTTCTTGGTCTGATGCATCTACGGATCGCGTTGCGGctc                                                                                                >  1:235787/1‑251 (MQ=11)
                                                                                                                                                         cACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGGGCCACTGCTTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGt                                                                                    >  2:384080/1‑249 (MQ=18)
                                                                                                                                                                           gcgcTGTGGCGCTCCATGGTAGCCGGCATCGGCTGAGACAAAGTGCTCCTCTCCATGAAGCAGATTAGCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAatga                                                                <  2:496003/251‑1 (MQ=21)
                                                                                                                                                                            cgcTGAGGCGTTCCTTGGTAGCCGGCAGCGGCTGAGAAAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAAGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCCCTCTAGGAATAGACTCCAATGTGGGCCTTCATGCCAAAGTGCAACTGATTGCCTTTCTTGGTCTGATGCATCCCCGGATCCACTTGAGGCTCTTTGTTCTTGGTAGAGCTGGGCGCCTccttgct                                                               >  2:406580/1‑245 (MQ=11)
                                                                                                                                                                                  ggCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGc                                                         <  1:528292/251‑1 (MQ=21)
                                                                                                                                                                                              tatccGGCATAGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGAATATCCAGCTGATTGAGGACCTGCTCGTTGGCCGAGGTGGTGAACAGGCTGAGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGGGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATAGCGTTTCTGCTCTTTTATCTTGGTTGAGCTGGGTGCCTCCatgatgtggtccactaccaaagt                                             >  2:315858/4‑234 (MQ=11)
                                                                                                                                                                                                      aTCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTggg                                     <  1:25182/251‑1 (MQ=38)
                                                                                                                                                                                                              agaCAAATTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATg                             <  1:502997/251‑1 (MQ=38)
                                                                                                                                                                                                                  aaaTTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGTGGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTGTGTTCTTGGTAGAGCTg                                                                             <  2:78379/199‑1 (MQ=255)
                                                                                                                                                                                                                  aaaTTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTg                                                                             >  1:78379/1‑199 (MQ=21)
                                                                                                                                                                                                                   aaTTGCTCCTCTCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGATCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCAATCTTGGCATCGACACCAATGTGGGCCTTCATGCCAACGTGCCACTGATTGCCTTTCTTGATCTGATGCATCTACGGCTAGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGAGCCTCAATGCTGGTGGCATCCACCAAGGTGCCTTAGGTCATCCTAACGcc                        >  2:20314/1‑251 (MQ=11)
                                                                                                                                                                                                                           ctctCCATGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCgg                >  1:393655/1‑251 (MQ=38)
                                                                                                                                                                                                                                 aTGAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGgccagcc          <  1:25281/251‑1 (MQ=38)
                                                                                                                                                                                                                                   gAAGCAGATTATCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGgccagcc          >  1:478577/1‑249 (MQ=38)
                                                                                                                                                                                                                                         gATTATCCAGCTGATTGAGGGCATGATTGTTGGGCGCGGTGGTGACCAGGCGGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGAGATTg  <  1:63769/251‑1 (MQ=11)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     
GAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTG  >  NC_000913/1427869‑1428352

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