Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,571,567 T→C 8.6% S14G (AGC→GGC)  U712_18320 ← putative allantoin permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,571,5670TC 8.6% 252.6 / 11.4 93S14G (AGC→GGC) U712_18320putative allantoin permease
Reads supporting (aligned to +/- strand):  ref base T (42/43);  new base C (6/2);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGAAAGGCCGATCGCAATTAATCCTCCCACGGTCGCGTAGGTCGGTATGTTATGTATACATCCCATCCAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGTTTTTTGACATTAGATGTCACTTAATATAACTGAAGAGTCAGACAATAACA  >  CP006881/3571429‑3571703
                                                                                                                                           |                                                                                                                                        
gCGAAAGGCCGATCGCAATTAATCCTCCCACGGTCGCGTAGGTCGGTATGTTATGTATACATCCCATCCAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCACAGGCACCAGATCTTCATTGCt                                                                                                                                          <  7:69260/139‑1 (MQ=255)
gCGAAAGGCCGATCGCAATTAATCCTCCCACGGTCGCGTAGGTCGGTATGTTATGTATACATCCCATCCAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCt                                                                                                                                          <  5:225151/139‑1 (MQ=255)
                    aaTCCTCCCACGGTCGCGTAGGTCGGTATGTTATGTATACATCCCATCCAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATtgctgc                                                                                                                       <  2:73859/138‑1 (MQ=255)
                                              tatgttatgtATACATCCCATCCAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGc                                                                                            >  6:222984/1‑139 (MQ=255)
                                                  ttatgtATACATCCCATCCAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCACAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTcc                                                                                         <  8:53413/138‑1 (MQ=255)
                                                  ttatgtATACATCCCATCCAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTcc                                                                                         <  1:92590/138‑1 (MQ=255)
                                                  ttatgtATACAACCCATCCAAATGGAGGCAAAGTGCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTcc                                                                                         <  5:27493/138‑1 (MQ=255)
                                                   tatgtATACATCCCATCCAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTcccc                                                                                       >  6:146110/1‑139 (MQ=255)
                                                   tatgtATACATCCCATCCAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCACAAAGGCACCACATCTTCATTGATCAGCCTGTTGGCTTTCTGCTGACTCGCTTTTAATTTCATAGACGGCTccc                                                                                        >  6:16876/1‑138 (MQ=255)
                                                   tatgtATACATCCCATCCAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTACCAAAGGCACCAGATCTTCATTGCTCATCCTGGTGGTTTGCTGCTGACTCTCTTTTAATTTCATAGCAGGCTccc                                                                                        >  6:27493/1‑138 (MQ=255)
                                                   tatgtATACATCCCATCCAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTACTGTCCCAAAGGCAACAGATCTTCATTGCTCAGCCTGTATGATTTCTGCTGACTCTATTTTAATTTCATATCCGGCTacc                                                                                        >  6:204087/1‑138 (MQ=255)
                                                   tatgtATACATCCCATCCAAATGGAGGCAAAGTTCAATGCTTTCCACGTCCGTTTTTGCTGTACCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTGACTTTTAATTTCATAGCAGGCTcccc                                                                                       >  2:144436/1‑139 (MQ=255)
                                                       tataCATCCCATCCAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAATGCACCAGATCTTAATTTCCCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGTCTCCCCTTTc                                                                                   >  5:152545/1‑139 (MQ=255)
                                                       tataCATCCCATCCAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCttt                                                                                    <  4:304350/138‑1 (MQ=255)
                                                       tataCATCCCATCCAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTc                                                                                   >  8:299583/1‑139 (MQ=255)
                                                        ataCATCCCATCCAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTc                                                                                   >  8:43418/1‑138 (MQ=255)
                                                        ataCATCCCATCCAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGAGGTTCATTGCCCCGCCTGTTGGCTTGCTGCTGACGATCGTTTAATTTCATATCCGGCTccccagtcc                                                                                  >  6:4429/1‑134 (MQ=255)
                                                        ataCATCACATCCAAATCGACGAAAAGATAATTGCTTTCCACGTCCGTTTTTCCTGTACCAAAGTAACCAAAACTACAAACCTCAAACAGTTGGATTGCTGCAGCCTCACTTTAAATTTCATAGCCGGCTCCCGTTTc                                                                                   >  4:296868/1‑138 (MQ=255)
                                                          aCATCCCATCCAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCa                                                                                 >  6:80363/1‑138 (MQ=255)
                                                           catcccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGTATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCata                                                                               <  7:308292/139‑1 (MQ=255)
                                                           catcccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCata                                                                               <  5:79996/139‑1 (MQ=255)
                                                           catcccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCata                                                                               <  1:144436/139‑1 (MQ=255)
                                                           catcccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACATGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCata                                                                               <  7:164967/139‑1 (MQ=255)
                                                           catcccatccAAATGGAGGCAAAGGGCATTGCTTTCCACGTCCGTTTATCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCata                                                                               <  2:95026/139‑1 (MQ=255)
                                                           catcccaaccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCAAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCata                                                                               <  2:270695/139‑1 (MQ=255)
                                                            atcccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCt                                                                                           >  4:94248/1‑126 (MQ=255)
                                                            atcccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCt                                                                                           <  3:94248/126‑1 (MQ=255)
                                                            atcccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCatat                                                                              >  3:98607/1‑139 (MQ=255)
                                                            atcccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCatat                                                                              >  2:18997/1‑139 (MQ=255)
                                                            atcccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCatat                                                                              >  3:22981/1‑139 (MQ=255)
                                                            atcccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCGGTTGGATTGCTGCTGACTCTCTTTTTATTTAATAGCCGGCTCCCCGTTCCatat                                                                              >  4:7701/1‑139 (MQ=255)
                                                            atcccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCCCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCata                                                                               <  2:191617/138‑1 (MQ=255)
                                                            atcccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCCCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCata                                                                               >  1:191617/1‑138 (MQ=255)
                                                            atcccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCGTCATTGCTCAGCCTGTTGGATTGCTGCTGACTAGCGTTTAATTTCATAGCAGGAGCCCCTTTCCatat                                                                              >  8:103604/1‑139 (MQ=255)
                                                            ataccatccAAATGGAGGCAAAGTTAATTGCTTTCCACGTCCGTGTTTCCTGTCACAAAGGCAACAGATATTCATTGCCGAGCCTGATGGATTGGTGCTGGCTCTATTTTTATTTCATAGACGGCTCCCATTTCCatat                                                                              >  4:115120/4‑139 (MQ=255)
                                                              cccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCata                                                                               <  7:266465/136‑1 (MQ=255)
                                                              cccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCata                                                                               <  3:154486/136‑1 (MQ=255)
                                                              cccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCata                                                                               >  4:154486/1‑136 (MQ=255)
                                                              cccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCata                                                                               >  8:266465/1‑136 (MQ=255)
                                                              cccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCCCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTg                                                                            >  8:171610/1‑139 (MQ=255)
                                                              cccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCACAAAGGCACCAGATCTTCATTGATCAGCCTGTTGGATTGCTTCTGGCTCTCTTTTATTTTCATTGCAGGCTACACTTTCCATAttt                                                                            >  8:164967/1‑138 (MQ=255)
                                                               ccatccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGt                                                                           >  8:88383/1‑139 (MQ=255)
                                                               ccatccAAATGGAGGCAAAGTTCATTGCTATCCACGTCCGTGTTTACATGACCCAAGGGCACCCGAAGTTAATTGGTCCGCCGGTTGGTTTGCTGGTGCGGTGGTTTTATTTTCATCGCCGGACCACCTTGCAATAttt                                                                            >  2:20155/1‑138 (MQ=255)
                                                                 atccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCa                                                                                                                                        <  8:202066/76‑1 (MQ=255)
                                                                 atccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTTCTGCTGACTCTCTTTTAATTTCATATCCGGCACCCCGTTCCAAATTTTc                                                                          >  6:210579/1‑138 (MQ=255)
                                                                 atccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTcc                                                                                         <  4:25550/123‑1 (MQ=255)
                                                                 atccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTcc                                                                                         <  3:23627/123‑1 (MQ=255)
                                                                 atccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTcc                                                                                         >  3:25550/1‑123 (MQ=255)
                                                                 atccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTcc                                                                                         >  4:23627/1‑123 (MQ=255)
                                                                 atccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTc                                                                          >  6:288125/1‑138 (MQ=255)
                                                                 atccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCATCATTGCTCAGCCTGTTGGATTGCTTCTGACTCTATTTTAATTTCATAGCCGGCTCCCCTTTCCAAAATGTc                                                                          >  8:308292/1‑138 (MQ=255)
                                                                 atccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCACATCTTCATTGCTCAGCCAGTTGGATTGCTGCTGACTCTCTTTTAATTTCATATCCGGCGCCCCGTTACATAtttact                                                                         >  1:245119/1‑135 (MQ=255)
                                                                 atccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACAAGATCTTCATTGACCAGCCTGTTGGATTGATGCTGACTCTCTTTTAATTTCATATACGGCTCCCCTTTCCATAATGTc                                                                          >  4:85829/1‑138 (MQ=255)
                                                                 atccAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTGTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTTATGCGGACTCTGAATTAATTTCATAGCAGGCTCCCCGATGCTTATTGTc                                                                          >  5:60144/1‑138 (MQ=255)
                                                                 atccAAATGGAGGCAAAGTTCATTGCTTTCCACGACCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCa                                                                                                                                        >  7:202066/1‑76 (MQ=255)
                                                                      aaTGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGCTCTTCATTGCTCAGCCTGGTGGCTTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCGGTCCAGATTGTGTATATg                                                                    >  5:276524/1‑139 (MQ=255)
                                                                             gCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCCCCAGAACTTCATTGCTCAGCCTGTAGGATTGCTGCTGACTCTCTTTTAATTTCATATCCGGCTCCCCATTCAAAATTGTCTATATGTCGTca                                                              >  6:185962/1‑138 (MQ=255)
                                                                                 ttttCATTGCTTTCAAAGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGGTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCGCTTTCCATATTGTCTATATGTCGTCACACGt                                                         <  3:296868/137‑1 (MQ=255)
                                                                                   ttCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGAct                                                                                                                  >  7:233437/1‑80 (MQ=255)
                                                                                   ttCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGAct                                                                                                                  >  5:266455/1‑80 (MQ=255)
                                                                                   ttCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGAct                                                                                                                  <  6:266455/80‑1 (MQ=255)
                                                                                   ttCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGAct                                                                                                                  <  8:233437/80‑1 (MQ=255)
                                                                                     cATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCata                                                                               <  1:116465/113‑1 (MQ=255)
                                                                                     cATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCata                                                                               >  2:116465/1‑113 (MQ=255)
                                                                                          cTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCGCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGttttt                                                <  3:53987/139‑1 (MQ=255)
                                                                                          cTTTCCACGTAGTTTATTCCTGTCCCAAAGGCACAAGATCATCATTGTTAAGCCGGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACTCGTTATGttttt                                                <  2:103749/139‑1 (MQ=255)
                                                                                            ttCCACGTCCGTTTTTCCTGTCGCAAAGGCACCAGATCTTCATTGCTAAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCAGGCTCCCCGTTCCATATAGTCTATATGTAGTCAAACCTTATGTTTTTTg                                              >  4:107563/1‑139 (MQ=255)
                                                                                             tCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATgtcgt                                                                >  7:87466/1‑120 (MQ=255)
                                                                                             tCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATgtcgt                                                                <  8:87466/120‑1 (MQ=255)
                                                                                             tCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGtttttt                                               <  5:16876/137‑1 (MQ=255)
                                                                                             tCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGGTATGTTTTTTg                                              >  1:164227/1‑138 (MQ=255)
                                                                                                      gTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGTTTTTTGACATTAGATg                                    <  5:80363/139‑1 (MQ=255)
                                                                                                           ttCCTGTCCCAAAGGCTCCAGATCTTCATTGCTCCGCCTGTTGGATTGGTGCTGACTCTCTTTTAATTTCATAGCAGGCTCCCCTTTCCATATTGTCTATATGTCGTCACTCGTTATGTTTTTTGACATTAGATGTCa                                 <  2:245119/138‑1 (MQ=255)
                                                                                                           ttCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGTTTTTTGACATTAGATGTCa                                 <  5:288125/138‑1 (MQ=255)
                                                                                                           ttCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGTTTTTTGACATTAGATGTCa                                 <  5:255444/138‑1 (MQ=255)
                                                                                                           ttCCTGTCCCAAAGCCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCTTAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGTTTTTTGACATTAGATGTCa                                 <  3:107563/138‑1 (MQ=255)
                                                                                                           ttCCTGGCCCAACGGCACCAGATCTTCATTGCCCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGTTTTTTGACATTAGATGTCa                                 <  3:85829/138‑1 (MQ=255)
                                                                                                                 tCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGTTTTTTGACATTAGATGTCACTTaa                            <  5:222984/137‑1 (MQ=255)
                                                                                                                           aCCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGTTTTTTGACATTAGATGTCACTTaa                            >  2:244098/1‑127 (MQ=255)
                                                                                                                           aCCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGTTTTTTGACATTAGATGTCACTTaa                            <  1:244098/127‑1 (MQ=255)
                                                                                                                             cAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGTTTTTTGACATTAGATGTCACTTAATATAACTGAAGAGt              <  7:299583/139‑1 (MQ=255)
                                                                                                                             cAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGTTTTTTGACATTAGATGTCACTTAATATAACTGAAGAGt              <  5:210579/139‑1 (MQ=255)
                                                                                                                                aTCTTCTTTGCTCAGCCTGTTGGATTGCTGCTGAAGCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGTTTTTTGACATTAGATGTCACTTAATATAACTGAAGAGTCa            <  6:126135/138‑1 (MQ=255)
                                                                                                                                   ttCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTcccc                                                                                       <  6:90701/60‑1 (MQ=255)
                                                                                                                                   ttCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTcccc                                                                                       >  5:90701/1‑60 (MQ=255)
                                                                                                                                     cATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCta                                                                        <  3:214242/73‑1 (MQ=255)
                                                                                                                                     cATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCta                                                                        >  4:214242/1‑73 (MQ=255)
                                                                                                                                     cATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGTTTTTTGACATTAGATGTCACTTAATATAACTGAAGAGTCa            <  2:164227/133‑1 (MQ=255)
                                                                                                                                       ttGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGTTTTTTGACATTAGATGTCACTTAATATAACTGAAGAGTCAGACAATaa    >  2:118672/1‑139 (MQ=255)
                                                                                                                                          cTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTcc                                                                                         <  8:306181/51‑1 (MQ=255)
                                                                                                                                          cTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTcc                                                                                         >  7:306181/1‑51 (MQ=255)
                                                                                                                                          cTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGTTTTTTGACATTAGATGTCACTTAATATAACTGAAGAGTCAGACAATAaca  <  5:146110/138‑1 (MQ=255)
                                                                                                                                          cTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGTTTTTTGACATTAGATGTCACTTAATATAACTGAAGAGTCAGACAATAaca  <  3:7701/138‑1 (MQ=255)
                                                                                                                                          cTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGTTTTTTGACATTAGATGTCACTTAATATAACTGAAGAGTCAGACAATAaca  <  1:18997/138‑1 (MQ=255)
                                                                                                                                          cTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGTTTTTTGACATTAGATGTCACTTAATATAACTGAAGAGTCAGACAATAaca  <  4:98607/138‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                        
GCGAAAGGCCGATCGCAATTAATCCTCCCACGGTCGCGTAGGTCGGTATGTTATGTATACATCCCATCCAAATGGAGGCAAAGTTCATTGCTTTCCACGTCCGTTTTTCCTGTCCCAAAGGCACCAGATCTTCATTGCTCAGCCTGTTGGATTGCTGCTGACTCTCTTTTAATTTCATAGCCGGCTCCCCTTTCCATATTGTCTATATGTCGTCACACGTTATGTTTTTTGACATTAGATGTCACTTAATATAACTGAAGAGTCAGACAATAACA  >  CP006881/3571429‑3571703

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