Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,353,758 T→A D254V (GAT→GTT)  sapC ← putrescine ABC exporter membrane protein SapC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,353,7580TA98.5% 223.6 / NA 65D254V (GAT→GTT) sapCputrescine ABC exporter membrane protein SapC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (39/25);  minor base C (1/0);  total (40/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GACCGCGGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCATCACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCAT  >  NC_000913/1353517‑1354002
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
gACCGCGGATTTCACCTTCGGTTAAGGGCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGc                                                                                                                                                                                                                                              >  1:396414/1‑250 (MQ=255)
     cgGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACTGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCt                                                                                                                                                                                                                                        <  1:375768/251‑1 (MQ=255)
       gATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTcc                                                                                                                                                                                                                                      >  1:167791/1‑251 (MQ=255)
          ttCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTcccc                                                                                                                                                                                                                                    <  1:438668/250‑1 (MQ=255)
                         ggTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAgggg                                                                                                                                                                                                                      >  1:115768/1‑249 (MQ=255)
                              tACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTtgtg                                                                                                                                                                                                               >  2:282070/1‑251 (MQ=255)
                                   aCCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCg                                                                                                                                                                                                          >  1:117688/1‑251 (MQ=255)
                                    cccGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCg                                                                                                                                                                                                          >  2:518173/1‑250 (MQ=255)
                                        gTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCg                                                                                                                                                                                                          >  1:109463/1‑246 (MQ=255)
                                                        cACTCATCATCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTACCAGCGCAACTCCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGGGCGCCGAGATCGAGAAAGCCCAGcgcg                                                                                                                                                                                     <  1:268569/251‑1 (MQ=255)
                                                         aCTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCgcgg                                                                                                                                                                                    <  2:156494/251‑1 (MQ=255)
                                                                   ggTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCg                                                                                                                                                                            >  2:360142/1‑249 (MQ=255)
                                                                    gTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATAAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAACCAAGCGCGGCGATATCg                                                                                                                                                                            >  1:380829/1‑248 (MQ=255)
                                                                        tAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAAtt                                                                                                                                                                      <  1:601326/250‑1 (MQ=255)
                                                                                      aGGTTACGAATATCGAGTAATGACATTATTCCACCCCCGCAATAATAGCACGCCGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTgcc                                                                                                                                                       >  2:616940/1‑250 (MQ=255)
                                                                                        gTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGcgcg                                                                                                                                                     >  1:354091/1‑251 (MQ=255)
                                                                                           aCGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGAGCCGAGCTCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCgggg                                                                                                                                                  >  2:569616/1‑250 (MQ=255)
                                                                                             gactatCGAGTAATGGCATTATTCCACCCCCGCCATAATCGCACGACGAACGCCATCGCAAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGAGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTACCCATTCAGGCGCGTGGAGTTGTGCGCCGAGATCGAGACAGCCCAGCGCGTCCATATCGAGAATTGCCATCGACAGTGCGCggggg                                                                                                                                                 >  2:219921/4‑248 (MQ=255)
                                                                                               atatCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCGCGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCAGGACAGTCCACGGGGCGACATAAATCAGTTCCAGCACAACACCGAGCATCGCTCCACATGCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCTGCGATATCGAGAATTGCCATCTACAGTGGTCGGGTGAtc                                                                                                                                              <  2:427301/251‑1 (MQ=255)
                                                                                               atatCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGACCGCCATCGCAAAGGAGGACACCCAACAACACGCTAATCAAAATTGCCGCCCCTGGCAGCATCACAGTCCGCGGGGCGACATAAATCAGTTCCTGCGCAACCCCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGACCGGGAAAGCCCAGCGCGTCGATATCGTGATTTGCCATCGagaaggcgcggttgacc                                                                                                                                              >  2:189763/1‑235 (MQ=255)
                                                                                                          tGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGAcc                                                                                                                                     >  1:520206/1‑249 (MQ=255)
                                                                                                           ggCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGc                                                                                                                                  >  2:463161/1‑251 (MQ=255)
                                                                                                           ggCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGAGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAg                                                                                                                                   <  2:7326/250‑1 (MQ=255)
                                                                                                              attattCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTg                                                                                                                                        >  1:326990/1‑242 (MQ=255)
                                                                                                              attagtCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTg                                                                                                                                        <  2:326990/242‑1 (MQ=255)
                                                                                                                 attCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCcggcg                                                                                                                            >  2:67270/1‑251 (MQ=255)
                                                                                                                     cACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTg                                                                                                                         >  1:471231/1‑250 (MQ=255)
                                                                                                                          ccGCAATAATCGCACGGCGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGATAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGttt                                                                                                                   <  2:471686/251‑1 (MQ=255)
                                                                                                                           cGCNATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTg                                                                                                                  <  2:595711/251‑1 (MQ=255)
                                                                                                                            gCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTgg                                                                                                                 <  1:360142/251‑1 (MQ=255)
                                                                                                                              aataatCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGc                                                                                                                >  2:120453/1‑250 (MQ=255)
                                                                                                                               ataatCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCAt                                                                                                              <  2:117688/251‑1 (MQ=255)
                                                                                                                                     gCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGCTCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGc                                                                                                        >  2:179733/1‑251 (MQ=255)
                                                                                                                                     gCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGc                                                                                                        >  2:222850/1‑251 (MQ=255)
                                                                                                                                         gacgaACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAAc                                                                                                    >  2:504439/1‑251 (MQ=255)
                                                                                                                                         gacgaACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCAACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGCAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAAc                                                                                                    >  1:653247/1‑251 (MQ=255)
                                                                                                                                                  cATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAat                                                                                            >  2:238632/1‑250 (MQ=255)
                                                                                                                                                           ggaggTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGtt                                                                                  >  1:593875/1‑251 (MQ=255)
                                                                                                                                                            gaggTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTg                                                                                 <  1:222371/251‑1 (MQ=255)
                                                                                                                                                                ttAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATgc                                                                             >  1:414932/1‑251 (MQ=255)
                                                                                                                                                                           cacGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGc                                                                                                                                                                                             <  1:183171/128‑1 (MQ=255)
                                                                                                                                                                           cacGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGc                                                                                                                                                                                             >  2:183171/1‑128 (MQ=255)
                                                                                                                                                                              gCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCGCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACggg                                                               <  1:67270/251‑1 (MQ=255)
                                                                                                                                                                              gCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACggg                                                               >  1:58545/1‑251 (MQ=255)
                                                                                                                                                                                     aTAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGAt                                                        >  2:456941/1‑251 (MQ=255)
                                                                                                                                                                                     aTAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGAt                                                        <  2:575958/251‑1 (MQ=255)
                                                                                                                                                                                              gCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTActc                                               <  2:414932/251‑1 (MQ=255)
                                                                                                                                                                                                   tGGCAGCATGACAGTCCACGGGGCGACATCAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAg                                                                                                                                                           >  1:532108/1‑138 (MQ=255)
                                                                                                                                                                                                   tGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAg                                                                                                                                                           <  2:532108/138‑1 (MQ=255)
                                                                                                                                                                                                   tGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTTGCATCACAGCAAACCAGAGAATATTCAGCGt                                                                                   >  2:20232/1‑210 (MQ=255)
                                                                                                                                                                                                   tGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGt                                                                                   <  1:20232/210‑1 (MQ=255)
                                                                                                                                                                                                    ggCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTTCTCTTTTTc                                         <  1:308034/251‑1 (MQ=255)
                                                                                                                                                                                                      cagcaTGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGAGCCATCCCGACGGGCGGCGATAACGTACTCTTTTTccc                                       >  2:574382/1‑250 (MQ=255)
                                                                                                                                                                                                        gcaTGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGAGCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGt                                     <  2:109463/251‑1 (MQ=255)
                                                                                                                                                                                                           tGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCACCACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGAGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCCTCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGTTCg                                  >  2:526229/1‑251 (MQ=255)
                                                                                                                                                                                                                  ccACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGc                            <  1:120453/250‑1 (MQ=255)
                                                                                                                                                                                                                   cACGGGGCGACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCAc                          >  1:650460/1‑251 (MQ=255)
                                                                                                                                                                                                                           gACATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGt                  >  1:57279/1‑251 (MQ=255)
                                                                                                                                                                                                                            aCATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCCTGCTGt                  >  1:249860/1‑250 (MQ=255)
                                                                                                                                                                                                                            aCATAAATCAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTa                 <  2:640281/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cAGTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAAc         >  1:363783/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gTTCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGt       <  1:238632/251‑1 (MQ=255)
                                                                                                                                                                                                                                        tCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTAc     >  2:176248/1‑251 (MQ=255)
                                                                                                                                                                                                                                        tCCAGCGCAACACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGAGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTAc     >  1:360261/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aGCGCAACACCGACCATCGCTCACCAATCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGGTGAGCGCGGCGATATCGAGAATTTCCATCGACAGGGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTGGATGCGCCAGCCAGACGGGCGGCGATAACGTACTCTTTGTCCAGTTCGTCATGCACCAGGCTGTAAATCGAACGTACCAt  <  2:75718/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
GACCGCGGATTTCACCTTCGGTTAAGGTCATACTTACCCGGTCGACGGCTTTAACCCACTCATCACCGGTTTTAAATTCAATGGTCAGGTTACGAATATCGAGTAATGGCATTATTCCACCCCCGCAATAATCGCACGACGAACGCCATCGCCAAGGAGGTTAACCAACAACACGCTAATCATAATTGCCGCACCTGGCAGCATGACAGTCCACGGGGCGACATAAATCAGTTCCAGCGCATCACCGAGCATCGCTCCCCATTCAGGCGAGGGGAGTTGTGCGCCGAGATCGAGAAAGCCCAGCGCGGCGATATCGAGAATTGCCATCGACAGTGCGCGGGTGATCTCGGTGACCAGCCCGGCGGTGATGTTTGGCATCACAGCAAACCAGAGAATATTCAGCGTTGATGCGCCATCCAGACGGGCGGCGATAACGTACTCTTTTTCCAGTTCGTCATGCACCATGCTGTAAATCGAACGTACCAT  >  NC_000913/1353517‑1354002

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