Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,790,458 G→A 19.8% W60* (TGG→TGA uhpB → signal transduction histidine‑protein kinase/phosphatase UhpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077793,790,4580GA19.8% 160.9 / 34.5 75W60* (TGG→TGAuhpBsignal transduction histidine‑protein kinase/phosphatase UhpB
Reads supporting (aligned to +/- strand):  ref base G (9/51);  new base A (3/12);  total (12/63)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

AAAACTGGGCGTCAGTAACGACGTAGAGCTGGCGCGCCGCATGTTTGATGGCTGGTGATGAAGACGTTGTTCTCCCGCTTAATTACCGTTATTGCCTGCTTTTTTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGGCGGCGGCGTTGTTGCAGTCGCTGCCCTGGCTTTGGCACGGCAAAGAG  >  NC_007779/3790222‑3790686
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                    
aaaaCTGGGCGTCAGTAACGACGTAGAGCTGGCGCGCCGCATGTTTGATGGCTGGTGATGAAGACGTTGTTCTCCCGCTTAATTACCGTTATTGCCTGCTTTTTTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGc                                                                                                                                                                                                                                     <  1:1127681/238‑1 (MQ=255)
                aaCGACGTAGAGCTGGCGCGCCGCATGTTTGATGGCTGGTGATGAAGACGTTGTTCTCCCGCTTAATTACCGTTATTGCCTGCTTTTTTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGgcg                                                                                                                                                                                                                      <  1:405965/237‑1 (MQ=255)
                   gacgTAGAGGTGGCGCGCCGCATGTTTGATGGCTGGTGATGAAGACGTTGTTCTCCCGCTTAATTACCGTTATTGCCTGCTTTTTTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGACCCGTATTGCTGGGCGCGGAg                                                                                                                                                                                                                 <  1:1537367/239‑1 (MQ=255)
                        agagCTGGCGGGCCGCATGTTTGATGGCTGGTGATGAAGACGTTGTTCTCCCGCTTAATTACCGTTATTGCCTGCTTTTTTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGgc                                                                                                                                                                                                             <  2:1051183/238‑1 (MQ=255)
                           gCTGGCGCGCCGCATGTTTGATGGCTGGTGATGAAGACGTTGTTCTCCCGCTTAATTACCGTTATTGCCTGCTTTTTTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGACCCGTATTGCTGGGCGCGGAGTGgctgc                                                                                                                                                                                                          <  1:434260/238‑1 (MQ=255)
                                 gcgcCGCATGTTTGATGGCTGGTGATGAAGACGTTGTTCTCCCGCTTAATTACCGTTATTGCCTGCTTTTTTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGACCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTa                                                                                                                                                                                                   <  1:584583/239‑1 (MQ=255)
                                 gcgcCGCATGTTTGAGGGCTGGTGATGAAGACGTTGTTCTCCCGCTTAATTACCGTTAATGCCTGCTTTTTTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTa                                                                                                                                                                                                   <  1:1152397/239‑1 (MQ=255)
                                           tttGATGGCTGGTGATGAAGACGTTGTTCTCCCGCTTAATTACCGTTATTGCCTGCTTTTTTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGc                                                                                                                                                                                                                                     <  1:1211108/195‑1 (MQ=255)
                                                    tGGTGATGAAGACGTTGTTCTCCCGCTTAATTACCGTTATTGCCTGCTTTTTTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGcggt                                                                                                                                                                                 <  2:564398/238‑1 (MQ=255)
                                                    tGGTGATGAAGACGTTGTTCTCCCGCTTAATTACCGTTATTGCCTGCTTTTTTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGcggt                                                                                                                                                                                 <  2:1276983/238‑1 (MQ=255)
                                                             aGACGTTGTTCTCCCGCTTAATTACCGTTATTGCCAGCTTTTTTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAg                                                                                                                                                                                                                 <  1:1451529/197‑1 (MQ=255)
                                                                      tctcCCGCTTAATTACCGTTATTGCCTGCTTTTTTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACAGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCg                                                                                                                                                              <  2:351228/239‑1 (MQ=255)
                                                                          ccGCTTAATTACCGTTATTGCCTGCTTTTTTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTAt                                                                                                                                                          <  1:230498/239‑1 (MQ=255)
                                                                          ccGCTTAATTACCGTTATTGCCTGCTTTTGTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCGGTTATTTCCGTTTGGTCTGCGTCGGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCACGTTCCGTTAt                                                                                                                                                          <  2:587768/239‑1 (MQ=255)
                                                                                 atgcgcGTTAGTGCCTGCTTTTTTAGCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTgatgat                                                                                                                                                   <  2:11904/234‑1 (MQ=255)
                                                                                 aTTACCGTTATTGCCTGCTTTTTTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTgatgat                                                                                                                                                   <  1:205782/239‑1 (MQ=255)
                                                                                 aTTACCGTTATTGCCTGCTTTTTTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTgatgat                                                                                                                                                   <  1:801901/239‑1 (MQ=255)
                                                                                 aTTACCGTTATTGCCTGCTTTTTTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTgatgat                                                                                                                                                   <  2:1131471/239‑1 (MQ=255)
                                                                                      cGTTATTGCCTGCTTTTTTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTa                                                                                                                                              <  2:1661713/239‑1 (MQ=255)
                                                                                                                       aTGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGAtc                                                                                                              <  2:1248332/238‑1 (MQ=255)
                                                                                                                       aTGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGAtc                                                                                                              <  1:95423/238‑1 (MQ=255)
                                                                                                                             ttGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCt                                                                                                       <  1:1023874/239‑1 (MQ=255)
                                                                                                                                          aTCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGc                                                                                           <  2:864139/238‑1 (MQ=255)
                                                                                                                                                  gCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGgc                                                                                  <  1:1263950/239‑1 (MQ=255)
                                                                                                                                                    aTCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGgcg                                                                                 <  1:969241/238‑1 (MQ=255)
                                                                                                                                                    aTCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGgcg                                                                                 <  1:878388/238‑1 (MQ=255)
                                                                                                                                                    aTCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGgcg                                                                                 <  1:128976/238‑1 (MQ=255)
                                                                                                                                                    aTCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGACCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGgcg                                                                                 <  1:390452/238‑1 (MQ=255)
                                                                                                                                                    aTCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGACACGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGgcg                                                                                 <  2:1265815/238‑1 (MQ=255)
                                                                                                                                                     tCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCa                                                                               <  1:1412917/239‑1 (MQ=255)
                                                                                                                                                      cTGGTTGGGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCAc                                                                              <  2:702442/239‑1 (MQ=255)
                                                                                                                                                      cTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCAc                                                                              <  1:628773/239‑1 (MQ=255)
                                                                                                                                                      cTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCAc                                                                              <  1:30068/239‑1 (MQ=255)
                                                                                                                                                      cTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCAc                                                                              <  1:575857/239‑1 (MQ=255)
                                                                                                                                                          ttGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGACCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCtt                                                                           <  1:630877/238‑1 (MQ=255)
                                                                                                                                                          ttGACCGCCCTGGTAGGCCGGTGCTGTTATTTCGTTTTGGTTTGCGTTTGGGGCTAATGCTGCGATGCCCGGGCGGATACGGGGCCGTATTGCGGGGCGCGGAGTGGCTGCTGATTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTg                                                                          <  2:452036/238‑1 (MQ=255)
                                                                                                                                                           tGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTg                                                                          <  1:761384/238‑1 (MQ=255)
                                                                                                                                                           tGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTg                                                                          <  1:1096514/238‑1 (MQ=255)
                                                                                                                                                           tGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTg                                                                          <  1:142992/238‑1 (MQ=255)
                                                                                                                                                           tGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTg                                                                          <  1:279614/238‑1 (MQ=255)
                                                                                                                                                           tGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTg                                                                          <  1:250468/238‑1 (MQ=255)
                                                                                                                                                              gcgcCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTg                                                                       <  1:161581/238‑1 (MQ=255)
                                                                                                                                                                  ccTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGACCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATtgatg                                                                                                                                                     <  1:1235310/156‑1 (MQ=255)
                                                                                                                                                                        atGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTGATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAggggg                                                             <  2:1496602/238‑1 (MQ=255)
                                                                                                                                                                                    ttATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAgcg                                                                                                                    <  1:126406/171‑1 (MQ=255)
                                                                                                                                                                                    ttATTTCCGTGTGGTCTGCGTCTGGGGCTAATGCTGCGATGCCCGCGCGGATACTGAACCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAAGCGTTCGGTTTAACCCAGTGTCCGTTATTGATGATCGGTAGTTTCCTGACGTTACTGCCCGTAGCGCGGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAAcggc                                                <  2:778247/239‑1 (MQ=255)
                                                                                                                                                                                         tCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAAcggcggcgg                                           >  1:762890/1‑239 (MQ=255)
                                                                                                                                                                                         tCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGGCGGTTTTAACCATTTTCCGTTATTTATGGTCGGGAGATTACTGAAGTTTATGCCCGTCGCCCTGCTCTCGCGCCATCCCCCTCCGGGTGGCTGGCCCACCCTGCCTTTTCAAGGGGGGGCGCTAAcggcggggg                                           >  2:1301315/1‑239 (MQ=255)
                                                                                                                                                                                         tCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGACCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTTACCCCTTTTCCGTTATTGATTATCGGTAGGTTACTTACGGTACTTCCCGTAGCGCTGATCTCCCGCTATCGCCATCAGCGTGACTGGCGCACCTTTCTGTTTCACGGGGGGGGGTTAAcgggggggg                                           >  2:27185/1‑239 (MQ=255)
                                                                                                                                                                                         tCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGACCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAAcggcggcgg                                           >  1:719209/1‑239 (MQ=255)
                                                                                                                                                                                           cGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGACCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGgcggcggcg                                         <  1:1432176/239‑1 (MQ=255)
                                                                                                                                                                                              ttGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGGCGGCGGCgttg                                      >  1:1541722/1‑239 (MQ=255)
                                                                                                                                                                                                ggTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGACCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGGCGGCGGCgttgtt                                    >  1:496564/1‑239 (MQ=255)
                                                                                                                                                                                                 gTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTGACCCCTTTTCCGTTTTTGTTGATCGGTAGTTTACTCACGTTACTGCCCCCAGCGCTTATCTCTCCCTATCGCCCTCCGCCCGGCTGGCGCACCTCTCCGTTTCAGGGGGGGGGCTGAACCgcggcggcggggt                                     >  2:755656/1‑233 (MQ=255)
                                                                                                                                                                                                       gTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGGCGGCGGCGTTGTTGCAGTCg                             <  1:1449399/239‑1 (MQ=255)
                                                                                                                                                                                                         cTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTAc                                                                                                                                        <  1:526385/130‑1 (MQ=255)
                                                                                                                                                                                                              gCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCAc                                                                              <  1:1393736/183‑1 (MQ=255)
                                                                                                                                                                                                              gCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCt                                                                            <  1:1812375/185‑1 (MQ=255)
                                                                                                                                                                                                                tAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGGCGGCGGCGTTGTTGCAGTCGCTGCCCTgg                    >  1:1176926/1‑239 (MQ=255)
                                                                                                                                                                                                                   tgctgcAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGGTTTACCCCTTTTCCGGTATTTATGATCGGGAGGTTACTGACGTTACt                                                                                                                              >  2:526385/1‑130 (MQ=255)
                                                                                                                                                                                                                    gcttcAATGCCCTGGGGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGGCGGCGGCGTTGTTGCAGTCGCTGCCCTGGCtt                 <  2:1847009/238‑1 (MQ=255)
                                                                                                                                                                                                                       gcAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGGCGGCGGCCGTGTTGCAGTCGCTGGCCCGGGTTTgg              >  1:1284970/1‑238 (MQ=255)
                                                                                                                                                                                                                        cAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCAc                                                                              <  2:1100808/173‑1 (MQ=255)
                                                                                                                                                                                                                        cAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCAc                                                                              >  1:1100808/1‑173 (MQ=255)
                                                                                                                                                                                                                         aaTGCCCGGGCGGATACTGGCCCGTATTGCTGGGAGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTCCTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGGCGGCGGCGTTGTTGCAGTCGCTGCCCTGGCTTTGGCa            <  2:1130772/238‑1 (MQ=255)
                                                                                                                                                                                                                         aaTGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGGCGGCGGCGTTGTTGCAGTCGCTGCCCTGGCTTTGGCa            <  1:1210615/238‑1 (MQ=255)
                                                                                                                                                                                                                         aaTGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGGCGGCGGCGTTGTTGCAGTCGCTGCCCTGGCTTTGGCa            <  1:18612/238‑1 (MQ=255)
                                                                                                                                                                                                                         aaTGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGGCGGCGGCGTTGTTGCAGTCGCTGCCCTGGCTTTGGCa            <  1:304004/238‑1 (MQ=255)
                                                                                                                                                                                                                          aTGCCCGCGCGGATACTGACCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGGCGGCGGCGTTGTTGCAGTCGCTGCCCTGGCTTTGGCACg          <  2:496564/239‑1 (MQ=255)
                                                                                                                                                                                                                              ccGCGCGGATTCTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTGCTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGGCGGCGGCGTTGTTGCAGTCGCTGCCCTGGCTTTGGCACGGCa       <  2:1797333/238‑1 (MQ=255)
                                                                                                                                                                                                                              ccGCGCGGATACTGACCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGGCGGCGGCGTTGTTGCAGTCGCTGCCCTGGCTTTGGCACGGCa       <  1:1583462/238‑1 (MQ=255)
                                                                                                                                                                                                                               cgcgcgGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGGCGGCGGCGTTGTTGCAGTCGCTGCCCTGGCTTTGGCACGGCaaa     <  1:1508151/239‑1 (MQ=255)
                                                                                                                                                                                                                               cgcgcgGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGGCGGCCGCGGTGTTGCAGGCGCTTCCCCGGCCTTTGCACGGCaaa     >  1:1505441/1‑239 (MQ=255)
                                                                                                                                                                                                                               cgcgcgGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGAGCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGGCGGCGGCGTTGTTGCAGTCGCTGCCCTGGCTTTGGCACGGCaaa     <  1:17097/239‑1 (MQ=255)
                                                                                                                                                                                                                                  gcgGATACTGACCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGGCGGCGGCGTTGTTGCAGTCGCTGCCCTGGCTTTGGCACGGCAAagag  <  1:803302/239‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                    
AAAACTGGGCGTCAGTAACGACGTAGAGCTGGCGCGCCGCATGTTTGATGGCTGGTGATGAAGACGTTGTTCTCCCGCTTAATTACCGTTATTGCCTGCTTTTTTATCTTCTCTGCCGCATGGTTTTGCCTGTGGAGTATCAGCCTGCATCTGGTTGAGCGCCCTGATATGGCGGTGCTGTTATTTCCGTTTGGTCTGCGTCTGGGGCTAATGCTGCAATGCCCGCGCGGATACTGGCCCGTATTGCTGGGCGCGGAGTGGCTGCTGATTTACTGGCTAACGCAGGCGGTCGGTTTAACCCATTTTCCGTTATTGATGATCGGTAGTTTACTGACGTTACTGCCCGTAGCGCTGATCTCGCGCTATCGCCATCAGCGTGACTGGCGCACCTTGCTGTTACAGGGGGCGGCGTTAACGGCGGCGGCGTTGTTGCAGTCGCTGCCCTGGCTTTGGCACGGCAAAGAG  >  NC_007779/3790222‑3790686

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