Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,321,603 T→C I32V (ATC→GTC)  phnG ← carbon‑phosphorus lyase core complex subunit PhnG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,321,6030TC98.0% 167.5 / NA 51I32V (ATC→GTC) phnGcarbon‑phosphorus lyase core complex subunit PhnG
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (18/32);  minor base A (1/0);  total (19/32)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGGCGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGCGCACGCACAGCAGCGGCGACATATTCGGGATTTCAAGACGCTGGCACTCTTTGGCCTGGGCGCGGCGGGCGCTGATCCGCGTCTGGCTGCGGC  >  NC_000913/4321313‑4321884
                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                         
ggCAATGCGTGGCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAgcg                                                                                                                                                                                                                                                                                   <  2:84183/299‑1 (MQ=255)
ggCAATGCGTGCCATATGATCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAGGTGGCGAGATTGCTGCAGCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTGATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGGCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAgcg                                                                                                                                                                                                                                                                                   <  2:84840/299‑1 (MQ=255)
 gCAATGCGTGCCATACGGTCAGCGTCTAGCTGTGCAATAATGGTTTCTGATAAGTTTTGAAGGTGGCTAGATTGTTGCATCAGCGCGTCAGTCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGTCGAGTGTGCTTTCAGTCAGGCGCACTGCGACGCGGGTCAGCGTGGCGTCGCCGTCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCATGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAgcg                                                                                                                                                                                                                                                                                   <  2:84322/298‑1 (MQ=255)
 gCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAATGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCATGCCAGTTTCAGCGGCGCGGATCACCTCATATTCGGCGGTGACGTTTAgcg                                                                                                                                                                                                                                                                                   <  2:172672/298‑1 (MQ=255)
        gTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCCGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGc                                                                                                                                                                                                                                                                          >  2:583352/1‑300 (MQ=255)
        gTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTATGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCCCCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTcc                                                                                                                                                                                                                                                                             >  2:610959/1‑296 (MQ=255)
                     aGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGtt                                                                                                                                                                                                                                                              <  1:269623/299‑1 (MQ=255)
                     aGCGGCCAGCTGGGCAGTAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTTGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGTCGTTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTc                                                                                                                                                                                                                                                             <  2:352687/300‑1 (MQ=255)
                      gCGTCCAGCGGGGCAATAGGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCCCTCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGNCAGTCCGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGACGGCAAAAAAACGTTCGCAGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGAAGTTCAGCGAGTTCAGGCGTGCTGCAAGTTc                                                                                                                                                                                                                                                             >  2:507707/1‑299 (MQ=255)
                      gCGTCCAGCGGGGCAATAAGGGTTTCTGGTAAGTTTTGAAAGTGGCGAGATTGCTGCATCGGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGNCAGTCAGACGCACGGCGACGCGGGTCAGCGTGGCGTCGCAGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGACAGTTTAAGCGGCGCGGATTACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCCGtt                                                                                                                                                                                                                                                              >  2:507526/1‑298 (MQ=255)
                      gCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTCCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTAGGCGTGTTGCCAGTTc                                                                                                                                                                                                                                                             >  2:424797/1‑298 (MQ=255)
                      gCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTAATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCa                                                                                                                                                                                                                                                            >  1:166819/1‑300 (MQ=255)
                                 ggCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTAATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCt                                                                                                                                                                                                                                                 >  1:71872/1‑300 (MQ=255)
                                                          tGAAAGTGGCTAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGc                                                                                                                                                                                                                        <  2:147680/300‑1 (MQ=255)
                                                          tGAAAGTGGCGAGATTGCTGCAGCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGc                                                                                                                                                                                                                        <  1:696091/300‑1 (MQ=255)
                                                               gTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCg                                                                                                                                                                                                                   <  1:416832/300‑1 (MQ=255)
                                                                 gtcGAGATTGCTGCAGCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTCATCACGCCCCTGCACCCAGCTGTAGCTGAGCGTGCCTTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGCTCGCCGGTGCCGCCCATCCGCGCCTAAACCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCGTAGTCGGCGGTGACGTTTAgcg                                                                                                                                                                                                                                                                                   <  2:172704/232‑1 (MQ=255)
                                                                                           aTCAGCGCGCAGCGTTCTGCGTACTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGTCGACGCGGGTCAGCGTGGCGTCGCCTGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGc                                                                                                                                                                                       <  2:191966/300‑1 (MQ=255)
                                                                                                        gTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAAt                                                                                                                                                                           <  2:137440/299‑1 (MQ=255)
                                                                                                                    gTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTAATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTAAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTc                                                                                                                                                                >  2:579238/1‑298 (MQ=255)
                                                                                                                              ccccTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTc                                                                                                                                                      >  1:462391/1‑298 (MQ=255)
                                                                                                                                          gCTGTAGCCGAGCGTGCTGTCAGTCAGGCGCACGGCGACGCGGGTCAGTGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGTCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGGCGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCCGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCCGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTa                                                                                                                                        >  2:177497/1‑300 (MQ=255)
                                                                                                                                          gCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTa                                                                                                                                        >  1:5043/1‑300 (MQ=255)
                                                                                                                                            tGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTAc                                                                                                                                       <  1:478732/299‑1 (MQ=255)
                                                                                                                                                               agtcagGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGNATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCAcc                                                                                                                   <  1:511371/300‑1 (MQ=255)
                                                                                                                                                               agtcagGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCAcc                                                                                                                   <  2:686534/300‑1 (MQ=255)
                                                                                                                                                                              gacgcgGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAg                                                                                                    <  1:617154/298‑1 (MQ=255)
                                                                                                                                                                                 gcgGGTCATCGTTGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAggg                                                                                                  <  2:433514/299‑1 (MQ=255)
                                                                                                                                                                                 gcgGGTCAGCGTTGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCGGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCAGAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAggg                                                                                                  <  2:624942/299‑1 (MQ=255)
                                                                                                                                                                                 gcgGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAggg                                                                                                  <  1:457704/299‑1 (MQ=255)
                                                                                                                                                                                   gggTCAGCTTGGCGTCGCCGGCAAAAAAACGTTCGCTAGTGCCACCCATCCGTGCCTAAATCTGCACCCTGCCGTTTTCGTCGTCGCGGATCCACTAATGGTCTGCGGTAACGTTTAGCGCTTTCAGGCGTGCTTCCAGTTCGGCCGTTTTGCTGTGCGCCAGCACGGACATCCAGTGCTGGCAGGTCGCGGTATCTGTGTGTATTCGGTGTTCTATGGTGAATTCAGTTATTTCGGAGCGCGTCGGGCTGACGGAGTATTCCGCTGGGTTGTTTTCACTGTCACGGTGGTTAAggg                                                                                                  <  2:447898/297‑1 (MQ=255)
                                                                                                                                                                                      tccgcttggtgtcgccttcaaaaggaCGTGCGCCGGTGCCGCCCATTCGCGCCTGAATCTGTACCATTCCAGTTTCATCGGCGCGGATCACCTCATCTTCTGCGTTGACATTTAGCGCGTTCAGGCGTGCTTCCAGTTCAGCCGTTTGGTTGTGCGCCAGCACGGACATCCAGTTCTTGCGGGTCGCAGTACCTGCGTGCAGTCAGTGCTCTATAGTGAATTCGATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGcgca                                                                                            <  2:88576/275‑1 (MQ=255)
                                                                                                                                                                                      tcagctttgCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCTTTAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCTGCGGTGACGTTTAGCGCGTTCAGGCGTGCGGCCAGTTTAGCCGGTTGGCTGTGCTCCAGCACGGACATCCAGTGCTGGCGGTTCGCAGTATCTGCGTGCATTCGGTGCTCCATGGTGAATTCAATCATGTCTGCGCGCTTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCATCACGGTGGTGAAGGGTGcgca                                                                                            <  2:429780/292‑1 (MQ=255)
                                                                                                                                                                                      tCAGCTTGGCGTCGCCGGCAAAAAAACATTCGCCGGTGCCGCCCATCCGCACCTGAATCTGTACCAGGCCAGTTGCAGCGGCGCGGATCACCTCATAGTCGGCGTTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCCGTTGGCTGTGCGCCAGCACGGACATCTAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAACTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAggg                                                                                                  <  2:293623/294‑1 (MQ=255)
                                                                                                                                                                                      tCAGCGTGGCGTCGCCTGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACATTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATATTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGcgca                                                                                            <  2:429772/300‑1 (MQ=255)
                                                                                                                                                                                         gCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCATTCAGGCGTGTTGCCAGTTCAGCCGGTTGGCTGTGCTCCAGCCCGGACATCCAGTGCTGGCGGGTCGCGGTCTCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCAGGGCTGCTTTCACCGTCAAGGTGGTTAAGGGTGcgcccgc                                                                                         >  2:369145/1‑300 (MQ=255)
                                                                                                                                                                                           gtgtcGTCGCCGTCAAGAAAACTTTCGCCGGTTCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTTACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGTGTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAATGGTGCgcacgcac                                                                                       <  2:627560/296‑1 (MQ=255)
                                                                                                                                                                                                   gCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGCGCACGCACAGCAgcgg                                                                               <  1:527815/300‑1 (MQ=255)
                                                                                                                                                                                                                         gTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCAGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGCGCACGCACAGCAGCGGCGACATATTCGGGATTTCacga                                                         >  2:395803/1‑297 (MQ=255)
                                                                                                                                                                                                                         gTGCCGCCCATCCGCACCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAATGCTCCATAGTGAATTCAATCATGTCGGCGCACGTCAGGCTGACGGAGCACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTACGGGTGCGCACGCACAGAAGCGGCGACATATTCGGGATTTCAAGa                                                         >  2:395798/1‑300 (MQ=255)
                                                                                                                                                                                                                              gccCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGCGCACGCACAGCAGCGGCGACATATTCGGGATTTCAAGACGCt                                                     >  1:638230/1‑299 (MQ=255)
                                                                                                                                                                                                                                  aTCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGCGCACGCACAGCAGCGGCGACATATTCTGGATTTCAAGACGCTGGCAc                                                <  1:458435/300‑1 (MQ=255)
                                                                                                                                                                                                                                  aTCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGCGCACGCACAGCAGCGGCGACATATTCGGGATTTCAAGACGCTGGCAc                                                <  1:573150/300‑1 (MQ=255)
                                                                                                                                                                                                                                           tGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATGTTCGGTGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGTTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGCGCACGCACAGCAGCGGCGACATATTCGGGATTTCAAGACGCTGGCACTCTTTGGc                                        <  2:638230/299‑1 (MQ=255)
                                                                                                                                                                                                                                            gAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGTCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGCGCACGCACAGCAGCGGCGACATATTCGGGATTTCAAGACGCTGGCACTCTTTGGCCt                                      <  1:417043/300‑1 (MQ=255)
                                                                                                                                                                                                                                             aaTCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGCGCACGCACAGCAGCGGCGACATATTCGGGATTTCAAGACGCTGGCACTCTTTGGCCt                                      >  1:111688/1‑299 (MQ=255)
                                                                                                                                                                                                                                             aaTCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGAAGGAGTACTCCGCAGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGCGCACGCACAGCAGCGGCGACATATTCGGGATTTCAAGACGCTGGCACTCTTTGGCCt                                      >  2:305116/1‑299 (MQ=255)
                                                                                                                                                                                                                                                 tGTATCAGGCCAGTTTCAGCGGCGCGGATCATCGCATAGTCTGCGTGGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCAAGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGCGCACGCACAGCAGCGGCGACATATTCTGGATTTCAAGACGCTGGCACTCTTTGGCCTGGgc                                  <  2:138616/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gcgcGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGCGCACGCACAGCAGCGGCGACATATTCGGGATTTCAAGACGCTGGCACTCTTTGGCCTGGGCGCGGCGGGCGCTGATCCGCGTc            >  2:614794/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                          ggATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCGGTGCGCCAGNACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGCGCACGCACAGCAGCGGCGACATATTCGGGATTTCAAGACGCGGGCACTCTTTGGCCTGGgcctgtctctgatacacatctcgcag          <  1:141609/298‑25 (MQ=255)
                                                                                                                                                                                                                                                                                ccTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGCGCACGCACAGCAGCGGCGACATATTCGGGATTTCAAGACGCTGGCACTCTTTGGCCTTGGCGCGGCGGGCGCTGATCCGCGTCTGGCTgcggc  >  1:560273/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                         
GGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGGCGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGCGCACGCACAGCAGCGGCGACATATTCGGGATTTCAAGACGCTGGCACTCTTTGGCCTGGGCGCGGCGGGCGCTGATCCGCGTCTGGCTGCGGC  >  NC_000913/4321313‑4321884

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